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1.0  INTRODUCTION 

The  flutter  problem  can  be  solved  with  either  a  collocation  or 
normal-mode  formulation.  The  collocation  approach  is  attractive  if  an 
accurate  stiffness  and  aerodynamic  influence  coefficient  matrix  can  be 
formed  for  the  system.  The  normal-mode  method  merits  consideration 
when  mode  shapes  and  natural  frequencies  for  the  structure  are  known. 
The  normal-mode  method  requires  the  aerodynamics  to  be  presented  as 
generalized  forces.  The  normal-mode  formulation  is  generally  presented 
in  the  following  manner 


Where 


=  Generalized  Mass  Matrix 


=  Generalized  Aerodynamic  Force  Matrix 
«  Artifical  Structural  Damping 


=  Flutter  Frequency 
-  Natural  Vibration  Frequency 


=  Generalized1  Coordinate 


This  volume  presents  three  computer  programs  that  calculate  the 
generalized  aerodynamic  forces  for  the  three  flight  regimes:  Subsonic 
Flight,  Transonic  Flighty  Supersonic  Flight  and  may  be  used  with  the 
Modal  Flutter  Analysis  Program  of  Vol.  III.  The  subsonic  program  is 
based  upon  Kernel  Function  formulation.  The  transonic  and  supersonic 
programs  are  based  upon  the  Mach  Box  technique.  (  , 


1 


2.0  UNSTEADY  AERODYNAMIC  GENERALIZED  FORCES 


The  generalized  force  Q 


ij 


is  defined  as  the  work  done  in  mode  i 
ii 


by  the  pressures  due  to  motion  of  mode  j.  If  z  (x,y)  is  the  .deformation 
pattern  of  mode  i  and  A?.  (x,y)  is  the  pressure  distribution  of  mode  j 
then  ■* 


z^Cx^y)  APj (x,y)dxdy 


2.0.1 


Oscillation  of  a  lifting  surface  creates  a  phase  lead  (or  lag)  between 
the  pressure  and  motion  of  the  surface  thereby  making  Q. ,  a  complex 


quantity. 


ij 


The  frequency  dependency  of  the  generalized  force  makes  the 
flutter  solution  a  trial  and  error  procedure.  Different  frequencies 
must  be  tried  until  the  flutter  frequency  and  the  generalized  force 
frequency  coincides.  Fortunately,  the  flutter  frequencies  are  usually 
near  the  natural  frequencies.,  This  is  especially  true  for  a  stable 
structure  with  modes  which  remain  fairly  uncoupled;  there  is  little 
drifting  of  oscillatory  frequencies  with  the  presence  of  aerodynamic 
forces. 

Various  lifting  surface  theories  were  employed  to  determine 
the  pressure  distribution  APj(x,y)  for  generalized  force  calculations. 
This  was  required  to  correctly  account  for  the  high  degree  of  chordwise 
and  spanwise  deformation  associated  with  mode  shapes  of  low  aspect  ratio 
surfaces.  Also,  it  was  necessary  to  include  aerodynamic  interaction 
effects  between  tandem  surfaces. 


The  methods  are  based  oh  the  linearized  equation  of  fluid  motion 
which  describe  the  flow  patterns  for  a  compressible,  inviscld-,  isentropic 
and  irrotational  fluid.  The  boundary  conditions  are  consistent  with  thin 
wing  theory.  The  effects  of  high  angles  of  attack  and  missile  body 
influence  on  the  aerodynamics  of  the  lifting  surfaces  were  not  considered 
in  the  analysis. 

The  linearized  equation  for  time  dependent  disturbances  is 


*XX  +  ^YY  +  ^ZZ  M  ^XX  +  U  ^XT  +  U 


2.0.2 


where  $  is  the  perturbation  velocity  potential  at  the  point  (X,  Y,  Z) . 

M  and  U  are  the  free  stream  Mach  number  and  velocity,  respectively.  Impos¬ 
ing  the  requirement  of  harmonic  motion,  this  time  dependency  can  be  expressed 


where  $  is  the  nondimens ional  complex  amplitude, of  the  velocity  potential 
and  b  is  a  reference  sem±-chord.  The  equation  may  be  cast  in  a  nondimerisior.al 
form  by  substitution  of  the  following  dimensionless  parameters: 


x  =  X/b 
y  =  Y/b ' 
z  =  Z/b 
k  =  aib/U 


This  yields 


<f>  + 

xx 


yy. 


$ ,  =  M 

zz 


+  2ik 


XX 


2.0.3 


It  is  this  relationship  which  must  be  satisfied  for  M<1,  M=  1,  and  M >T.  A  brief 
description  of  the  analytical  techniques  employed  for  each  speed  regime  is 
presented  in  the  following  sections.  The  numerical  and  computational  schemes 
used  are  similar  to  those  which  are  described  in  detail  iV.  References  5,  6,  and  7<. 


3.0  SUBSONIC  UNSTEADY  AERODYNAMICS  PROGRAM 


3.1  Theoretical  Derivation 

Subsonic  aerodynamic  loads  were  derived,  by  the  kernel  function 
method  and  the  resulting  loads  were  then  used  to  compute  generalized 
forces.  The  subsonic  kernel  function  relationship  for  the  dovnwafch 
(x,y)  on  a  surface  in  terms  of  the  pressure  over  the  entire  surface  Is 


w(x,y)  «  -  -  ff— - 2  K<x-S>  y* 

8  JJ  1/2PU 


7))dSdT) 


3.1.1 


This  relationship  is  deriyed  by  noting  that  the  potential  equation  for 
subsonic  flow. 


(l-M2)*  +  <*>  +0  =  M2  (2ikd>  -  k Z<t>) 

xx  - yy  ^zz  x 

is  satisfied  at  ,(x,y,z)  by  the  pulsating  pressure  doublet 

f  =  -A  _L,-i Hx-Or^fr  exp  (X-HR')j  <1 

where  ■ 

R'  =  ^fX2  +  (1-M2)  (y-7j)2  +  (z-O2 
and  the  pressure  doublet  strength  is  ;given  by 


A  4  V 


1/2 Pii 


3.1.2 


3.1.3 


3.1.4 


3.1.5 


The  velocity  normal  to  the  lifting  surface  is  given  by 
w(x,y)  =  lim 


3a. 6 


Substituting  equation  (3.1.3)  into  the  equation  above  and  integrating  over 
the  area  yields  the  kernel  function  relationship,  equation  (3.1;1).  For 
tandem  surfaces,  the  integration  extends  over  both  surfaces: 


w(x,y)  =  -  p. 


.  +s  t.e. 

•fe  /  / 


1/2PU 


cs  (g,n) 


K(x-£,  y-r,)d£dn 


-s  l.e. 

control  surface 


1/2 


K(x-£,  y-n)  d§dn 


t 


The  kernel  function  at  any  point  is  given  by 


K{x-S,y-n)  -  £  £ 


.  x-4  = 

jy  i.  f  ’  _x ,  i  —  o-hr'):J 

dz  as  J  S'  1-m2  ax 

„ .  ^CO  L  -J  ; 


(3.  .1.8)  r. 


Equation  (3.1.7)  then  represents  the  integral  equation  wherein  given  the  downwash 

v(x,y)  over  the  wing  and  control  surface,  AP  and  AP  must  then  be  determined. 

W  '  cs 

The  pressure  distribution  is  approximated  as  the  sum  of  a  series  of 
functions  which  have  the  proper  behavior  as  inferred  from  steady  state  and 
two-dimensional  solutions.  The  pressure  on  each  surface  can  be  approximated 
in  the  form 


where 


,  _  J.2  2  V  M 

b(n)  n=  o  iu=o 


fo  0 


VP* 

V-l  +  £ 


fn  C5)  -  C  Tr_i(C)'3  n>l 


(3.1.9) 


Po  (n)  -  1.0 


PTn  =  **  1  (h)  •)  m>  1 

m  —  —  m— X  —  — 


(3.1.10) 


T  (x)  =  1.0 
o 


Tx(x)  =  2x 


Tk  ^  =  2xTk-l^^  “  \-2^  ’  k  -  2 


5  -  % 
b  (n) 


«  ’  1/2  <  «1...  +  5 


s  is  the  starboard  coordinate  of  the  surfaice  tip  and  b  '07)  is  the  local 
semi-chord.  The  functions  T  are  Chebyshev  polynomials  and  are  introduced 
for  purposes  of  convenience. 


Substituting  equation  (3.1.8,  3.1.9,  and  3.1.10)  into  equation  (3.1.7) 
yields  a  set  of  equations  relating  the  pressure  coefficients  for  the  wing, 

W  c  s 

a  ,  and  the  control  surface,  a  .  to  the  dcwnwashes;  In  matrix  hotat ion 
run  nm 

this  gives 


(3.3.11) 


The  integrals  involved  in  evaluating,  the  D’s  were  solved  by  the  methods  of  Hsu 
(Reference  5).  In  this  procedure  the  integrals  are  numerically  evaluated  using 
the  Gauss~Mehler  quadrature.  Upon  determining  values  for  the  D's,  the  pressure 
coefficients  are  found  by  direct  inversion  and  the  generalized  forces  are  found 
-by  application  of  equation  (2.0.1). 


3.2  Program  Description. 

The  Subsonic  Generalized  Force  Unsteady  Aerodynamics  Program  calculates 
generalized  forces  for  up  to  10  deformation  modes.  The  computer  solution, 
which  is  based  upon  the  kernel  function  formulation,  is  applicable  to  a- 
variety  of  configurations.  The  various  configurations  which  can.be  ana¬ 
lyzed  are  shown  in  Figure  3.2.1  and  Table  3.2.2.  The  analysis  includes 
interaction  effects  between  tandem  surfaces  and  wake  effects  on  the  trailing 
surface.  The  number  of  integration'  stations  are  chosen,,  and  they  are  auto¬ 
matically  located.  The  collocation  stations1  are  then  interdigitated  between  the 
integration  stations,  according  to  Hsu  (Reference  5),  internally  in  the  program. 
The  solutions  at  the  collocation  stations  are  then  matched  to  terms  in  the 
downwash  series  by  a  least  squares  method  and  the  surface  pressure  are  deter¬ 
mined.  The  method  of  solution  programmed  does  not  allow  for  a  single  surface 
to  be  analyzed  separately;  however,  an  option  to  isolate  and  eliminate  interfer¬ 
ence  effects  between  the  two  surfaces  is  available.  Thus,  a  single  surface 
can  be  analyzed  by  using  the  option  IS0LAT  and  inputing  a  dummy  second 
surface. 


The  solution  for  the  generalized  aerodynamic  forces  requires  the 
input  of  the  deformation  modes  due  to  vibration.  The  program  considers 
the  modes  to  be  expressed'  as  analytic  functions  of  the  form: 

K  n 

w.(x,  y)  =  £  £ 
in— U  ir-0 


'(n-m). 


m 


(n-m)  m 

x  y 


To  meet  this  requirement  only  the.  coefficients  "c"  are  required  as  input 
into  the  program.  These  coefficients  can  be  obtained  in  several  ways, 
the  most  common  way  is  to  surface  fit  the  modes  by  the  least-square 
technique. 
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TABLE  3.2.2  -  OPTIONAL  -CONFXGlUlATlCNS 


V't 

4  •: 


■PT 

12 


J 


II  j 

~  j 


if  J 

i  i 

w  a 

n. 


~  j 
fli 
rh 

-*■  J 

*  / 

Hi 

Hi 


Configuration 

,  Chordyiso  Coordinates 

,  Spanwise  Coordinates 

Rectangular 

X(l)  =0.0 

X(2)  =0.0 
;  x(3)  >  o.o 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y(l)  =  0.0 

Y(2)  =  0.0 

Y(3)  >  0.0 

Delta 

‘  X(l)  -  0.0 

X(2)  >  0.0 

X(3)  =  X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y(l)  'm  0.0 

Y(2)  =  0.0 

Y(3)  >  0.0 

Trapezoidal 

. 

X(l)  -  0.0 

X(2)  >0.0 

X(3)  =  X(2) 

X(4).  >X(3) 

X(5)  >  X(4) 

Y(l)  =  G.O 

Y(2)  >  0.0 

Y(3)  >Y(2) 

Trapezoidal 
• (Cropped) 

X(l)  =  0.0 

X(2)  >  X(l) 

X(3)  >  X(2) 

*(4)  >  X(3) 

X(5)  >  X(4) 

Y(l)  =  0.0 

Y(2)  >0.0 

Y.(3)  >Y(2) 

.Delta 

(Cropped) 

X(l)  -  0.0 

X(2)  >  0.0 

X(3)  >X(2) 

X(4)  >  X(3) 

X(5)  >  X(4) 

Y(l)  =0.0 

Y(2)  =  0.0 

Yt(3)  >  Y(2) 

3.3  TSLinr  INSTRUCTIONS 


Instructions  for  preparing  input  data  £>r  the  subsonic  computer  program 
arc  presented  here.  The  field  location  and  format  fox  each  input  quantity  is 
specified.  Any  set  of  units  may  be  used  for  geometric  dimensions  and  acoustic 
velocity  as  long  as  they  are  consistent  c.g..,  if  feet  is  used  for  length  then 
the  acoustic  velocity  must  have  dimensions  of  feet  per  second. 

1.  Strcamwise  Coordinates  (6E12.5  format) 


Column 

Kane 

Item 


(1)  X(] )  King  root  leading  edge  streanwise  coordinate  (See  Figure  3.3.1) 

(2)  X(2)  Wing  tip  leading  cage  streamsisc  coordinate 

(3)  X(3)  Wing  trailing  edge  streamwise  coordinate 

(4)  X(4)  Control  surface  leading  edge  streamaise  coordinate 

(5)  X(5)  Control  surface  trailing  edge  streamoise  coordinate 

The  technique  for  generating  various  configurations  is  shove  in  Table  3.2 

The  origin  for  the  pianfom  and  AIC  station  coordinates  mast  be  at  the 
leading  edge  root  of  the  uing  therefore  X(l)  and  Y(l),  described  bel ea, 
must  always  be  0.0. 

2.  Spanvjise  Coordinates  and  Acoustic  Velocity  (6E12.S  format) 


Column 

Name 

Item 


(1)  Y(l)  Wing  root  spanwise  coordinate 

(2)  Y(2)  Wing  leading  edge  spanv7ise  coordinate 

(3)  Y(3)  Wing  (and  control  surface)  tip  spanwise  coordinate 

(4)  SOUND  Acoustic  velocity  for  altitude  at  which  analysis  is  performed 

(5)  RHO  density  of  fluid  -at  1000.0  (M/L^) 


1-12  : 

13-24 

25,-36 

37r4S 

49-60 

61-72 

Y(l) 

Y(2) 

Y(3) 

SOUND 

RHO 

(1) 

cai _ 

(3) 

(5) 

1-12 

13-24 

25-36 

37-4S 

49-60 

1  61-72  * 

X(!) 

X(2)„ 

X(3). 

X(4) 

X(5)  . 

JD _ 

(3) 

_ 

_j5 ) _ ; 

1 

_ _ l 

10 


General  Information  (6112  format) 


Column 

Saas 

Item, 


13-24 


1  25-35  !  37-48 

49-60  i 

1  61-721 

1  XH&BES  !  XC0LL 

1  o)  :  i  m 

IPRifSH 

LPR a 

(6) 

(1)  NXACH  Number  or  Kach  numbers  (max  6) 

(2)  NFREQ  Number  of  frequencies  at  each  Kach  cumber  (max  1G) 

(3)  NM0DES  Number  of  deformation  nodes  (max  10) 

(4)  LC0LL  Print  collocation  stations;  1  — Yes;  0  '-'No 

(5)  LPRWSH  Print  pressures  and  downvashes;  1  '-—Yes;  0  — '  No 

(6)  LPRC0  Print  pressure  coefficients;  1  --  Yes;  0  —  No 

4.  General  Information  (6112  format) 


Column 

Same 

Item 


■  1-12  4  13-24  !  25-36 

-  NVCX  i  NWPX  I  NCCX 

(1)  I  (2)  !  (3) 


37-48  r  49-60  1  61-72 

»■■■.■  _  I. 

NCPX  NIgNCX  8 _ 

(4)  i  (5)  1 (6) 


(1)  NNCX  Number  of  chord  collocation  stations  -  wing  (max  10) 

(2)  NWPX  Number  of  chord  pressure  series  terms-  -  wing  (max  10) 

(3)  NCCX  Number  of  chord  collocation  stations  -  control  surface (max  10) 

(4)  NCPX  Number  of  chord  pressure  series  terms  -  control  surface (max  1G) 

(5)  NI0NCX  Choose  a  value  of  NI0NCX  such  that  Jfl^NCX  *(S»CX  or  NCCX) 

equals  the  number  of  cnordwise  integration  stations 
/„  40  40  ] 

H^KHCX  °r  NCCX/ 

5.  General  Information  (6112  format) 


Column  |_1 

Name  i_S 

Item  l_J 

(1)  NIY 

(2)  NWCY 

(3)  N?Y 

(4)  INWTS 

(5)  IS0LAT 


1-12 
S  NIY. 

I  (1) 


13-24 

25-36 

37-48 

49-60 

61-72  ! 

NWCY 

NPY 

INWTS 

ISdIAT 

- 3 

_£2 1_ 

_ (3) 

(4) 

_ (3) 

-j 

Number  of  spaiiwisa  integration  stations  (max  11) 

Number  of  spanwise  collocation  stations  -  wing  (max  11) 
Number  of  spanwise  pressure  series  terms  (max  10) 

Read  downwasb  error  weighting  factors;  l'-'  Yes;  0— -  No 
Isolate  wing  and  control  surface;  l'^-'Yes;  0-“'So 


Weighting  Factors  (6E12.5  format) (omit  these  cards  when  INWIS  =  0) 


Column 

1-12 

13-24  I 

25-36 

Name 

WT(1) 

WT(2) 

WT(3) 

Item 

(1) 

(2)  _ 

1  43) 

Continue  on  successive  cards  until  WT(i)  =  WT(KWTS) 
Where  NWTS  =  NWCY  *  NWCX  +  NIY  *  NCCX 


7.  X0HTt  (6112  format)  (omit  these  cards  when  NWCY  =  NIY) 

Cobras 
Kame 
Item 

Continue  on  successive  cards  until  X0H(i)  =  K0M(N0MIX)  where  K0MIT  * 
HIY-KWCY.  For  the  solution  to  be  carried  out,  KWCY  must  equal  NIY. 
When  an  excessive  number  of  collocation  stations  exist-,  they  must  be 
eliminated  by  spanwise  rows  K0M(i)  is  the  spanwise  row  number  to  be 
eliminated .  KjSffiT  =  NIY-KWCY  is  the  number  of  spanwise  rows  to  be 
eliminated. 


(1)  FKACE(l)  Mach  number 

(2)  FKACH(2)  Mach  number 


(HMACH)  FK&CH(KHACH)  Mach  number 

XMACH  values  of  Mach  number  must  be  input;  Mach  numbers  must  be 
greater  than  zero  and  less  than  0.95. 

9.  Frequencies-  (6E12. 5  format) 

Column 
Name 
Item 

(1)  FREQ(l)  f(cps) 

(2)  FREQ (2)  f(cps) 


(NFREQ)  FREO (NFREQ) 

For  NFREQ  6,  continue  input  on  new  card. 


13 


Repeat  the  following  cards  i  =  1,.  .  . ,  NM0DES 

t 

10.  Number  of  Deformation  Mode  Polynomial  Coefficients  to-be  Read  for 
the  i£l2  Mode  First  Surface 

Column 
Name 
Item 

(1)  NZC0(l,i)  Number  of  polynomial  coefficients  to  be  read  for  the 
first  surface,  the  mode. 


Column 
Name 
Item 

O'0(i,  j)  Deformation  polynomial  coefficients  for  the  first  surface 
in  the  order:  QjO;  1,0;  0,1;  2,0;  1,1;  2,0;  etc. 
where  the  first  integer  is  the  power  of  "x"  and  the  second 
is  the  power  of  "y" 


Column 
Name 
Item 

(1)  NZC0(z,i)  Number  of  polynomial  coefficients  to  be  read  for  the 
second  surface  and  the  i-^1  mode. 


Format (6112) 


1-12 

13-24, 

NZC0(2,i) 

(1)  -  (2) 

Format  (6E12.5) 


1-12 

13-24..  „ 

25-36 

37-48 

49-60 

61-72  | 

C0(O,O) 

C0C1.O)  . 

.  C0'O.l) 

C0(2,O) 

.  c0a.i) 

-  -  1 

C0CO,2)i 

_ 

-  (2) 

.  (3> _ 

_ _ 

_ (5) 

_ C§LJ 

Format (61 12)  , 


1-12 

13-24 

NZC0Cl,i) 

_ 

_ £2) 

Column 

Name 

Item 


C0(i,j)  Deformation  polynomial  coefficients  for  the  Second  surface 
in  the  order:  0;0;  1,0;  0,1;  2,0;  1,1;  2,0;  etc.;  where 
the  first  integer  is  the  power  of  "x"  and  the  second  is 
the  power  of  "y:‘ 


1-12 

.  13-24 

25-36  . 

37-48 

49-60  j 

61-72 

C0 (0,0) 

C0(1,O)  l 

C0CO.D  - 

C0(2.O) 

C0C1.1)  1 

C0(O,2) 

All  . 

_ (21_ _ L 

_(3) . - 

:  (4) 

1 

_J5) _ 

_C§) _ 
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‘.f  I  ft 


3.4  'SAMPLE  PROBLEM 

The  generalized  forces  ere  calculated  f.or  the  configuration  below. 

The  flight  parameters  and  pertinent  input  data  are  presented  on  the  first 
page- of  the  computer  print  out.  , 

The  coefficients  of  the  deformation  modes  for  the  forward  suface  are  shown 
On  the  second  page  of  the  computer  print  out,  and  for  the  aft  surface  oh  the 
third  page  of  the  computer  print  out. 

0.0  c 


Ml 


:  ! 


53.500 


57.55" 


;  f 


63.55’ 


Figure  3.4.1 
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s 

s 

s 

CHAIN 


UPIIUH  FORTRAN 
FORTRAN  HLSTUU, DECK 
1NC0DE  18HF 
MAIN 

COMPLEX  A, AA , ANH, CZERU, GFGRC, DELP* HASH, APR, BHASH*CMASH,P* 
DIMENSION  DWASH(90,lQ),PR(9i,ll), CRASH ( 99 , li  ) 

DIMENSION  GF0RC(lfl»10,3),  OELPdC  }i» NASH ( It ) 

CQMHGN/C1/ A ( 60 ) >AA ( 50 , 60 ) , AMM( 59 , 11  ) * CZERO 


SUB7S02S 

SUB7B32B 

SU£7»33« 

SUF7I349 

SUB79358 

SUB7B3M 

SUB79378 

SUB703SB 


C0HH0N/C2/CLEN,  SN,  N I Y, NWCX,  NCCX, NHPX,  NCPX*HKER(21  )»2KEtC29)*N6SKlNSUB7839« 
CUHH0N/C3/NPY, SOUND, NKACH,FHACH(6),NFREU,FREO(li),HAUC,MO|ICX,*HO  SUB79491 
Cl)MM0N/C4/NH0nES» LCGLL, LPRHSH, LPRCQ,MOHl5},IIY,I IX* NSURF* ISOLAT  SUB7«4l> 
COMMON/C5/FW, FC* NCOLS*  NtiH IT*  ALPHA (18  )*  It (59) ,HC0R(6), ZC0R(6)*HACH  SUB7042* 


COMMON/ C6/ HXCMN ( 11 )  ,  HMCNdl  ) »  WBI N  (11  )'* HT (98  )  *  XCOLL,  YCOLL*Pl *U 


SUB7B43B 


C0MMQN/C7/C0( 10  *28, 2 X, NZCO (19  >  2 ) *  EM* EK* N2* NHI X*NCIX* W30*CB0M*MHCY  5UB7844I 


C 

4 


KOUT(2) 


7 


CUMH0N/C8/IFR,XE(5),YE(3>*UXUI),UYUi),HXIHN(ll)*El,E2 
C0HH0K/C9/S I X ( 4  0,2) , SCX (10 , 2 ) . Y (11  )  *  ETA  dll *  8MY, QHHX, QHCX*CXHM 
COMMON/ CPR/APk (90,60),  IHOD,  IRON 
EUli  1  VALENCE  (  OF  ORC,  A  A  ) ,  C  A,  HASH  )  ,(  HASH*  DfcLP  ) 

EgU  I  VALENCE  ( DWASH  ( 1 . 1 ) ,  APR  ( 1 , 51)  ) ,  (  PR  (1 ,1 )  ,  APRd,ll)  ) 

euu i Valence  (cuash(i,i ),apr<i.4D) 

UHI I £ ( 6*  66  ) 

FORMAT  LlHl ) 

CALL  KFBA 

Fk  =  2*NH1X  ♦  1 

FC  s  2*NC I X  ♦  1 

QHHX  =  2«0*PI/FH 

QWCX  =  2.a*PI/FC 

QHY  -  PI/FL0AT(2*NIY) 

CALL  GEOh 

DO  100  MACH=1;,  NMACH 
MACh  =  MACH 
EM  =  FMACH(MACH) 

CALL  KOUT(l) 
l)  s  EM*SOUND 

IF ( LCOLL. HE. 0 )  CALL  KOUT(2) 

BUU  =  HBO/U 
B2  =  1.0  -  EM«EM 

nu  loo  ifr=i,nfreu 
IF R  =  IFR 

EK  -  2.0*PI»FREUUFR)*BOU 
NSURF  =  1  (WING)  OR  2  (CONTROL) 

NCX  =  NUCX 
MOM  IT  =  1 

HAUG  =  NCOLS  *  NHUDES 

DO  4>  J=l* NCOLS 

IL(J)  ?  0 

DU  4  K=l* MAUG 

AA ( J, K )  =  CZERO 

I  ROM  =.  1 

DO  15  HSURF =1 , 2 

NSURF  =  NSURF 

CALL  K0UT(6> 

KOUT  (6)  PRINTS  COEFFICIENTS  OF  DEFLECTION  SERIES 
DO  14  I Y*l*  N I Y 

I  I  Y  =  IY 

IF(NOMIT-MOMIT.LT.O)  go  TO  7 
I F (I Y-NOM (MUM IT ) . EOvO )  GO  TO  13 
YCOLL  =  SN*  Y  C I Y  > 

DU  12  I X=1 , NCX 

I I  X  =  IX 


COEFFICIENTS  OF  DEFLECTION  SERIES 


^US7945I 
SUB7B46C 
SUB78479 
SUB7S48S 
SUB7B49S 
SUB785BI 
SU378518 
SUB7I52I 
SUB7853I 
SUB/9549 
SUB7B559 
SUB7856B 
SUB7I57B 
SUB7B580 
SU378598 
SUB7I69S 
SUB79618 
£087*626 
SUB7I638 
SUB78648 
SU87I658 
3U878668 
SU879673 
SUS7968B 
SUB73698 
SUB7I798 
SUB7371B 
SUB79729 
SUB7B730 
SUB78749 
SU870758 
SUB7I760 
SUB78770 
SUS7078B 
SUB70790 
SUS79890 
SUB70810 
SUB7Q820 
SUB 7 0830 
SUB70840 
SUB70850 
SUB/0860 
SUB/ 0870 
SUB70883 
SUB70890 
SUB7090U 


&  .  w  issi  -KKa  iwi  isk 


XCOLL  =  XS(1#MSURF* IX* I Y) 

CALL  CORD 

•«•  CCR0  FILLS  OUT  A  ROM  OF  THE  BCWMWASh  MATRIX  EACH  TIME  CALLED 
00  5*  M  *  lrXHOOES 
CALL  ZBZXtK, SLOPE* DEFL) 

HMC  =  MCOLS  *!i 

5«  HASHCHMCT  ~  CrtPLX(5L0PE*DEFL«EK/WR0> 

SO  41  K=1,MAUG 
XI  =  K 

IFCX.GT.MCOLS)  XI  =  5*  ♦  X  -  MCOLS 
4«  APRC IROW*Kl )-  *  A(X) 

DO  61  11=1* MCOLS 
6t  C AIL  CGREDCA.HsH) 

DO  7«  L=l?>MMODES 
ML  =  MCOLS  +  L 

7*  AL?3IA(L>  =  SQHTC ALPHA CLT**2  +  CABSC  A  (ML-)  >**2) 

JHOH  =  1KOH  *  1 

12  COMIIMUE 

CO  TO  14 

13  MOM IT  =  HUHIT  *1 

14  CtiKi I HUE 
MCX  =  NCCX 

15  CONI I HUE 
CALL  XLSO 


IFCLPRCC.EQ.8i  CO  TO  17 


00  16  IHOOsJ^UMOnES 
{MOD  »  1MO0 

16  CALL  XGUTC3) 

C  K0U1C3S  PRINTS  COEFFICIENTS  OF  PRESSURE  SERIES  FOR  EACH  MODE 

c  *********************,**************•**»************•**••**** 

17  IFfEPRHSH.Ea.M)  GO  TO  69 

C  ***************** ******************************************* 

c  COMPUTE  AND  SCORE  THE  PRESSURE  DOHNHASHES 
MUM1T  =  1 
IKON  =  l' 

MCX  =  MHCX 
00  115  MSURF  =1 *2 
DO  114  LY=£*HIT  . 

IECKOMIT.LT.MOHIT)  GO  TO  197 
IF C I Y.EO.NOHCHOMI D )  GO  TO  113 

107  DO  112  IX=1*NCX 

DO  140  IH=1,NH0DES 
AUM)  =  CZERC 
DO  140  JC=1, MCOLS 

140  A (  I.K J  =  AYIM)  +APRC IROR,  JC)*ARM(  JC,  IM) 

DO  150  I H=l»  NHODES 

150  CHASHCIROH,  IH)  =  A.CIH) 

IriOk  =  f ROM  +1 

112  CONTINUE 
GO  TO  114 

113  HOHIT  =  HOMIT  *1 

114  C0N1IHUE 
HCX  =  NCCX 

115  CONTINUE 


_  I*  " 

*  "  69  CALL  FORC  ( NR  I X*  NC  1'X*  N I Y  *  ETA*  SIX*WBIN»G) 

*  IF CLPRHSH.EO.O )  GS  TO  77 

CALL  F0RC( NWCX* NCCX, N I Y»  Y, SCX*  RBCN*  LPRHSH ) 

J  DO  9,9  l  M0I)=1 »  NHODES  ,  23 


SUB799I0  $ 
SU370920  | 

SOB7993D  | 
SUB70940  | 

SUB7  095(9 
SUB78969  fj 
SUB7897B  f" 
SUB70980  I 
SUB70990 
SUB/1000 
SUB71C10 
SUB71020 
SUB71030  9 

SUB/1040  H 
SOB  71 058 
SUB  74.060 


SUB71070 
SUB 71 080 
SUB71U90 
SUB 711 00 
SUB/lllO 
SUB/1120 
SUB/1130 
SUB71140 
►SU871150 
SUB7li60 
►SUB 7 1170 
suB7iiW 

SUB7119U 
SUB712Q0 
SUB/1210 
►SUD71220 
SUB71230 
►S0871240 
SUB71250 
SUB71260 
SUB71270 
SUB71280 
SU871290  . 

SUB7J300 
SUB71310 
SUB71320 
SU871330 
SUS71340 
SUB/1350 
SUB71360 
SUB71370 
SU871380 
SUB71390 
SUB/1400 
SUB/1410 
SUB71420 
SUB71430 
SUB71440 
SU8/1450 
SUB/1460 
SUB 7 1470  , 
SUB71480 
SUB71490 
SU871500 
SUB71510 


PRESSURE  AMS  BCk'NifASH  ARRAYS 


C 

c 

c 

c 

c 

c 


I  MOD  =  THOD 
NCX  =  NWCX 
HONlT  =  1 
I  HOW  =  1 
no  99  NSUF/F=1,2 
NSURF  =  NSURF 

CALL  KflUT  (7)  i 

KOUl  (7)  PRINTS  fHE  HEADER  FOR  THE 
9n  no  v8  iY-iiNir 

IFtNOMIT.LTiMOHlT)  GO  TO  91 
IF { I Y.FU»NOM{HOMI T) )  GO  TO  97 
93  YCOLL  =  W80*SN*Y ( I Y ) 

DO  96  lX=l,NCX 

XCOLL  =  WBO»XS(l,NSURF, IX, If) 

CALL  KOUT (9 ) 

K U II 1  (9)  PRINTS  THE  X  AND  Y  COORDINATES  OF  EACH  COLLOCATION 
POINT,  THE  LOCAL  PRESSURE,  THE’  DOMNWASH  CREATES  THERE  BY  THE 
PRESSURE  FIELD,  AND  THE  DOHNUASH  OF  THE  SURFACE  AT  THE  POINT* 

THE  DEGREE  TO  WHICH  THE  TWO  SETS  OF  BOHNWASHES  HATCH  IS  A  MEASURE 

of  i he  Accuracy  of  the  solution  of  ihe  houndary  value  problem. 

THE  ROOT- MEAN- SUUARE  OF  THE  ERRORS  IS  S1VEL  IN  THE  ARRAY  ALPHA. 
IROW  =  IROW  *1 
GO  TO  98  ' 

HON  IT  =  NON  IT  +1 
CONTINUE 
NCX  =  NCCX 
CONTINUE 


96 

97 

98 

99 

77 

78 
100 


00  78  L=l,3 
NSURF  =  L 
CALL  KOUT (8 ) 
KOIIT  (8)  PRINTS 
CONTINUE 
GO  10  1 
END 

FORTRAN 

1NCUDE 


THE  GENERALIZED  FORCES 


SU37152S 
SUB71538 
SUB7154* 
SUB7155* 
SUB7156C 
SUB7157I 
SUB7156S 
SUB/159* 
SUB71688 
SUB7161B 
SU87162C 
SUB7163* 
SUB7164S 
5UB71658 
SUB71668 
SUB71678 
SUB71658 
SUB7169* 
SUB717C8 
SUH71718 
SUB71728 
SU871738 
SUB7174C 
SUB71758 
SUB7176S 
SUR71778 
SUB71788 
*SU87179§ 
SUB7189S 
SUB71818 
SUB71828 
SUB71838 
SUB71848 
SUB71858 
SUB71868 
SUB71838 
SUB 7 8 *4 8 
SUB79B5I 
SUB7886B 
SU870878 
SUB7G0SS 


f  FORTRAN  NLSTOU, DECK 

%  1NCUDE  1 BHF 

CC0NS1  C0NS1 

HLOCK  DATA 

COMPLEX  A,  AA ,  ANM,  CZERll 

C0HM0N/C1/A(  6II  Y,  AA  C  50  ,  60  )  ,  ANH(50 ,10  >,  CZERO 
COMMON/ C2/CLEN , SN, N 1 Y, NWCX, NCCX, NWPX, NGPX, HKER (20 ) , ZKERC2B ),NGSKRNSUB70898 
COMMON/ C3/NPY, SOUND, NMACH, FMACHC6) , NFREu, FREQtlfl ) ,MAUG, NIOHCX, RHO  SUB?BI80 
COMMON/ C 4/ NMODES, LCOLL, LPHWSHi LPRCO, N0ML5 ) , 1 1 Y, 1 1 Xi N5URF , I SOLA!  SUB/0118 
C0MH0N/C5/FW,FC,NC0LS,NdHIT,ALPHA{lU),  II.  (50),HCOR{6),ZCOR(6),HACH  SIIB78J.2D 
ClJMMON/C6/WXCHN(ll),WHCNUl),WBIN(tl),WIC90),XCOLL,YCULL,PI,U  SUB78130 

COMMON/C7/CO(lO,28,2),NZCO(10,2),EM,EK,H2,NWIX,NCIX,WBO,CBON,NkCY  SUB/0140 
CUMM0N/C8/  I FR,  XE( 5 )  ,  YE (3.)  ,  UX (10  ) »  UY (10  ) , WX 1 MN (11 ) >E1» E2  SUB7B150 

C0MMflN/C9/SlX(4U,2),SCX<l0,2),Y(ll),ErA(ll)>QWY, QWWX, QWCX, CXMN  SUS70160 
nATA  P 1/3 • 14159265/  5U870170 

RATA  HCOR/O. 08566225, 0.18038079, 0.23395697, 0.23395697, 0.18038079,  SUB70180 
10.08566225/  SUB7019U 

DATA  ZC OR/ Q. 03376524 , 0. 16939531, 0. 38069 p 41, 0.61930 959, 0.830604 69,  SU870200 
10.96623476/  SUB70210 

(]  *#*«*  ####«*####  #*#***#******##**#**##*#**#»##***  **#♦#.*  #***#**#***»SUB  7  0228 

C  NUSKRN  SHOUl  D  BE  COMPATIBLE  WITH  HKER  A*D  ZKE.R  LISTS  SUB70230 

DATA  NGSKRN/8/  SUB70240 

DATA  ( HKER  (  I  ) ,  1=1 , 8  )/ 0 . U5U61427, 0 . 11 119i.i52, 0 .15685332, 0 .18134189  SUB/0250 

X, 0.18134189,0.15685332,0.11119052,0.05061427/  SUB70260 

DATA  (ZKER(  IJ,  l  =1 , 8  )/  0 . 01985507 , 0 .1  01t>6ft76 , 0 .23723380 , 0 .40828268  SU870270 

X, 0.59171732, 0.76276620, 0.89833324,0.98014493/  24  SUB70280 


k  hi 

<■  a 


*  f 

f  I 


s 

s 

CKFDA 


>|| 


►SUB7Q290 

DATA  El/8. 8100881/, E2/8«009Q861/,CZERO/(0. 0,0.5)/  SUB70300 

END  '  SUB 7 0310 

FORTRAN  NLSTOU,DECK  SU 8718 7.0 

IHCODE  I8MF  SU87186G 

KFDA  SUB71890 

SUBROUTINE  KFDA  SU871900 

COMPLEX  A,AA,ANM,CZERO  *  SU871910 

C0KHUM/C1/A(6U  )  , AA(50 ,68 ) , AhH(5C  ,18),CZER0  SU871920 

COMKOK/C2/CLEM,SN,NLY,NWCXiNCCX,NMPX,NCPX,HKER(20>,ZKER(20),N6SKRNSUB71930 
C0MH0H/C3/NPY, SOUND, NHACH, FMACH ( 6 ) , NFREu, FREQ ( 10 ),HAUG»MIONCX»  RHO  5UB71940 

C0HMGN/C4/NH0n£S, LCOLL, LPRHSH,LPRC0,N0H(5),  ILY,  I IX*NSURF,  ISOLAT  SU871950 


CUKHGN/C4/NHOnES,LCOLL,LFRHSH,LPRCO,NOH(5>, II Y, IIX,NSURF, ISOLAT  SU871950 
CaNH0N/C5/FH,FC,NC0LS,N0H IT, ALPHA C10>, I L (50 ) , HCORI 6 ) , ZC0R(6 > , MACH  SUB7196U 


COMttOH/C6/HXCMN( 11 ) , HBCN ( 11 ) , MB IN(ll) , HT{ 90 ) » XCOLL,  YCOLt*  PI »  U 


SUB? 1970 


C0KH0N/C7/C0  (Xfi,  28,2 )  ,NZCO(  18,2 ) , EH»JEX» B2, NHIX,NC1X,H8Q,CB0N, NWCY  SUB71980 


COKhON/C8/IFR^XEC&),YE(3),UX(10),UY(10),HXIHN(ll),El,E2  SUB71990 

C0HH3N/C9/Si X (40, 2 ) , SCXClp, 23 »Y(1T),ETA(11),QWY, QWWX, QHCX, CXHN  SUB7200  0 

R£AD(5,11)(XE( I), 1=1,5)  SU872010 

READC5,11)  (  YE(  F)  ,  1=1, 3  )  ,  SOUND^RHO  SU872020 

RHO  —  RHO/IOOU.O  SUB7203U 

KEAb(5,12)  NHACH, NFKEQ,NHOD£S»  LCOLL >LPRkSH,LPRCO  SUB72Q40 

READ (5, 12)  NUCX, NHPX, NCCXi NCPXiNIONCX  SU8/2050 

RfcAD{5,12)  NIY,NUCY,HPY, INHIS, ISOLAT  SUB72060 

HMTS  =  HHCY*NWCX  +  NLY*NCCX  SUB72070 

DO  40  I=1,NHTS  SUB72080 

40  Ui.(I)  =  1.0  SUB/2090 

n  (INHTS.NE.O)  READ  (5, 11 )  (HT  (  I ),  I=1»N>«TS)  SU8?2100 

HCOLS  =  HPY  *  (NNPX  +  NCPX)  SUB/2110 

HCIX  =  NCCX*N10NCX  SUB7212Q 

M*IX=  NHCX*NIOHCX  SUB/2130 

HUM IT  =  0  SUB/2140 

00  4  1=1,5  SUB/2142 

4  :  Nfljjt  I  )  =  0  SUB/2144 

IF (NUCY.SE^NI Y)  GO  TO  5  SUB7215Q 

MOMIT  =  NIY-HHCY  SU872160 

‘READ (5, 12)  (NOH( I ), ]=1,H0HI1 j  SU872170 

5  READ (5, 11)  (FHACHCI), I=liNMACH)  SUB72183  , 

00  7  1=1, NHACH  SUB/219U 

IMI  HACH(T  ).IF.U.99)  GO  TO  7  SUB7220.0 

HKITE(6,13)  SUB72210 

13  FURMA1C71H  A  MACH  NUMBER  GREATER  THA«>  0.99  HAS  BEEN  READ  IN - SUB72220 

1CASE  TERMINATED)  SU872230 

CALL  EXIT  SUB72240 

7  CONTINUE  SUB/2250 

HEAu(5ill){FRFQU ),  I=1,NFRE0)  SUB72260 

DO  HO  I  =  3',  NHOOES  SUB/2270 

no  20  L.=  1,2  SU972280 

DO  TO  K  =  1,28  SUB72290 

in  Cu( 1 , K, L)  =  0.0  SU872300 

READ ( 5, 12 )  NCO  SUB72310 

N/CU( I , L)  =  NCO  SUR72320 

20  READ ( 5*11 )  ( CO (I , K , l ) , K=1 , NCO >  SUB/2330 

DU  30  1=2,5  SUB72340 

30  XE (  I  )  =  XEU)  -XE(1)  SUB7235U 

YE(2 )  =  YE(?)  -YE(i>  SUB72360 

YE ( 3 )  =  YE(3)  -YE(1).  SUB72370 

11  FORMA  T ( 6E12 .8 )  SUB/2380' 

12  FORMAT(6li2 )  SUB/239U 

IFfNClX.6T.40.DR.NHIX.GT.40)  GO  TO  86  SUB72400 

I F ( NHCX*NHCY+NCCX#N TY . GT  e  9  0)  GO  TO  86  SUB72410 

IK (NPY*(NHPX+NCPX) •GT.50)  GO  TO  86  9c  SUB72420' 


n  a  n  n  n  n.n 


$ 

S 

CGEOM 


Ih(  NWPX.GT  .1  0.  flR.NCPX.GT,10»OR»N'PY.ST*lU  } 
RETURN 

EM  =  FHACH(l) 

CALL  KOUT(i) 

CALL  KOUT !5) 

RETURN 

END  „  « 

FORTRAN  NLStrtUi  DECK 
INCODE  1BHF 
GEOM 

SUBROUTINE  GR»K 
COMPLEX  A,AA,ANM,CZERO 

COMHON/Cl/AtAiDf  AA  ( 50,60  }>AhML 50,10  J^CZERO 


TO  36 


SU972439 

SUB7244B 

SUB/2450 

SUB72460 

SUB7247B 

SU8724S0 

SU872490 

SUB72560 

SUB72510 

SUB72520 

SUB72530 

SU87.2545 

SUB72550 


CUHM0N/C2/CLEN»SN,N1Y,NKCX»NCCX»NWP.'.,  NCPX,  HKERC20  )  ,ZKER  128  }  ,  H6SXRnSUB72560 
CUMM0N/C3/NPY,  SOUND,  NHACHiFHACH  ( 6 ) ,  NFSEO,FREQ{10  )  ,  HAUG,  XIQ.NCX,RHO  SU8/2570 
CUMM0N/C4/NMnuES,LC0LL,LPRHSH,LPRC0,NUH(5),  I!.y,  IIX>NSURF,iSOLAT  SUB72580 
CuHMON/CS/FM, ^C, NCULS, NOH IT, ALPHACIU },IL(5G), HC0RC6 J , ZC0RC6 ),MACH  SUB/2590 


COMMON/ C6/HXCHN! It ) ,HBCN(ll),HBIN(ll),Mr(90}, XCGLL, YCQLL,PI ,U 


SUB7260U 


COMMON/C7/CO(tO,?8,?), NZCOd  0,2), EM,EK, 02,NHIX,NCi X, KBO, CBOX* NHCY  SU872610 
COHMON/C8/IFR,XE!5),YE!3),UX<:iO),U7!10),HXIHN{ll),£l,E2  SUB7262D 
CUHHON/C9/SIX(40,2),SCX(lUi2), Y(il),ETA(il),QHY,OHHX,CHCX,CXMN  SU872630 
WHO  =  MING  ROOT  SEMI-CHORD  SUB?264p 
S  =  SEMI-SPAN  SUB72&5D 
MIC N  =  MING  TIP  CHORD  -  NORMALIZED  ON  MHO  SU37266D 
WHEN  =  MINI)  HP  Lit.  -  NORMALIZED  SUB72670 
Sn  =  SEMI-SPAN  -  NORMALIZED  SUB72680 
can  =  CONTROL  SEN  I -CHORD  SUB72<f90 
Fm  =  2,*NMIX  +  i  SUB72700 
FC  =  2*NC;I  Xtl  SUB72710 
WHO,  =  X  E !  J  )  /  2 . 0  SUB72720 
CLFN  =  Xfc(4  )/HBO  SU972.73U 
S  =  YE  (3)  SIJB72740 
WTCN  =  (XE(3)-XE(2))/MB0  SUB72750 


MUfcN  s  XE(2)/MH0 
SN  =  S/W8U 

CBO  *  !  XE!5)-XE!4><>/2.0 
F1  =  FM 

F2  =  Fl*Pl/2. 

J  =  NMIX 

C  COMPOTE  CHORUM  1  Sfc  INTEGRATION'  AND  :COLLOCaT ION  STATIONS 
C  FIKS1  ON  THE  MING  SURFACE 
DO  5  I  si, NMIX 
F<2  =  F2  -2.*PF 
SIX(J.l)  =  5  I N! F2/FL) 

II  =  FLOA  I  (.!  ) /FLOAT  ( NJ'ONCX)  0,99 
SCX1  11,1 )  =  -S I  X C J,1 ) 

5  J=J-i 
F  L=FC 

F2  =  Fl*Pl/2. 

«T  =  NCIX 

C  THFN:  ON  THE  CONTROL  SURFACE 

00  6  I  si, NCIX 
F2  =  .F2  -2.#P! 

S I  X ( J, 2  )  =  S  I N !  F 2/F 1 ) 

II  =  FLOAT! I ) /FLOAT ! NIONCX )  ♦  0.99 
SCX! 11,2)  =  -SIX! J,2) 

6  .J=J-1 

FI  =  4*Nl  Y 
F2=ti  •  0 

C  COMPOTE  SPANWISF  INTEGRATION  AND  COLLOCA I  ION  STATIONS 


SUB/2760  ! 
SUB72770  j 
SUB72780  | 
SU8/2790  I 
SUB72800 
SUB/281 0 
SUB72820 
SU872830 
SUB72840 
SUB7285Q  \ 
SUB72860  | 
SUB7287|>  l 
SUB72880  | 
SUB7289b  'i 
SUB7290U  | 
SUB72910  | 
,SUB7292b  j 
SUB72930  j 
SUB/294  b  I 
SUB72950  I 
SUB72960  [ 
SUB72970  I' 
SUB/2980  ['■ 
SU87299U  I 
SUB73000>  I 
SUB/3010  I 
SUB73020  L 


26 


c 

c 


12 

13 

14 


16 

66 


T 

S 

cxs 


80  fl  1=1, HIY 
YCH  =  SINCF2/F1) 

F2  =  F2  +PI 
ETACl)  -  SIk(F2/Fl5 
F2  =  F2  *PI 

COMPUTE  WING  SEMI-CHORDS  AND  HID-CHORD  LOCATIONS  AT  THE 
SPANWISE  COLLOCATION  AND  INTEGRAT ipN  STATIONS 
PI8  =  Y£(2)/YE.(3  > 

POB  =1.0-PI8 

CBON  =  (XEC5 )-XE(4 )  )/ (2.0*W80 ) 

CXHPi  =  C80N  <-XFM)/WBQ 
DO  f6  1=1, NIY 

IFCETACD.Lt.BI8j  GO  TO  12 
FI  =  WTLEN©C£IaCI)-PIB)/POB 
If  CYCD.LE.tl8)  GO  TO  13  ' 

F2  =  HTLEN*C  Y  C I J-PIB )/POB 
GO  10  14 
FI  =  0.0 
F2  =  e.o 

WBIkCl)  =  0.5*(2.U-F1) 

NXIHNCD  =  NB INC  D  +F1 


in 

20 

30 

40 

50 

60 


S: 

cbesl 


NLS.TOU,  DECK 
IBHF 


SUB 73030 
SUB73040 
SUB73050 
SUB73060 
SUB73070 
SUB73080 
SUB/3090 
SUB 7 31 00 
SUB7311-0 
SUB/3120 
SUS73130 
SUB73140 
SUB73150 
SUB7316U 
SUB/3170 
SUB73180 
SUB73190 
SUB73200 
SUB73210 
SUB73220 
SUB73230 
SUB 73240 
SUB73250 
SUB73260 
SUB7327J) 
SUB;7328U 
SUB 7 3290 
SUB  7,.33G  U 
SUB73310 
SUB73320 
SU873330 


HHCNCI)  =  0.5»C2.U-f2) 

(CXCHN(.f)  =  HBCN(J)  +  F2 
RETURN 
End 

FORTRAN 
INC OUt 
XS 

FUNCTION  XSCL,NS,13, J3) 

CUMPLEX  A,AA,ANH,CZERO 

CUHHON/Cl/AC6'»},/ACbd,60),ANH(5_0,iO),CZtRO 
C0HH0N/C2/CLEN,SN,N1 Y,NHCX,NCCX,NHPX,NCBX,HKERC20),ZKERC20),NGSKRNSUB73340 
COHKON/C3/NPY,SOUND,NMACH,FHACH{6>,NFREO,FREOC10),HAUG,NiGNCX,RHO  SilB/3350 
CUHHON/C4/NHOi)FS,LCOLL,LPRWSH>LPRCO,NOH(5),  11%  I  IX,NSURf  ,  IS0LAT  SUB73360 
CUHH0N/C5/FW.  FC,  NCOl  S,  NOHi  T,  ALPHA C 1U  )  ,  1 1.  C50  )  ,  HCOR  C6 )  ,  ZCORC 6  )  ,  HACH  SUB/3370 
CUHHON/C6/WXCnNCil),WHCN(ll),WBINCll)>WK90),XCOLL,YCOLL,PI,U  SUB73380 

CUHH0N/C7/CU ( 1 0 ,28, 2  ) , NZCG ( 10 , 2 ) , EM, EK, 02, NH I X>  NC I X, WBO, CBON,NWCY  SUB73390 
COHHON/CB/If  R,XEC5 1, YET3 ) , UXC10 ) , UYC1U ), NX IMNCll ) ,E1, E2  SUB/3400 

C0HH0N/C9/S I  X( 4 0 , 2 ) » SCX ( 10, 2 ) »  Y ( 11 ), ETA ( 11 ), OHY, OWHX ,  QWCX, CXHN  SUB73410 

GO  10  (10 , 40 ) , L  SUB73420 

GO  TO  (20, 30), NS  SUB73430 

XS  =  HXCHNC 03 )  ♦  WBCN  ( J  .3 )  *  SCX(I3,1)  SUB73440 

RETURN  SUB73450 

XS  =  CXHN  +  COON  *  SCX( 13,2 )  SUB73460 

RETURN  SUB73470 

GO  50  (50,60 ),NS  SUB7148U 

XS  =  WXT  M.NC  J3  )  ♦  MBINCJ3)  «  S I X  C  13,1)  SUB73490 

RETURN  SUB73500 

XS  =  CXHN  ♦  CHON  *  SI X( 13, 2 )  SU87351Q 

RETURN  SUB73520 

END  SUB73530 

FORTRAN  NLSTOO, DECK  SU8/3540 

INCODE  IBHF  SU8/3550 

BESL  SUB/3560 

FUNCTION  BESUX)  SUB73570 

If (X.GT.2.0 j  GO  TO  50  SUB73580 

T=X/3.75  SUB73590 

T=T#T  SUB73600 

«S 11=0. 5+ T*<0. 87890594+ r*<0.5149886y+T« (0.15084934+1*10. 02658733+TSUB73610 
1*  CO  . '0  0381532 +>f*  CO.'  0  0  032411  )  ))•)))  SUB73620 

RS11=BSI1*X  ,27  SUB73630 


'4 


Y=X/2..  0 

8SX1=X*AL0G( YXBSU+1.0 

y=y*Y  • 

B$K1  =  BSK1 +  Y*(  II.  154431 44+Y*(-0. 67275579* Y«(-0.18i56397+Y» 
1 ( -0. 01919402+ Y*( -0.00 110404* Y«(-O.Q 0004686 ))))) ) 

RfcSL  =BSK1/X 
60  10  60 
50  Y=2 . 0/X 


SUB73640 
SUB73650 
SUB73660 
SUB73670 
SUB73680 
SUB73690; 
SUB73700 
SUB737-1 0 


BSX1=1. 25331414 *Y*(0. 23498619+ Y«(-0.03655620+Y*(0.01504268+Y*  SUB7372G 

i (-0.0 0780353+ Y*( 6. 00325614+Y* (-0.00 068245) ))) ))  SUB73730 

RfcSL  =BSK1/(SURT(X)*EXP(X))  SU373740 

60  RETURN  SUB/3750 

END  ,SU873760r 

T  FORTRAN  NLSTUU.DECK  SUB73770 

V  INCODE  1BMF  SUB73780 

CCRNL  CRNL  SUB73790 

COMPLEX  FUNCTION  CRNL ( CK , X. Y , CM',  B2 )  SUB7380  0 

C0HH0N/C2/.CLrN,$N,NlY,NKCX,HCCX,NWPX,NCPX,HKER(26>,ZKER(2Q},NGSKRNSUB7381U 


R=ABS(Y) 

R2=K*R 
CKl  =  CK*H 
Gl<  =  0.0 
03  =0.0 
G4  =  0.0 

5  S2  =  X»X  +  U?*R2 
S  =  SORl(S2T 
Ul  =  (CM*S-X)/(B?*R) 

UK  ?  CK1*U1 

HU  20  I  =  i ,  Nl'SKRN 

Ul  =  U1*ZkERU) 

UZ2  =  UZ**2 
0  =  UK*ZKER( !  ) 

F  =  HKFR(l)/  SORT(1.0+UZ2)*UZ*U1 
G3=03*F*CUS(G) 

G4=G4+F*S I N ( G ) 

V  =  1.0  -  ZKFR( I ) **2 

F  =  HKER(  lc)*2 . 0*  V*  fcXP(-CKl*V>/  SORT(t.«*V) 
20  Gi=Gl+F 

G7  =  61  ♦  03 


SUB73820 
SUB73830 
SUB73840 
SUB73850 
SUB7386U 
SIJB73870 
SUB/3880 
SIJ873B9U 
SUB/3900 
SUB/3910 
SUB/3920 
SU8/3930 
SUB73940 
SUB/3950 
SUB  7.3960 
SU873970 
SUB/3980 
SUB7.399Q 
SUB74000 
SU874010, 
SUB74020 


XS  =  X/S 

IF (CK.NE.0.0)  GO  TU  22 
F14  =  1.0 
GU  TO  23 

22  F14  =  CK1*RF SL(CKt) 

23  Gl  =  CKl#G4,-F14~.XS»COS(UKi) 

G2=UK1*G7+XS*SIN(UK) 

XK  =  CK#X 
CO  =  COS  ( XK) 

S!  =  S!N(XK) 

CRNL  =  CMPLX< (CO*G1+S1*G2)/R2»TCO*G2-SI#G1)/R2> 
rejurn 

FND 

f  F  ORTRAN  NLSTtlll*  KfcCK 

5  INCOOF-  I BMF 

CCORO  CURD 

SUBROUTINE  CURD 
CUMPLEX  A,AA, ANM.CZERO 
COMPLEX  AK»  Il2 »  TRM*  01  *  CRNL 

ClJMMON/Cl/A  (  6(1  )  #  A  A  <  50,60  ),ANH<  50,10  ),CZt  RQ 
COMMON /C2/CLFN,  SN . N  l  Y , NW&X » NCCX , NWPX , NCPX , HKER ( 20 
28  CUMM0N/C3/NPY. SOUND, NMACH , FM ACH ( 6 ) * NFREu, FREQ (10) 


SUB/4030 
*  SU874040 

SUB74050 
SUB7406JJ 
SUB74070 
SUB74080 
SUB74090 
SUB74100 
SUB74110 
SUB 7 4 120 
SUB74130 
SUB74140 
SUB74150 
SUB74530 
SUB74540 
SUB74550 
SUB74560 
SUB74570 
SUB74580 
SUB745,90 

)>ZKEPi20),NGSKRNSUB74600 

#  MAUG I ONCX » RHO  SU874610 


CUHMON/C4/HHnDES,LCOLLiLPRHSH,LPRCO,NOH(5),  JH.  I!X,NSURF,  ISUlAT 
C0MH0H/C5/F-Jf»  FC , NCOLS*  NON IT,  ALPHA C 10 ) ,  I L (58  >, HCOR ( 6 ) ,  ZCOR  C  6  )  ,  HACK 
CUHH0N/C6/HXCHNC 1 1 )  ,  WHCN  C  IT)  ,  K8  IN(li),  KT{90  )*  XCOLL,  YCOL.L,  P I ,  U 
CUHHbK/C7/eO{10;28,2),NZCO{18,2),EM,EK,H2,NWlX,NCIX,HBO,C8ON,NllCY 
CG8H0N/C8/IFR,XE(5l),  YEC3),UX(lfl),UY(i8)>HXIHN{il),El»E2 
CUHH0H/C9/S 1 X ( 40 , 2 ) , SCX { 10 , 2  }  ,  Y  C 11 ) , ETA  C 11 ) ,  GWY,  Qk'ifX  ,  QHCX,CXMN 
THIS  SUBROUTINE  CONSTRUCTS  A  HOW  qF  THE  DQWNRaSH  MATRIX 


THE  PRESSURE  SERIES 


PRODUCT  OF  CHEBYSHEV  POLYNOMIALS  IN  THE 


13 


NEGATIVE  OF  PERCENT  SEMI-CHORD  FROM  THE  MID-CHORD  AND  PERCENT 
SEMI-SPAN  FROM  THE  ROOT. 

DO  6  JC=1,NCQLS 
A (JC)  =  CZEKO 
ict  =  o 

NIX  =  NWI X 

QKX  =  -0WWX»SH*«2/(8.fl*PI7  * 

NPX  =  NWPX 

THE  DO  14  LOOP  COHPOTES  THE  NON-SINGULAR  PORTION  OF  D(N,H) 

OUF  TO  BOTH  SURFACES 
DO  14  HSUKF=1„2 

IKhSURF.NE.NSURF.ANO.ISOLAT.NE.fl  i  GO  13 
OU  12  IY=1,NIY 
EfAT  =  SN*'ET  A  ( I.YT 
E 12  =  ETAJIY)**2 

IF ( NPY.GT .1 )  CALL  CHEH ( NPY-1,ETA { I Y > , UV(2 ) ) 

ilYU)  =  1.0  -FT2 

DO  3  K=2,NPY 

UY(A)  =  E  T2-*ilY  C 1 )  *UY  { K  ) 

DO  10  I  X=1 ,  N.I  X 

XI  =  XSC2, HSURF  <1 X, TY ) 

XID  =  XCOl.'L  -XI 

AK  =  CRNL(EK,XID»YC0LL-ETAI,EN,B2>  ♦  CRnL C EK, XI D, YCOLL+ETA I , EH, B2 
IC  =  1C1  *1 
HZ  =  AK*OHX*(JWY 

IF(NPX.GTil)  CALL  CHF.H(NPX-1,-SIXC  IX, HSURF  )  ,UXC2)  > 

UX(  1 )  =  1.0  -SIXCIX, HSURF) 
no  4  K=2»NPX 

J)X(K)  =  Cl. 0  *SIX(  IX, HSURF)  )*UX< 1 >*UX C K ) 

«*  ADD  AN  ihCRFMENT  TO  EACH  ELFHENT  OF  I  HE  ROW  FOR  CXI, ETA!)  ** 

DO  10  NY=1 ,  NPY 
TRM  =  H2  »  UY(NY) 

00  30  NX=l,hPX 

AUC)  =  A  (TO  +  TRM*UX (  NX  ) 

IC  =  ic+i 

**  IC  EQUALS  NPY*NWPX+1  AT  1  HE  END  OF  TmE  FIRST  PASS  #* 

CONTINUE 
NIX  =  NCTX 

OWX  =  -QWCX*SN**2/<8.0*PI ) 

NPX  =  NCPX 
IC1  =  NP Y*NWPX 
I  Cl  -  0 
NPX  =  NWPX 

XCOIS  =  <XS(  1 }  HSURF,  T  I>X,  MY) 

Y2  ?  Y(IiY)**2 

CALL  CHEBC  NPY-1 >  Y ( I  I Y ) , UY ( 2 ) ) 

II Y Cl )  =  -2.0 
l)U  15  K=2,NPY 
IJY(K)  =  -2.0«Y2*IIYOO 
|i0  40  MSURF  =  1,NSURF 

»*  THIS  LOOP  ADDS  THE  CONTRIBUTION  OF  THE  SINGULAR  INTEGRAL 
ALONG  THE  LINE  FROM  THE  WING  L.E.  TO  TME  COLLOCATION  POINT 
IF ( HSURF • NE. NSURF. AND. I SOLAT .NE. 0 )  GO  Tu  23  29 


SUB74628 
SU874630 
SUB74640 
SUB74650 
SUB7466B 
SUB74670 
SUB74680 
SU874690 
SUB74700 
SUB74710 
SUB74720 
SU874730 
SUB74740 
SUB74750 
SUB74760 
SUB74770 
SUB74780 
SUB/4790 
SUB7480U 
SUB74810 
SUB7482U 
SUB7483U 
SUB74840 
SUB74850 
SUB74860 
SUB74870 
SUB7488Q 
SUB74890 
SUB  7  4.9  00 
SUB74910 
JSUB74920 
SUB74930 
SU874940 
SUB74950 
SUB74960 
SUB7497U 
SUB74980 
SUB7499U 
SUBZ5000 
SUB 750 10 
SUB75020 
SUB75030 
SUB75U4U 
SUB75050 
SUB7506U 
SU875070 
SUB/5080 
SlJB/5090 
SUB/5100 
SilB75110f 
SUB 7512 6 
SUB75130 
SUB/5140 
SUB75150 
SUB75160 
SUB75170 
SUB75180 
SUB75190 
SUB/5200 
SUB/5210 
SUB/5920 


IFINSURE.LE.HSUHT )  GO  T 0  16 
UPLIH  =  PI 
QU  10  IB 

16  XI  =  SCXU  IXiNSURF), 

IIPLIH  =  -ATARI  SORT  ( 1  .'0-XT**2  )/XT  ) 

I F ( UPL I M. IT . 0*0  >  UPL l K=UPL1 K*P I 
18  QwSNG  =  FLOA T( 2*N I Y ) *UPL I  M/8 » 0  , 

liu  22  N=l»  6 
IC  =  ICl+1 

C  **  THIS  LOOP  CONSTRUCTS  0(0*H)  . H=D  ,4  , . . . .NPX-1 

V1N1  =  UPL1H#ZC0R(N> 

C  =  COS(V.lNf) 

CALL  CHE8 { NPX- 1 »  C>OX<?)~) 

U X ( 1 )  =  1.0  *c 


SUB75230 

SUR75240 

SUB75250, 

SU8752vO 

SUB7527U 

SU875280 

SUB75290 

SUB7530U 

SU8/531Q 

SUB/5320 

SUB^5330 

SU875340 

SUB75350 

SUB/536U 


no  19  K=2VNPX 

19  UX(K)  =  <1.0  -C)*UK<i  )>UX<K) 

AMR  =  EK*  (  XCOl.S  -HXCHNUIY)  ♦C*WBCN  (  1 1 Y  )  ) 
IMHSURF.L0.2)  AKG=EK*(XCOLS»C)fMN  +  C*CB0;O 
Cl  =  COS(ARU) 

SI  =  SORT  Cl  .fi  -Cl* *2) 

1)1  »  CMPlxtCl  ,-Sl)«HCOR(N) 

no  22  NY=i,NPY 

TKM  =■  OHSNG*UY(NY)«l!l 


SUB75370 
SU87538U  ' 
SUB7539C 
SUB  7539,5 
SUB/54D0 
SUB 75410 
SUB7542Q 
-SUB/5430 
SUB7544U  , 


2? 

23 

4n 


% 

* 

cznzx 


t)U  22  NX=1  j  HPX- 
AUC)  =  A  (,;TC )  ♦IRH*UX(NX) 

rc  =  ic*i 

l  Cl,  =  NPY*NHPX 
NPX  =  NCPX 
RETURN 
END 

FORTRAN  NLSTtUl.DbCK 
1 NCOOE  IBHF 
ZDZ* 

SUBROUTINE  ZnZXtMiiOfc,  OZiZ) 

COMPLEX  A » AA >ANM> CZERO 
OUMMOn/CI/ A  (  6ft  ) »  A  A  <  bO »  6(1 ) »  ANH  ( 50  »  10  )  »  CZERU 


SUB/5450 
SU675460 
SUB 75470 
SUB75480 
SUB/5490 
SUB 75500 
SUB7551 0 
SUB74160 
SUB?4170 
SUB74180 
SiJB74l9u 
SUB/4200 
SUB74210 


CUMMUN/C2/CLFU.SN.NI Y , NWCX, NCCX/NWPX, NCPX, HKFP ( 20 ) , ZKER < 20 ) , NGSKRNSUB/4220 
CUMM0N/C3/ NP Y, SOUND, NMACH, FM ACH < 6 ).NFREU,FREO (id), »MAUG.NIONCX.RhO  SUB74230 

C0MM0N/C4/NM0HES,  LCULL»  L PRMSH#  LPR.CO#  N.OM < 5 ) » 1 1  Y »  1  lfX»  NSURI  ,  1SOLAT  SUB7424Q 
COMMON/C5/FW,  FC,  NUO|.S»  NOMIT#  ALPHAdO  )  #  II  < 50  ) ,  HCOB< 6  ) ,  ZCOR < 6  ) ,  MACH  SUB74250 
C0MM0N/C6/HXCMN  <  l‘i  )  #  WBCN(  11 )  #  HBJ  N  ( 11 )  I  ( 90.)  »XCULL»  YCOLL,PI»U  SUB  74260 

CUMM0N/C7/CUt ■10»28.2),NZCO(10#2).EKJEK»m2#NMIX'#'MPIX#HBO*CBOH.NHCY  SUB74270 
CUMM0N/C8  /  1 1  R  ,  XE  (b  ).,  Yf: ( 3  ) »  UX ( 10  ) ,  U Y <  10  ) ,  HX’I MN ( 11 )  »  El » E2  SUB/4280 

COMMON/ C9/S  1  X  ( 40 »  2  ) »  SCX,<  1 0># 2  ) »  Y (1 1 )  » ET A <  11 )  #  OH V,  QWWX,OWCX,CXMN  SUB/4290 

PX  =  XCOLL  *  WBO  SU8743C0 

IFtNSURF  .GT.l)  P Xr P X - C L E N * Vt B 0  SUB/431;', 

PY  =  YCOLL  *  (BO  SUB/ 4320 

Z=  COtHODE.l ,nSURF)  SUB7433Q 

nz  =  6.U  SUB74340 

K  =  2  SUB74350t 

YX  =  PY/PX  SUB74360 

PPX  =  PX  SllB/4370 

nu  40  N  =  2,7  SUB/4380 

PX  Y  =  PPX  SUB/4.390 

F  =  N-l  SUB74400 

no  30  M  -  1,N  SUB/4410 

1 F ( K » G T • NZCU ( MODE » NSUHF ) F  G0‘  TU  50  SUB/4420 

,ZP  =  PXY  *  CO(MOnE»K»NSURF)  SUB/4430 

Z  =  Z  ♦  ZP  SUB/4440 

nz  =  DZ  ♦  ZP  *  F  SUB/ 4450 


30 


PXY  =  PXY*YX 
F  =  F  >r  1.0 
30  K  =  *  * 1 
40  PPX  =  PPX  <*  PX 
50  1)2  =  DZ/PX 
RETURN 

END  , 

%  FORTRAN  NLSTiUl/DHCK 

$  INCODE  i BMF 

CCGRFi)  CGREDr 

SUBROUTINE  CORFU (V,  IRf  JC) 

COMPLEX  ANH, CZFRO 
DIMENSION  V(P.l) 

COMMON/ Cl~/ A  {  2 »  601 » AA ( 2j#  50  >  60  ) , ANH(5U, 10 ) , CZERO 


SUB74460 

SUB/4470 

SUB/4480 

SUB74490 

SUB74500; 

SUB/4510 

SUB/4520 

SUB 75520 

SUB75530 

SU875540 

SUB75550 

SUB75560 

SUB7557U 

SUB/5580 


CUMMUh/C2/CLEN,SN,NJ Y;NWCXiNCCX,NWPX,NCPX,HKER(20)*ZKER(20)>NGSKRNSU875590 
CUMMON/C3/HPY»SUUND»NRACH#FMACH<6)*NFREu,FREQ(io)»  HAUG*  N I ONCX » RHO  SUB/5600 
COMMON/ C4/NM0 OF S.«  LCOLL, LPRRSH, LPRCO, NOMI 5  ) » IT  Y ,  I  IX,NSURF  ,  I  SOLA!  SUB/5610 
CUMHON/C5/FR, F  C»  NCGLS*  NUM  LT#  ALPHA  (10) »  IT  (50  )  ,  HCUR  ( 6  > >  ZCUR(  6  >  ,  HiCFI  /SUB7562U 


CUMMGN/C6/WXCHN(ll),HUCN(li),  MBIN(  11  )  »  N  T (  90  )  ,  XCOLL,  YCOLL*  P  l»  U 
COMMON/C7/COMO  #  28, 2  ) » NZCOdfl »2  ) » EH» EK» l*2» NHI  X,  NC I  X»  WBO»CBON,NWCY 
C0HM0N/C8/  IF  R ,  XE ( 5  ) »  YE('3  V<»  UX  ( 10) ,  UYdOJ,  MXIHNdl  hr  El.,  E2 
C0MM0N/C9/S 1 X <  40 , 2 > , SCX .10 , 2  > ,"(11) » E  TA (111 , QWY, URWX ,QWCX, CXHN 
RrtN  =  SUK1TA, A <  1 ,  IR,  10**2  ♦  AAT2,  IR*  ICJ**2) 

IF  ( A  A  ( 2 ,  ;1R  >  1 C )  .1 E  •  E? )  GO  TO  30 
Ck  =  A  A  ( 1 ,  IR,  IO/KMN. 

Cl  =  AA(2, Ik, Icj/RMN 
1)0  20  N=  I C ,  M  A  lift 

Tl  =  CH*AA(1. IR/N)  +  C I *AA (2  * IR*  NT 
AA(?,JR,N)  =  CR«'AA(?,1R,N)  -  Cl*AA(l,  IR,N) 

AAd.lR.N)  =  II 

RAN  ~  SORT ( V( 1 , 1C )**?  *  V(2, IC>**2  ) 

Lf4rAN.LE.E2,)  .  GO  TO  60 
RAN  =  SORT (RAN**?  ♦  RMN**2) 

CR  =  y ( 1 # 1C j/RAW 
Cl  =  V  ( 2 , I C ) /RAN 
RMN,.  =  RMN/RAN 
110  50  N=  I C »  MA  U(t! 

AIR  =  RMN*AA  (1  »>I  R,>N )  ♦  CR*Vd,N>  *  CI*V(2,N) 

All  =  RMN* AA  f  2 , 1 R»  N )  ♦  CR*VE2,N)  -  Cl*V(i,N) 

VH  =  RMN* V ( 1 »  N )  -  CR*AA (1* IR*  N )  +  C 1  * AA 1 2 , 1 R» N > 

VI  =  RMN*V(2,N)  -  CR*4A  ( 2  >  I  R  »  N  )  -  C 1*  A  A  ( 1,.I  R,  N,) 


VI  =  RMN* V( 2 ,  Ni)  -  CR*4A  ( 2 »  I R ,  N  )  -  C 1  *  A  A  d ,1 R , NT  SUB/6850 

A  A  d  ,  I  R ,  N  )  =  AIR  SUB75860 

A  A ( 2  *  I R  *  N  j  =  All  SUB75870 

V(1,N)  =  VR  SUB75880 

50  V!  2  »‘N  )  =  VI  SUB/5890 

6 d  RETURN  .  SUB/5900 

END  .  SUB  759,10 

$  F  ORTRAN  NLSTOU,  DECK  SUB7592U 

!S  INCOUE  I  BMF  SUB75'930 

CXESO  XL SO  SUB/5940 

SuBkOlJT  I  Nt  XL  SO  SuB/5950 

COMPLEX  AVAA, ANM.CZLRO  SU875960 

COMMON/CIV  A  ( 6il ) »  A  A  ( 50  »  60  ) ,  ANM(50 , 10  > ,  CZbRU  SUB  75970 

COMMON/C'2/CI  EN,SN,NIY.NWCX,NCCX,NWPX,NCFcX,HKER<2U>,iKfcR(2Ii»»NGSKRNSUB/598U 
CUMH0N/C3/NPY, SOUND, NMACH,FMACH<6>,  NFRE*,,  FRfcO (1 U  >  ,  MAUG.  N !  ONCX ,  RhO  SIIB75990 
COMMON/C4/NMOMFS,  LCOI.L ,  LPRMSN*  LPRCO, NOM ( 5  ).,  1  I  Y,  I  I  X,  NSURF  ,  ISOLAT  SUB/6000 
C0MM0N/C5/FW,  EC,  NCOLS,  NOM I  T»  ALPTIA.CI  0  >  ,  Ii.  (bO  1,  HCOR ( 6  ) ,  ZCUR  (  6  ) »  MACH  SUB/6  01U 
C0MM0N/C6/WXC1N(U),WHCN(11),WBIN(U)/WI  ( 90  ) ,  XCOLL,  YCOLL ,  P I »  U  SUB/6020 

COMMON/ C7/CU  (10, 2d, ?),  NZCOdO,  2  ) ,  E  M,  EK ,  *2i  Nwl  X,  NC I  X,  WHO,  CBON,  NwCY  SuB/6030 
C0MM0N/C8/1 F  R,  XE  (  5  )  #.'Yb  ( 3  )  ,  UX  ( 10  ) ,  UY  ( 1  U  ) ,  WX  I MN  ( 11 )  ,  El ,  E2  SUB76040 

C.UMMON/C9/S1X<40,?),SCX(10,2),  Y(ll  ),EIAH1  ) ,  QW  Y ,  OWWX ,  OWCX ,  CXHN  sub/6  050 


iSUB/5630 

'SUB/5640 

iSUB/5650 

SUB/5660: 

SUB/5670 

SUB75680 

SUB75690 

SUB/570U 

SUBZ5710 

SUB/5720 

SUB75730 

SUB75740 

SUB75750 

SUB75760 

SUB 75770 

SUB/5780 

SUB/5790 

SUR/5800 

SUB75810 

SUB / 5 82 U 

SUB/5830 

SUB/5640 

SUB/58f>0 

SUB75860 

SUB75870 

SUB75880 

SUB/ 5890 

SUB/5900 

SUB 75910 

SUB/592.U 

SUB 7 5930 

SUB/5940 

SUB/5950 

SU875960 

SUB  75970 


11  =  1 

SUB?6Q6U | 

nu  136  l=leNCOl s 

SUB76070 | 

Ri I  =  CAfiSlAAUl.I ) ) 

SUB/6880 I 

fF'Ckl'I.LE.E?)  BO  10  135 

SUB76090 | 

I L ( 1 )  =  11 

SU876180 

11  =  11  ♦  1 

SUB 7611 0 

GO  10  136 

SUB 76120 

135 

I L ( 1  )  =  -1 

SUB76130 

1 12  =  NCOLS  -  1  -  (1-11) 

SUB76140 

DO  1135  11=11 .1*2 

SUB 7615 G 

1135 

CALL  CGKfclJlAAl ll  +  l.i).ll.l*l) 

.SUB/ 6160 

00  2135  L=l,  NHOiJFS 

SUB 76170 

ML  =  NCOLS  ♦  L 

sub 7618 o 

2135 

ALPHA ( L  )  =  SORT ( ALPHA ( L ) »»2  ♦  CABS ( AA  ( ! J 2+1 , HL )) **2 ) 

SUB/6190 

136 

CONTINUE 

SUB76200 

C 

SOI  VE  FOR  THE  CUEFF  ICIFNTS  BY  BACK  SUBSTITUTION 

SUB76210 

14  0 

11  =  NCOLS 

SUB 7 6220 

nu  150  1  =  1, NCOLS 

SUB 7623 G 

no  350  l=i>nhhdes 

SUB76240- 

150 

AnH(1,L)  =  CZERU 

SllB/6250 

no  210  J=l. NCOLS 

SUB/6260 

IF  (  1L(  I  D.LE.M)  BO  TO  210 

SUB7627U 

J1  =  IL(M) 

SUR/-6280,' 

110  20 U  L=l.NHODES 

SUB76290 

HL  =  NCOLS  *  1. 

SUB/630  0  c 

imi-NCOLS)  170,1911,220 

SUB76310 

176 

IK  =  II  ♦'  1 

SUB76320 

no  180  K=IK, NCOLS 

SOB/633U 

lan 

SNH(II,L)  =  ANHUI, L)  -  AAUI,K)*ANH(K,L) 

SUB/6340 

190 

ANHIIL.L)  =  UNM(IIA)  +  AA.(J1,HL))/AA(JI,II) 

SUR76350 

2«o 

CONI INOE 

S0B76363 

210 

11=11-1 

SUB76370 

22  0 

RtTORN 

SUB763C9 

END 

SUB76390 

* 

FOkIRAN  NI.S.T  Oil,  DECK 

SUB 7 64 00 

s 

1NCOUE  1 BHF 

SUB76410 

CCHER 

GHE8 

SUB 76420 

SUBROUTINE  CHER <N1 » X , UX )  ; 

SUB76430 

III HENST  ON  UX(1) 

SUB/6440 

no  io  i =r, ni 

SUB  7  645,0 

in 

IJX(  1  )  =  0.0 

SOB7646.0 

ux(i)  =  l.o 

SUB764710 

Ox  ( 2 )  =  2 . n  *  x 

SU8/6480 

11  (Nl.LT* 3)  RETURN 

SU876490 

no  20  1=3 , Nl 

SUB76500 

20 

IJX(I)  =  2«0*X*UX(  1-1  )  -  U  X  (  I  -  2  ) 

SUB7651U 

RETURN 

SUB 76520 

END 

SUB 76530 

s 

FORTRAN  NLSTUU, DECK 

SUB  7  654  0. 

% 

INCODE  1 BMF 

SUB76550 

OP  ORO 

FORC 

SUB76560 

SUBROUTINE  F  ORC (  NW  I CX,  NC I CX,  N 1 CY ,  YE  1 A ,S  1 CX ,  WB 1 CN,  LPR  ) 

SUB 76570 

CUHPLEX  A,AA,ANH, C ZERO, GF0RC,9ELP, HASH, APR , UWASH,  CRASH,  PR 

SUB76580 

DIMENSION  Yl  TA  ( 1 ) ,  S 1  CX  ( 11 ,  WB'l  CN  ( 1 ) 

SUB  76590 

DI  HENS  ION  DWASH!9U>  1  O'),  PR(  90,10  ),  CRASH!  *0,10) 

SUB7660  •» 

D 1  HENST, ON  Gf  ORC  ( l  U  ,  1  0 , 3  ) ,  HELP  ( 1  0  ) ,  HASH  (  Ifl  ) 

SUB70370 

COMMON/Ci/A ( 60 ),AAf 50,60 ),ANM( 50,10 F/CZERO 

SUB76620 

CUHM0N/C2/CLFN.SN.NI  Y, NWCX, NCCX , NHPX, NCRX, HKER < 20 ) , ZKER < 20 ) 

, NGSKRNSUB/6630 

CUHM0N/C3/NF  Y,  SOUND,  NM  ACH ,  FM  ACH  (  6  ) »  NFREu  ,  F  REQ  <  11) ) .,  MAUG,  N I ONCX ,  RHQ  Si B76640 
CUHM0N/C4/Nh0UtS.l.CuLL,LPHWSH»LPRC0*N0M(5) ,  1.1  Y»  I  IX,NSURF  ,  I  SOLA!  SOB  76650 

32- 


CUMH0N/C5/FW,FC#  NCOLS# NON! T» ALPHA (10 ) #  II.  (58  ) #  HC0R(6 ) # ZC0f*(6)  »HACH  SUB/666.0 
:CUHHbN/C6/WXCHN(ll  ).#  WBCN(ll)#WBIN(ll)#WT  (90  )# XSQLL#  VCOLL#  PI #  U  SUB76670 
CUHH0N/C7/C0(  1  0»28»2>#  NZCOdO#  2}»EH#EX#  H2l#  NWI  X#NCIX#  WBO#CBON»NWCY  SUB76681J 
CUHH0N/C8/  If  R,  XEL5) *  YE(3)  #  UX{10  > #  UY( 10 ) #  MX  I HN ill  ).*E1»  E2  SUB76690 

C0Mh0M/C9/S  1 X( 40#2 ) »  SCX ( 10  # 2 )#.Y ( 11  )  #  ETA( 11 }  »QWY# QWWX, QWCX#  CXMN  SUB76700 

COHHON/.CPR/APH (90,600 #  I  HOD #1 ROW  SUB 7.671.0 

EUU1  VALENCE  (tiFORC# AA ),  ( A, WASH ),  (.WASH, Df  LF>  SUB7672U 

EQUIVALENCE  I HWASH(3.»1 )» APR(l#5l)  )# (PR(  1 » 1 )  #  APRM#1T) )  SUB76730 

EQUIVALENCE  CCWASH(l»l)#ApR(l#4lI>  SUB76740 

Q«F  =  tt.54RHn*Mj*HBU*SN)**2  *QWY  SUB76750 

TF(LPR.NE.OT  NO  TO  2  SUB76760 

DU;  i  l=i#NHODFS  SUB76770 

ftli  1  J=l# RHODES  SUB76780 

fiU  1  K=l-#2  •  SllB/67^0 

1  GFORC( T# J#K)  =  CZERU  '  SUB7680U 

?  ICJ  =  0  SUB76810 

NiCX  =  NWICX  ‘  SUB/6320 

NPX  =  HWBX  SUB76830 

QWFURC  =  UWF  *UWWX  SUB/6840, 

F  HOW  4  i  SUB76850 

LflNli  =  40  ;S*JB7&860 

If (LPR.NE.B)  I  MNT1  =  Lfi  SUB/687U 

no.  1000  NS-1 .2  SUB 7 6880 

NSUHF  =  NS  SUB76890 

Pf  AC  =  SN/CUON  SUB 7690 0, 

nil  900  1 1 Y  =  1 #  NICY  SUB76910 

YCOLL  -  YE-TA(.HY)  StJB/6920 

»*  YCOLL#  AND  LATER  XCOLL#  ARE  USED  HERE  TO  DENOTE  INTEGRATION  SU876930 
STATIONS  BECAUSE  WE  USE  SUBROUTINE  Zi>ZX  TO  COMPUTE  THE  SUB7694U. 

DISPLACEMENT  THROUUH.  WHICH  THE  PKESSuRE  ACTS  TO  DO  WORK  SUB/6950 

Y2  =  YC0LD»*2  SUB  7.6960 

If (NPY.GT.l)  CALL  CHER<NPY-1#YC0LL»UY(2>)  SUB7697IJ 

(I Y ( 1 )  =  1.0  -,Y2  SU876980 

IF  ( 1-PKiNE.  O  )  UYIl  j#S.ORT(i.O-Y2  )  SUB/69?U  ; 

00  3  K=2#NPY  50877000 

3  II Y  C  M  =  Y2»UY(1)*I»Y;(K)  SlJB77*JlU 

IFTNSURf  ,E0.1  )  PF:AG=SN/W8  I CN(  1 1  Y  )  SUB770^0 

YCOLL  =  SN# YCOLL  SU877025 

DO  80  0  .IT  X=l#  N I CX  SU8Z7030. 

LMN !  =  MX  tLMNii«  (?NSURf-l  )  SUB/7040 

XCOLL  =  SICX(LKNT)  SUB/7050 

I  F  (NPX.GT.l )  CALL  CHEhT NPX-l , tXCOLL# UX(? ) >  SUB77060 

U X ( 3  )  =  liO  -XCOLL  S.U877070 

IF  ( LPR.NE. 0,>  UX(1  )=SQR:T(UX(1  ),/(!. O  +  XCOLI  )  )  SUB77U8U 

Du  4  K=2-»HPX  SUB77090 

4  IJX(h)  =  (1.0  ♦XCOLL)»UX(l)«OX(K)  SUB/7100 

IQ  =  1  Cl  +1  ,  SilB77110: 

DO  10  .1=1#  NHODES  SUB77120 

10  nEI.PTj)  =  CZEHO  '"SUB7713U 

nu  20 U  NY  =  1#NI>Y  SUB77140 

1)0  200  NX=i#NPX-  SUB7.715.0 

DO  20  J=1 #  NMODES  SUB77160 

nEl.P(J)  =  DELP(J)  ♦OX(NX)»UY(rsY)*ANH(  IC,  J)  SUB/ 71 70; 

If  (LPR.NE. Ob)  PRUrOw.J)  =  DE!.P(;j)*PFAC  SUB/ 7 180 

20  CON  i  I  NUE  SUB7719B 

200  IC  =  1C+1  SUB/7200 

**  IC  =  N P Y * N W P X  +T  AT  THE  END  OF  THE  NS=1  PASS  SUQ'721.0 

AND  DP  CON  I  A  INS  DELTA  R/CI  At  ( XCOLL#  V  CULL >  SU07722U 

IRON  =  FROM  *  V  SUB77230 

If  (LPR.NE.O  )  GO  TO  800;  SUB77.240 

XCOLL  -  XS(2,NS,  I  TX>:i  IY.)  33  SUB/7245 


I, 


<5 


DO  40U  { =1, NHGDES 

CALL  ZDZX( 1  *  SLOPE* D ISP) 

DO  4011  J=l, RHODES 

40fl  GFOKCCi.J.NS)  =  GFOKCt I , J,NS )  *aUFOKC*D ISP*BELP(  J)*2.8 
800  CONTINUE 

900  CONTINUE 

N1CX  =  NCICX  « 

NPX  =  NCPX 
I  Cl  =  IC->1 
OWFORC  =  UWF »UVCX 
1000  CONTINUE 


s 

s 

CKOUT 


Return 

eng 

FORTRAN  NLSTnll,  DECK 
f NCOUE  sbhf 
KOUI 

SUBROUTINE  KOiiTt  1N0) 

COMPLEX  A,AA,>4NH,CZER0,GFURC,DELP,HASH,APR,DWaSH,CWASH,PR 
01  MENS  I  ON  CAR»S425,50  ) 


SUB77250 

SUB77260 

SUB/7270 

SUB7728U 

SliB/7293 

.SUB773Q0 

SUB77310  I 

SUB7732U  | 

SU877330 

SUB7734U 

SUB/7350 

SUB/7360 

S0B7737U 

SUB/7380  1 

SU87739U  § 

SUB/740U  | 

S0B7741U  1 

SUB7742U  I 

********  i 


DIMENSION  GT  n»C  (-10,10,3)  #  DELP ( 1 0  A  *  WASHflG  >  SUB70370 

D  IMENSION  DWASH( 911  >  TO  ) ,  PR( 90 , 10  )  >  CWASHfyQ,  TO  )  SUB/7440 

DIMENSION  SURF  f;2 » 3 ) »  XPR ( 50 )  SUB774$0 

COMMON/Cl/ A  ( 60  ) ,  AA ( 50 , 60  ) ,  ANH.(50r.>  1 0  )  > CZERU  SUB77460 

COMMON/ C2/CI.FN,  SN>  N IY, NHCX, NCCX,  NWPX, NCPX.HKERf 20  >  ,ZKER ( 20  )  »  NGSKRNSUB/7470 
COMHON/C3/MPY:,  SOUND*  NHACH, FMACH( 6  ) ,  NFRE'j,  T  REQ  { 10  )  »  HAUG,  N I GNCX, RHQ  SUB7748Q 
C0MH0N/C4/NHniiFS»  LCULL, LPKWSH, LPRCO » NUH  ( 57  » II Y#  1 1  X>.NSURF  ,  1SOLAT  SIJB77490 
COHMON/C5/FR,FC,NCOLS,HOH IT, ALPHA £1U),  It  (50 > , HCOR( 6 ) , ZCOR( 6 )  ,  MACH  SUB77500 
CUMM0N/C6/KXnHN(Ll),H0CN(ll),NBlN(ll),Hl ( 90  ) , XCOLL, YCOLL, PI » U  SU87751U 

COMMON/C7/CO  (10»28*2),NZCGilO,2:),EH,FK,h2»NWlX»NCIX»  WHO, CBON,  NHCY  SUB/7520 
CUMhON/C8/IFR,XE<5),YE<3),UX<10),UY(lU),UXIHMai>,El,£2  SU877>J0‘! 

CUHM0N/C9/.S  I  X  C  40*  2 )  ,.SCX (,1U,2  )  »  Y  ( 11)  , ETA { 11), QHY, QWWX,  QWCX, CXHN  SU877540 
CUHhON/CPR/APP(90,60),TM6D, IRON  SUB77550 

EOtll  VALENCE  r.»;FORG,  aA  ) ,  ( A,  WASH  )  >  ( WASH*  DELP  )  SUB/7560 

EUUI  VALENCF  1  »iWASH<  1  ,  1 ),  APR (1,51 ),),( PR(  1 )»  APR(1.  1W1  SUB77570 

FUU1  VALENCE  (CHASHd  »1)>  APR (1,41))  SUB77580 

EUIII VAIENCE  f  XPR,  It.)  '  SUB77590 

DATA  (SURF<1, I), I=l,3)/bHWING  ,  8HTA I L  *  1 1HHTNG  +  I  AIL  /  SU877600 

GO  10  (10, 20,30 >40,50, 60, 70,80, 90)>lND  SUB77610 

.*/*** *i  *#***.«***«:*#**  ******  ##«*##«*#******#** **«**#»*•***« #*«*#*«#**SUB 77620 
in  X  V  =  XF  (5-i  -r  XF(44  SUB7  763U 

*X  =  XE  ( 3 )  -  XF(2)  SUB7764U 

Ah  =  -2  •  0*XE  ( 3 )  *  YF  ( 3)  -  XE (2 >* ( Yt < 3 >-YE< 2 > )  SUB/7650 

AI  =  2«  0*XV* YE  M )  SUB/766U 

WK 14  E(  6,  IT) L  M>  SOUND,  RHfl»>XE  (1 ) »  XE  ( 4  )  >  XE (  3  ) ,  X  V,  YE ( 2  )  >  YET 3  )  >  YE ( 3 )  >  SUB7767U 

lYh(J), XX, XV, AW. AT,NWCY,NI Y , NWCX , NCCX, NW I  X, NC I  X , NPY , NPY „ NWPX, NCPX  SUB77680 
1*  FORHAf(lHl/////  3iX, 41HHAC/NAA  MISSILE  SUBSONIC  AIRLOADS  PROGRAM  SUB7769U 
1  ///37X»30HKI  fGHT  CUNUlIIONS  AND  GEUHETRY/1H0//15X,  13HMACH  NUMbERSUB/7700 
T  s„F 8.5,4X, 1 6USPEED  OF  SOUND  =F1Q.3,4H  l / T , 4X , 4HRH0= , E14 . 8//1H0/  SUB7771U 
X54X,4HWING,18X,  SUB/7720 

3-  4HlAIL///22X,16HL.fc.  STATION  ( L  ) , 2F 22.3//22X , 16HR00T  CHORD  ( L  ), SUB 7 7 730 

2F22.3//  22X.16HL.E.  SPAN  ( L ) , 2F22.3//22X i 16HT.E .  SPAN  (L), SUB/7740 

2F22.3//  22 X, 4  6HT I P  CHORD  ( L  )  , 2F22.0//22X .16HTOT AL  AREA  4L*L)  SUB7775U 

2F22.3//22X,  16HSPAN  COLL,  STA . , I IV,  122 , //22X, 16HCH0RD  CULL.  ST  A.SUB7776U 
119,  1 22//22X , 1 6HCHURD  I N 1 6 .  STA . ,4 19 , I 22//22X , 16HSPAN  PRES  MODESSUB/7770 


4 

5 
b 

7 

8  1 


MODES, 119,122) 
TO  15 


119,  122/ /?2X , 16HCH0R0  PRES 
I F ( T  MACH( MACH  V. LF • 0 .95 )  GO 

WKI  IE«(6,T4.) 

14  FURMA  T  ( 9.2H  A  MACH  NUMBER  GREATER  -THAN  0.95  HAS  BEEN  USED- 

1IJSE  CAUTION  IN  APPLYING  CASE  RESULTS) 

15  IKNOMI  I  .EO.O  )  RETURN 


SUB  77780 
SUB77790 
SU877800 
•SUB77810 
SUB  7 782 U 
SUB  77830 


34 


fcfHI TE(6,12>  CWiIH(  I ) »  I=1,H0HIT  ) 

12  F0RMAT<lflBU5X,51HTHE  SPANWISE  COLLOCATION  STATION! S)  OMITTED 
1ING#9I5> 

RETURN 


SUB77800 
ON  WSUB77850 
SUB778<5« 
SUB77870 


*'»*»  ****«»♦*  *«••*  •**•#*•*••••*»**«*•••*••••••»**••«••*•****••**  **»SllB7788ti 


20c  HCX  =  NWCX  S0877890 

NIX  =  HU  IX  «  SUB7790U 

DO  150  NS=1»2  SUB7791U 

OKI  TE( 6,22), f SURF .61',  NS),  1=1,2)  SUB77920 

22  FlJRBAT(lHl#31X,42HHISSILE  SUBSONIC  AIRLuaDS  PROGRAM  (C0NT-P)/1H  /  SHB7793U 

1  25X.39HCULi.0nA  1 1  ON  STATION  COORDINATES  ON  THE  2A6/1H0.12H  S  STSU87794U 
2A  N0,7X>2HYC.8X»7X,11:HXC  VALUES—)  SIJB7795U 

00  123,  IY=1,NIY  SUB77960 

YC  =  NB0«SH*Y(1Y)  .  ,  SUB/797U 

DO  120  IX=1,NCX  SUB/798IT 

120  XPRIIX)  =  HbO«XS(l»NS#,IX»  I Y )  SU87799Q 

123  WHI1£(6,124>  I Y. YC, ( XPRCI X ) . !X=3 ,NCX>  SUB78Q0U 

124  FORMAT (1H0, II?, 5E17.6/(1H  , 29X,4E17«6 ) )  SU8/8010 

RRITE(6,lUb)  (SURF ( 1 3  NS ) , I =1 »  2 )  SUB/8029 

105  FORMA T(1HU>24x,39H?NT£GRATI0NSTAIi ON  COORDINATES  ON  THE  2A6/1HG.  SU8/8030 

112H  S  STA  003.7X32HYl38X,7X3llHXI  VALUES—)  SUB78040 

DO  106  1Y=1,N!Y  SUB78050 

Y I  -  HBO»SN*FTA( IY)  SUB/8060 

DO  126  IX=1,NPX  SU8/8070 

126  XPRMX)  =  WBO*XS( 2#  NS.cXX.  I Y )  SUB/8080 

106  Wkl IE(6,124 )  fYVYi»(XPR( IX), IX=1,NIX)  SU8/8090 

NCX  =  NCCX  SUB7810# 

NIX  =  NCIX  SUB7811U 

150’  CON  I  I NUE  SUB/8120 

RETURN  SUB/8130 

»»•••«•*•*•*•*•**••« *••**••••*•«•*•*••**••*••**»»•**••*•*«**•** *#*SUB 7 8140 

30  no  34  NS  =  1,2  SUB78l‘iO 

MR  I  IE (6, 21)1  RtO( IKR),NMUnES,EK»EH  SUB/8160 

21  FORMATS  1H1  *c31  X,42HMISSILE  SUBSONIC  AIRLlADS  PROGRAM  (C0NT^D)//1H  /SUB/8170 
1  9X, 27HOSC ILl  ATURY  I  RbOllENCY  (CPS)»F12.b,13X,I2,17H  DEFLECTION  MODSUB78180 
2ES/1HU,8X,3HHRHJUCEI1  FREQUENCY  (SEMI  CHORD )* F9* b,14X, 23MFREE  STREASU8/8190 

3M  MACH  NUMBER. F9.3/1H  )  SUB7820U 

HR  I  I E ( 6, 31 )  T  MOD  SUB78210 

31  FURMAT(34X,34HPRESSURh  COEFFICIENTS  FOR  MODE  NO. 1 3//19X, 1H 1 11X. 10HSUB78220 

1R  CUEFF ( I )12X» 1  OH  I  COhFF(I)  9X, 9HSPAN  MuDE  3X,  KlHCHURD  MODE*)  SUB7823C 

WM 1 1 E ( 6., 32  )  (  SURF ( K 1 ,NS ) , KJ  =1  .2)  SUB/8240 

32  FORMA T ( 1H0 » 9X»  2A6// )  SU878250 


vb 

GO 

10.(2#  3)  ,  NS 

SUB78260. 

?  NL 

=  NWPX 

SU8/827  0 

ML 

=  NPY 

SU878280 

i 

i=K 

=  1 

SU8/829.0 

GO 

TO  4- 

SlJB/8300 

T 

3  NL 

=  NGPX 

SUB/8310 

1 

ML 

=  NPY 

Sua/8320 

* 

IK 

=  NWBX*NPYM 

SU878330 

4  DO  6  IM  =  1 » Ml 
00  6  IN  =  1 » Nl. 

WRITE (6, 33) IK,ANM< IK, I  MOD)# IM, IN 

33  FORMA T(1H0,  1 fc9 , 1P2E22 . 5 ,21 13  ) 
ft  IK  *  IK  ♦  1 

34  CUN  1 1 NUE 
RETURN 


SUB/834U 
>SUB/8350 
SUB/8360 
SUB/, 837  U 
SUB,/ 838  0 
SUB/ 8390 
SUB78400 


#**#♦»#**###»*# **##**#■**#*###**##**#*«*#**##*  ##*»##»* *«##*#**SU 8 784lo 


40  UR  1. 1 E  ( 6 . 41 ) 

41  FORMAK  1  HO  #  20  X»  38HERR0R 
1  TERM  1  NAT  I  ON  OF  CASE) 


SUB78420 

1NPU1  DA  1  A  (f^O  TAIL)  REGU I  RES// , 21X , 19HSUB  /843U 
35  SUB78440 


CALL  EXIT 


5 fl  WHHE(6,5l) 

51  FORMAT.<1HO,20X,63HNUH8EP  OF  COLLOCATION  OR  INTEGRATION  STATIONS 

1  PRESSURE  TF  RMS/Z21X, 25HEXCEEDS  ALLOWABLE  HAX1HUN///35X,18HCASE 

2  IERHINATED) 

CALL  EXIT 


SUB/8450 
►**SUB7846U 
SUS'/ 8,470 
0RSUB7848L* 
lSSUB/8490 
SUB/8500 
SUB78510 
^♦SUB/8520 


60  WR 1 1 E  ( 6»  21  )f  RfQCIFR),  NHODES*  EK<i  EH 
KHnE(6,i01)  SUHFt  1  * NSURF i 
101  FURMAT<1H0,27X,45HTNPUT  HODfc  SHAPE 
1//22X,62HREI FkENCED  TO  THE  SURFACE 
2C1I0N  /'/2X,4HHOD£*  2UX,  7HC0EFFS.) 

BO  69:  L=1  ,  NMOnES. 

NTH  =  NZCU(L, NSURF ) 

WK  I  I E( 6* 66  )  L » ( C0( L»  K, NSURF )  #'K=1,  NTH) 
FORMAT  Cl  HO,  14 . 4X» iP7El3»4/(9X, 1P7£13,4 ) l 
CUNUNUE 
RETURN 


66 

69 


SUB78530 

SUB78540 

POLYNOMIAL  COEFFICIENTS  FOR  ,A6SUB/8550 
LEADING  EDGE-CENTERLINE  INTEKSESUB78560 

SUB78570 
SUB78580 
SUB/859b 
SUB786Clb 
SUB78610 
Si;B78620 
SUB/863U 

f#*»********************SUB/8640 


70  WK M  E  ( 6» 21  )KR*;0 (  I>  R ) , NHODES,  EK,  EH 
WRI 1 E ( 6 » 35 )  ALPHA (I MOO ) 

35  .FORMA  J(24X, 47 HRMS  ERROR  OF  Duk'NWASHES  A1  COLLOCATION  POINTS 
X  .  6  ) 

HR  II E  ( 6, 36  )  I  HOI) 

36  FlJRHAT(lHQ*23x. 57  H  PRESSURES  AND  UPWASHES  AT  COLLOCATION 
X  MODE  NU • 13) 

WR  HE  <6,32  )(  SURF  (L,  NSURF  ),L=1,2) 

WK 1  I E ( 6 , 37 ) 

37  FURMAf(THU,/X,lHX,8X,  LHY,9X,8HR  P (  X*  Y )  ,  5X»  8HI,  P(X,Y),6X, 

1  9KR  CM( X, Y ) , 4  X, 9H I  Cw( X »Y ) , 6 X, 9HH  UW( X, Y ) » 4X»9H!  DR(X,Y>) 
RETURN 


2/  ) 

DO  78  I=l,NMnilES 
OH  78  J=1 ,  NMOIJFS, 

IF (NSURF. EO. 3)  GO  TO  76 
01  =  REAL ( GFONC ( I  *  J* NSURF ) ) 

02  =  AIHAG(GFURC( 1,J, NSURF)) 

OU  10  77 

76  01  =  REAL  (GKORC<  ?,  J.l ) )  +  REAL  <  G.FORCT  I  ,  J,  2  )=) 

02  -  A  1  MAG  (  OF I1RC  (  I ,  J,  1 )  )  ♦  A  I  MAG ( GF ORC ( I , J, 2  )  ) 
KKK=2*NM0DES 
NNN=2#J-1 
N  N  N  N  =-2  *  J 
CaRI)S(  I,NNN)=G1 
CARDS! I , NNNN ) =02 
n  03  =  SORT <61 **2+G2*»2 ) 

04  =  0.0 

IF  (G3.NE.U.U  )  04  =  57 . 2957795*  A  TAN2  ( G2 ,  Gl^) 

WR  I  I E ( 6 , 71  )  J , J, 01, 02,03, 04 
IF  (NSURF  .NfV.  3)  00  ?()  7  8 
IF  (I  .Nr.  NMODES)  GO  TO  78 
IF  (J  .NE.  UHUDES)  GO  TO  78 

PUNCH  6969,  < ( CARDS < I  I , J J ) , J J  =  1 , KKK ) ,  jl=l,NMODES) 
6969  FURMAI  (1P6P12.5) 

71,  F0RMAH1HU,  IQ,  I7.2X,1P3E19.5,0PF16.3,4H  DEG) 


SUB78650 
SUR7866C 
=, 1E13SUB7867U 
SUB78680 
SUB7869U 
P-UINTS  F ORSUB7870  0 
SU878710 
SUB  787.20 
SUB 78730 
SUB78740 
SUB7875U 
SUB7876U 


<mm<  *##*#*********e***##****#**»*********SUB7877Q 

80  WK  J  I E  ( 6 , 21 )  F  R|:0  <  TF'R  ),NMUDES,EK,EM  SUB78780 

WKIIF  <6,61J  <  SURF  (L,  NSURF  ),L=.l,  2.)  SUB78790 

61  FORMAT! 3 5X*23H GENERAL 1ZEU  FORCES  FOR  2A*/3 HO , 6X, 4HDEFL, 3X, 4HL0AD, 1SUB7880 0 
10  X,  9HRF-AL  PART  *  10  X,  9H I  MAO  PART  ,  10  X ,  9HABS  VALUE,  10X,  11HPHASE  ANGLE/SUB/881U 


SUB78820 
SUB 78830 
SUB/8840 
SUB  7885,0 
SUB78860 
SUB 7 88 70 
SUB7888D 
SUB788>6 
SUB 7 8900 

##*#****. 

i 

**»*#*»# 

##*###»# 

•*#*#*##« 

SUB789101 

SUB/8920 

SUB78930 

SUB7894U 

######## 

##**»**# 

*•**•«** 

SUB78950 
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78  CONTINUE  -  SUB/.8960 

RETURN  SU878970 

C  ••*••••••••••*••• ■  *«*»*»•«•«.»***«•*•<-**•••*•••«»•***«•«•****•#•  4  #*SU8  7  898  U 

9*1  HRI1E66,19)  XCOLl.*  YCOLL*PH(  IRON*  IHOD}*C«ASH(  iROH* 1HGB)»  SUB78990 

X  RWASH(IROW>IHOD)  SUB79000 

19  FURhAT{lHQ,?X,2F9.3.2X,2E13.4,2X,2E13.4,2X,2El3;4)  SUB7901U 

RETURN  ‘  t  SUB79020 

END  SUB79030 


i 
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4.0  TRANSONIC  UNSTEADY  AERODYNAMICS  PROGRAM 
4.1  Theoretical  Derivation 

Whan  the  flight  speed  approaches  H-  1.0,  the  velocity  potential  equation 
can  be  written  as  . 


$  .  +  «J>  «  M2  (2ik$  -  k2<*>) 

yy  Yzz  y  x  ~ 

which  is  valid  if  k»|M-lj.  The  linearized  equation  is  applicable  when  the 
lifting  surfaces  are  oscillating  rapidly  such  that  non-linear  disturbances  in 
the  flow  do  not  have  time  to  accumulate. 

Equation  (  4.1.1)  is  satisfied  by  a  pulsating  doublet  which  produces  a 
velocity  potential  at  (x,y,z)  given  by 


1  k  (z-g) 

2  (x-O* 


exp 


-1/2  ik 


[<- 


0  +  (y-n  )2  +  Cz-d2 
(x-?)  - 


where  the  doublet,  is  positioned  at  (5,n,£)'.  The  doublet  has  no  influence 
at  points  upstream  of  the  line  x  *  £  and,  consequently,  the  potential  is 
zero  in  that’  region. 


A  solution  to\  equation  (4.1.1)  may  be  obtained  by  superposition.  This 
solution  will  be  represented  in  the  form. 


<Kx»y,z)  « 


<KS,n,)$D(xiy,z,5,q,'o)d5  dn 


(4 


where  <{>(£, q  )  is  the  dpublet  strength  of  point  (£,n). 

To  compute  the  velocity  potential  distribution,  the  wing,  wake  and 
control  surface  is  divided  into  a  lattice  of  square  boxes  as  shown  in  Figure  4.1.1 
The  potential  function1  is  replaced  by  a  set  of  point  values  at  the  box  centers. 

The  potential  functions  and  downwash  value  is  assumed  constant  over  each  box  and 
equal  to  the  central  value. 

The  problem  reduces  to  imposing  boundary  conditions  and  determining  the 
doublet  strength  for  each  box  to  satisfy  the  boundary  conditions.  The  boundary 
value  problem  becomes 


1.  Tangential  Flow  Condition 


w(x,y) 


wing 


f  (£,  n)lim  d£dn 
z+o  3  z 


(4.: 


control 

surface 


•tt 


w<x,y) 


//. 

wing  + 

wake  + 

control 

surface 


<K£»n)lira 
z+o 


d£dn 

cjz 


(A. 1.5) 


2.  Zero  Pressure  Jump  in  Wake 


—  ik$  =  0  '  (4;li6) 

ox 

Equation  &.X.6  )  is  an  ordinary  homogeneous  differential  equation 
subjected  to  the  condition  $  is  equal  to  the  value  of  the  velocity  potential 
at  the  wing  trailing  edge  for  x  =  x  ^  t  e  *  This,  gives  the  solution 


\ 


^wake  ^wing  t.e. 


-ik(x-x  ,  .  ) 

e  wing  t.e. 


(4.1.7) 


Equations  (4.1.4),  and  (  4.1.5)  and  ("4. -1.7)  form;  a  system  of  equations 
from  which  the  point  values  of  the  potential  functions,  can  be  found  at  the 
box  centers.  The  pressure  distribution  is  then  determined  from  ;the  relationship 


AP(x,y)  =  -fpu2C^+  2ik$) 


(4.1.8) 


and  the  generalized  forces  are  found  from  equation  (2.0-1) . 


44  PROGRAM  DESCRIPTION 


.  The  Sonic  Unsteady  Aerodynamic^  Program  calculates  generalized 
forces  for  up  to  10  deformation  modes .  The  computer  solution  is  based 
upon  the  Mach  box  technique.  The  various-  configurations  which  can  be 
analyzed  are  shown  in  Figure  4.2.1  and  Table  4.2.2.  The  analysis  includes 
interaction  effects  between  tandem  surfaces  and  wake  effects  oh  the 
trailing  surface.  Single  surfaces  may  be  analyzed  by  imputing  a  second 
surface  with  a  zero  chord  length. 

the  transonic  box  method  calculates  the  unsteady  potentials  from, 
which  the  pressure  distributions  may  be  obtained  for  arbitrary  modes  of 
surface  motion.  The  method  used' was  suggested  by  the  successes  of  the  super¬ 
sonic  box  methods  of  Pines  and  others.  Reference  il.  The  potential;  is  gener¬ 
ated  by  a  doublet  distribution  rather  than  by  a  source  distribution 
because  the  latter  method  would  involve  diaphragm  regions  of  infinite 
extent,  whereas  the  doublet  distribution  is  confined  to  the  wing  and  its 
wake.  As  with  the  subsonic  problem,  the  differential  equation  solution  is 
an  integral  equation.  The  integral  equation  is  approximated  numerically 
by  a  matrix  equation  so  that  the  basic  step  in  the  box  method  is  the  solu¬ 
tion  of  the  system  of  simultaneous  equations  which  determine  a  set  of  values 
of  potential  bn  the  surface  from  a  corresponding  array  of  upwash  values. 

The  solution  procedure  obtains  the  velocity  potential  over  the  surface  one 
spanwise  row  of  boxes  at  a  time  until  the  trailing  edge  row  is  completed. 

The  numerical  complexity  is  hot  increased,  however,,  by  a  large  number  of 
box  rows  over  the  configuration  because  the  influence  coming  from  mor6 
than  15  rows  away  has  been -found  to  be  negligible.  .  The  results  are  valid 
for  high  reduced  frequency,  k,  -such  that  k  ))1  M  -1  j  where  M  is  the  Mach 
number,-,  I  I 


The  solution  for  the  generalized  aerodynamic  forces  requires  the 
input,  of  the  deformation  modes  due  to  vibration.  The  program-  considers 
the  modes  4va  tie  expressed' as  analytic  functions  of  the  form: 


w(x,  y) 


N 

-£ 

m=0 


n 

m=0 


C 


(n-m)., 


(n-m),  m 
x  y 

,m  J 


To  meet  this  requirement  only  the  coefficients  nc"  are  required  as  input 
into  the  program.  These  coefficients  can  be  obtained  in  several  ways, 
the  most  common  way  is  to  surface  fit  the  modes  by  the  least-square 
technique. 
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TABLE  4.2.2  OPTIONAL  CONFIGURATIONS 


configuration' 


RECTANGULAR 


CKORDWISE 

COORDINATE 


SPANWISS 

COORDINATE 


i(I)  *  X).0 
X(2)  =  0.0 
X(3)  >  0.0 


Y(l)  =0.0 
Y(2)  =  OoO 
Y(3)  >  0.0 


X(4)  >  X(3) 
X(5)  ?  X(4) 


■ 

X(l)  =  0.6 

X(2)  >0.0 

Y(l)  =0.0 

Y.(2)  =0.0 

DELTA 

X(3)  =  X)2) 

X(4>  £  X(3) 

X(5)  5  X(4) 

Y(3)  >0.0 

i 

X(l)  ■  0.0 

X(2)  >  0.0 

Y(l)  =  0,0 

Y(2)  >  0.0 

TRAPEZOIDAL 

X(3)  =  X(2) 

Y(3)  >  Y(2) 

s 

X(4)  5  X(3) 

X(5)  3  X(4) 

, 

TRAPEZOIDAL  (CROPPED); 

X(l)  =  0.0 

X(2)  >X(1) 

X(3)  >  X(2) 

Y(l)  =  0.0 

Y(2)  >  0.0 

Y,<3)  >  Y(2) 

- 

X(4)  >  X(3) 

X(5)  iX(4) 

• 

X(l)  =  0.0 

Y(T):  =  0.0 

DELTA  (CROPPED) 

X(2)  >  0.0 

X(3)  >  X.(2) 

Y(2)  =  0.0 

.  X(3)  >Y(2) 

X(4)  5  X(3)  ; 

*  '  X(5)  *  X(4) 
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4.3  INPUT  INSTRUCTIONS 

Instructions  for  preparing  input  data  for  the  transonic  . 
computer  program  are  presented  here.  The  field  location  and  format 
for  each  quantity  is  specified.  Any  set  of  units  may  be  used  for 
geometric  dimensions  and  acoustic  velocity  as  long  as  they  are  consistent, 
e<g.,  if  inches  is  used  for  length,  then  the  aeons tic  velocity  must 
have  dimensions  of  inches  per  second.  The  required  data  and  the  sequence 
in  which  the  information  is  entered  is  as  follows: 


1.  Scvearawise  Coordinates  (6E12.5  format) 

Column 
Name 
Item 

(1)  X(l)  King  root  leading  edge  coordinate  (See  Figure  4.3.1) 

(2)  X(2)‘  Wing  tip  leading;  edge  coordinate 

(3)  X(3)  King  trailing  edge  coordinate 

(4)  X(4>  Control  surface  leading  edge  coordinate 

(5)  X(5)  Control  surface  trailing;  edge  coordinate 

A  single  surface,  the  wing,  may  ;be  analyzed  by  setting  X(4)  and  X(5) 
equal  to  X(3) .  -The  various  configurations  are  generated  as  shown  in 
Table  4.2^1.  lT«e  origin  for  the  planrorm  and  AIC  station  coordinates 
rausL  be  at  the  leading  edge  root  of  the  wing,  therefore  X(l)  and  Y.(l) , 
described  below,  must  always  be  zero. 


Spanwise  Coordinates  and -Acoustic  Velocity  (6E12.5  format) 


Co] umn 

>12  ! 

13-24 

25-36 

,  37-48 

:  49-60 

61-72.  1 

Name 

L— - - : - ' — 

Y(l)  ! 

Y(2)  '  ' 

X(3) 

SOUND 

l  RHO 

_Q_I - - - 

! 

i 

Item 

a)  ! 

(2) 

(3) 

!  (5)'  . 

i 

j 

(1)  Y(l)  Wing  root  spanwise  coordinate 

(2)  Y(2)  Wing  leading  edge  spanwise  coordinate 

(3)  Y(3)  King  (and  control  surface)  tip  spanwise  coordinate 

(4)  SOUND  Speed  of  sound  at  altitude  for  which  analysis  is  performed 

(5)  RHO  density  of  fluid  *  1000 ..0  C-M./L^) 
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3.  General  Information  (6112  format) 

Column 
Name 
Item 

(1)  NMACH  Number  of  Mach  numbers  (max  6)- 

(2)  NFREQ  dumber  of  input  frequencies  (max  10) 

(3)  NM0DES  Number  of  input  modes  (max.  10) 

(4)  HBW  Number  of  chordwise  wing  boxes  (max  10 ) 

(5)  LVPIC  Print  velocity  potential  influence  coefficients;  0 

'l  Yes 

(6)  LSSVP  Print  upwashes;  0  ^No,  1  ''-'Yes 


4.  Mach  Numbers  (6E12.5  format)’ 

Column 
Name 
Item 

(1)  FMACH(i)  Mach  numbers  fGr  which  the  analysis  is  to  be  performed. 


5.  Frequency  (6E12.5  format) 

Colusa 
Name 
Item 

(1)  FREQ  Frequencies  for  which  the  analysis  is  to  be  performed. 

Continue  on  next  card  for  FREQ(i)  6. 

6.  Deformation  Modes.  Repeat  the  Following  Cards  NM0DES  Times 


Column 

Name 

Item 


(1)  NTMl(i)  Number  of  deformation  mode  coefficients  for  the  wing,  mode  i 

(2)  NET  Compute  generalize  forces;  Q^No,  1^  Yes 


If  NFI  =  0  the  program  will  compute  the  VPIC's  ahd  stop. 
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Column 

1-12 

jj.  ruimai, 

.  13-24  . 

23-36 

37-48 

49-60 

61-72 1 

Name 

C0  (2)  _ 

00(3) 

.  00(4)  _ 

C0(5)._ 

IK?!R3l 

Item 

w-.V 

l.-iai _ 

-  .  (3)  .. 

(6)  ! 

(1)  C0(i)  i  =  1,  HIM  deformation  polynomial  coefficients  to  be  input 
in  the  following  order:  0,0;  1,0;  0,1;  2,0;  1,1;  0,2; 

3,0;  where  the  first  integer  is  the  power  of  "x"  and  the 
second  is  the  power  of  "y".  Continue  on  successive  cards 
until  all  polynomial  coefficients  are  input. 


Column; 

Name 

Item 

(1)  NTM2(i) 

(2)  NFI 


j  1-12 

13-24  . 

25-36 

37-48 

49-60 

i  "  — 

1  NTM 

.  NEE  '  . 

• 

} 

(2> 

Number  of  deformatio-  mode  coefficients  for  the  control 
surface,  mode  (1) 

Compute  generalized  forces;  0 ''-'No;  l '■'-'Yes.  If  NFI  =  Q, 
the  program  will  compute  the  VPIC’s  and  stop. 


Column 

1-12  * 

13-24 

25-36  ~ 

37-48 

|  49-60 

Name 

caci) 

C0(2) 

C0(3) 

C0C4) , 

C0  (5) 

Item 

a) 

(2) 

(3) 

..  (4) 

(5) 

C0(i)  i  -  1,  NTH  deformation  polynomial  coefficients  to  be  input 

in  the  following  order:  0,0;  1,Q;  0,1;  2,0;  1,1;  0,2;  3,0; 
"  etc.  where  the  first  integer  is  the  power  of  "x1'  and  the 

second  is  the  power  of  "y".  Continue  on  successive  cards 
until  all  polynomial  coefficients  are  input. 
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SAMPLE  "PROBLEM 


The  generalized  forces  are  calculated  for  the  configuration  below. 
The  flight  parameters  and  pertinent  inpiut  data  are  presented  on  the  first 
page  of  the  computer  print  out.  *' 

The  coefficients  of  the'  deformation  modes  for  the  forward  surface, 
are  shown  on  Che  third  page  of  the  computer  print  put,  and  for  the  aft  surface 
on  the  fifth  page  of  the  computer  print  out. 
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C0HH0N/C7/C7F»0,P0l  »  PHI  f  E/ DPHI  ,SPHI<  SDN7009U' 

COHiiON/CB/Rl.O  ,  SON70100 

lift!  a  CZI-RO/C  0.n,«.0  )/, TP  1/6. 2831 853/ /F  N/0 .  /  1.  ,  0/.  / 2. / 1. ,  IF.  , 3 .  »  2.  /  SUN/0110 

1  1  .  #  (l>i  4»/3..2.«l.  >11  ./5a/4*/3. *2.-/1./  0  •  /  6  .»  5  •  *  4  .  /  3  •  >  2  •  *  1  *  *  0  ._/  SUN/Ul  ?U 
POO  SUN70130 

FPPTnAH  NL'STHII,  O-CK  S0N/lli4U 

HA  hi  SON7016I) 

COMPLEX  V!‘Ii.,i«S,n-J,U,PHlW,CK,CZbRO,PHl,PHITfc,DPH|,SPHI>ASU,&XF  S0N7U170 
m*JlNSloN  A!«(J  ( 4!i  *  4  0  )  SON/018U 

COOMlM/t  l /K:  nx(xil;iU),xF(5  >,  YE  ( 3  ) ,  AR  ( 3  )  *  XI*  X2,  X3 .  X4 *  Y1/.Y2, 8FT  A ,  NBS  SON? 0190 
r.OKh0N/C?/A-..:H.MMi.F  4ACH ( 6 ) , NFRl Q, FREO < 10  )  /  NMODE*  NSURF  *  LVP1  C>  LSS  YPSON702.U  U 
i:0MS.0M/i;3/n-w  n,2«i^?>,Nl  <  L0,2,>iirF(lG/i>/NTMAX(2)/FNl28)/DXE(6):,  IPISON7021U 
nmr  ON/(:4/Viii:(HU,i,j>/.bsi2!iu0),l)Q(28/2}/G(10,l0/3>/PHlW(5U)/CK(4  0)SON/64f>.u 
ClikHP«!/i;t>/Ks>»(l*in  I*||RL-QOU),FO  *  I  FR>  XL*  NJS*  N I H*  K*  QR»  01  *  0A8  *  QAN  SON/ 0  230 

MI.,hIH/1.6/X,  Y,IU,  0Y,  FH,tK,EK8/.EKR/NP,KP,Ntf,oB0X/K0Dfc,MbijE/N8N/NBr  SUN/ 024 u 
piihr  lM/C7/C/f:<0,p.<l  ,P.li]c.J)PHl/SPHi  SON70250 

COhlU'N/l  8/Hi-i)  SUN/0260 

CAM  oAlti  SON  7  027  U 

{•it  l(liil)  HA.Hi.si ,  SON/U280 

PH  =  H'A'ChC >  ACH)  SON/ 0290 

inARS(fcH-J.Q)«isT.!(i.9t>>  GO  10  100U  SON/03UU 

CA!  I  COliF  SON 7 0310 

CAM  i’llUT  { 1 )  *  SON  /  0320 

TPII=1HI/(AS*F  D  SUN7033U 

pro  =  ir-x  .  son/ 0340 

110  you  If  |J  =1  ,H\  RM1  SON/U350 

FKM  KFO(  1FR).-1P.<  '  SON/036.!J 

!  F  O  K.  f:0«  0  .  »•  )•  C"  Id  900  SON71P37I1 

FM’-  =FK*«F  ‘  S0N/U38U 

FK»  =  FK»XJ/5»il>  SON/0390 

CAM  POT2h  *  SON7040U 

lFnVPlC.NF.il)  CAM  PlHJl  (?)  '  .  SON? 0410 

AMl-FK*lrX  S0N71I42II 

F.XF  rC-lPLX  (  Ci'S  (  AkC  )  >-S I N<  AHO ) )  SON/ 0  430 

1)11  t» 0 il  HUI)F  =  1  ,  N-AO’ip  SUN/II44U 

HO  Hi  J  =  l/5'  SON/ 04511 

1)0 1  J>  l  ):=C2F»<0  SON/046U 

irnif  (HOHI-,  I  ).F-J.-i)  nil  TO  210  SON/11470 

x = n  i  *>  *  r«  x  son/  iMBo 

Nil  -  J  SON/U490 

!F(ISSVP.uF.f|)  CAI.L  PuNT(3)  SON/U5UO 

no  200  NPrl.N'iOX  SON/061U 

Pil  =  HPH(NP)  Suft/0520 

Y  =  0.f)  SUN/0530 

KliOF  =  KNOX (  H*»  >  SON/O'HU 

NS  =1  SON/055U 

nil  10  (1?»M»J2»i  i»M»120-)»hOlJF.  SDN/0560 

NS  =2  SON/U57U 

00  Ml  MP  =  1  /  R  SON/IJbflU 

SPM  =  CZhRP  ,  SON/0600 

I  F  <  i'P.  0  I « 1  )  CAM.  PFIIK  SON7  0  61 0 

r  A I  I  NASH  S0N/062U 

CKd’P)  =  0b(N-  )  SON/0  630 

os •  up )  =  iis.di,  )  -  sp»ii  Sun /ii 64 u 

Y  =  Y+liY  SUN?  U65U 
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2li  Kli  =  Nft*l 
NO  =  NH-Mk 
no  30  1 11  =  3  , 

no  ao  ju=ii*8 

J  JO  =  IA8S{  I  O-.l  J>»1‘ 

25  AS0(I0,JGJ  =  VPlC(lJO,l) 

I F  c  jQ.FO.i  )  on  1 6  3 n 

fJO=I(j  +  JO-l  1 

ASO(  l«J,  JO)  =  /“i'>(  lu.  Jf;)  +  VPlU(  1  JU,IJ 
an  CONTINUE 

ISO  =  MS  l  ME»'(40  >  h»>i  3  .  ASO»  l)S(  MB  )-> 

IF  (I  SO* F 0.1,)  no  T0>  39 
CAI  l  POUT (8  ) 

-no  io  90 u 
0‘J  CON  11  NOE 
Y  -  0.0 

If(liP.9fc.l  )  On  1 1)  5 H 
HO  -|5  MI?=l,~R 

45  nsf NP)  =  llS(rt»*-)»?  KV3. 1415927 
5.1  DO.  I  I  OIJF 

IF  (KO0fc. «E. 4  )  >H)  TO  80 
nu  on  hp=i,fr 

ns <  1-.P  )  =  PFI I  W ( H*’  ) R  FS( NH ) -Ph IN ( HP ) ) *2 . 0>'.3 .1 415927 

OH  N'Hrf‘R  +  1 
NH=NR-H8 
«■!  COM?  HUE 

00  jno  HP=i,sn 

If  (KOIJE.EO.3)  PHlri(.MP)=llS{NH)*EXF 
IF  (op.Ffi.NBOX-l  )  -’H I  U { Mf* )  =nS ( H8  ) 

Pll  ME  =  US  bin 

tn^P.FlJ.NHOX)  MHlff  =  PHI  Tt*  (Plil  lE-PI;Ih(HP)  )«DXE(5.)> 
PHI  =  ns (MB) 

f r f i  ssvp.ne.o  >  car.  poiikV) 

CAM  D01J 
NH  =  NH  +  1 
1  U I)  Y  =  Y  4  n  Y 
no  io  .20  0 

1  2n  no  1.3 u  MP  =  1 .  h>< 
nSf  f*fi)sphlH(MI») 

PH  I  K (  HP  )  =  F  XMH’Hi  W(  H|») 

CM  fr'P  >  =-CZlP>* 

If  (I  SSVP.OE.fi  /  CA'lL  Pull, I  ( 4  ) 

lan  n«  -  nb+1 

2U.I  x  r  X  +  0X< 

710  no  4  0  0  HO  =J  , '(MOIjl: 
no  30  0;  NS=1 ,  NSU^F 
C  (  MODE  ,  hO,  NS  )=C2EKO 
NTM  =  N  T  ( MO,  N.S ) 
no  300  H=i;  M 

3UI1  R{Honp,HO,N.:.)=(Hrlonh»MO,NS)+CU(HO»N,NS)«DO(N,NS) 

C  ( Munp  « Mfi,  3 )  =  C75RO 

4  0  0  Q(l:<inE>M0,3)=i)('if..iF,Hu,  j  )  ♦()  (  MODE ,  MO,  2 ) 

5  U  i|  COMIHI-I 

00  HO  0-  N S  =  1 , 31 

CAl  1.  POO  1(6) 

no  700  J-l,.‘MOlH- 

no  70  0  K  =  l,  .HOlit 

Ql<  =  F  0*  a  I  (.JFK,  J,FS)  ) 

0  I  =  R,  *  A  !  M  <"  n  (  0  <  M  J ,  N  S  )  ) 

C  AP  =  S0F  i  (  01  sQ  [c+  JK*OP) 
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SUN7066U 

SON/U67U 

SON70680 

SON/ 0690 

SON/O 700 

SDN70710 

SON7072U 

SON/D730 

SON/0740> 

SON 7 0750 

SON70760 

SON70770 

S0N7U78U 

SOf^/ 0)790 

SON/OHOO' 

SON/O&lO. 

SUN78O20 

SON/O830^ 

SON/064U 

SON>I)85U.: 

SuM/0860 

SoN70H7.lI- 

SUN40880 1 

SUN70b90 

SON/ 0900 

SUN/ 0.9-10 

SON/0920 

SON/094U 

SON7095U 

SCN/0960 ; 

S0N/097U; 

SON/0980 

SON/099U 

Son  /  moo 

SON/l OlU 

SON/1026 

S  f)N  / 1  il  3  0 

StlN/1  li  4  U 

SON/1060  o 

SON/1U70 

SON71080 

§0N/li)9U 

SON/13  00 

SUN'/lillJ 

SUN/1120 ■ 

SUN/1130 

SUN/1 140 

S  0  N  71 1 5  0 

SUN/ 11 60 

SON  / 1.17  U 

SUN/1180 

SON71190 

S0N7120U 

SON/12TO 

S0N/1220 

SON 71230 

S0N/1740 

SON/ 1250 

SON/126U 

son /l 27u 


i 


7Uf| 

80  0 
90(1 

1  non 


r,ii£  i  f 


in 


2n 


3  ii 

l  i 

l '? 


r.rau- 


oAfi=o.n 

1 t  (OAB.NE.O.O)  JAN=57.29578*AlAfi2(QI»CR> 
CAI  L  P0UK7) 

I F f  MSURF.EQ.1  )  CD  T(J  900 

COM-1  HUE 

CUM1INUE 

cominuf  - 

fid  TO  1 
FMI* 

M»RTKAN  NLSTnO,  D*-Ch 
DA  in 


SON7.1280 

S0N/1290 

SON71300 

S0N7131U 

SON713.20 

SON71330 

SOH/1340 

S0N71350 

S0^/1360 

SO 

SON/ 1390 


SORkOllTINE  i-AIN 

COHHnN/Cl/K,«fiX</UuO  ),XE<-5),  YE(3),AR(3),X1,X2,X3, 
C»mhON/C2/AS.  <MAC'I,  f  HACK ( 6  )  »'NF REQ,  FREQ MQ  )  ,  RHODE 
COMf  On/i;3/C- ( 10 ,23,2) , NT ( 10 ,? ) ,.NF  ( 10 , 2  >=,  NJ  H  AX  (2.) 
CUMO0N/C6/X,  Y  ,||  *,  :iY,  Frti EK , EK8  *  EKR »  NP,MP,NB,nBOX» 
CflH«ON/V.8/Ri*0 
R1AM5, 11)  (  Xl  f  1  )  .1  =  1,5) 

UFA  1.(0,  J  1 )  {  Yh  (  I  ),.  1=1,3),  AS,  RHO 
H*l(i  =  Ri|Cl/l .» 0 .1 .  ' 

HI  Ai.(*3,T2)  --Mac  » » “(FI'FO,  NMODE,  NftR ,  I.  VP  I.C rLSS VP 
*1 1:  A !•  <  *3 , 1 1 )  (  FhAi.h(  I  ),  )  =1,  NHACH  ) 

II  l:  A  f*  ( 5 ,  J'i  )  (FRF  KT),  !  =  1,  NFRfcQ ) 

NSOPF  =  2 

1 1'  (  XF  (4  ) .  LT»  XF  ( *>')  )  GO -TO  10 
nSORF =1 
<FM  )=X[.  (3  ) 

VI  ( *> )  =Xt  ( 3  ) 

N?fMAX(1.)  =  U 
N|KAX(2)=U 
DO  30  KOlJF.  =  l,..H  »Dl 
JiO  30  !  =1 ,  NMIRt 
no  20  J  =  l,2»- 


SON/140U’ 

X4,Y1,Y2,8ETA,NbS  SON7141U 
,NSURFiLVPIC,IU:SSVPSON/1426 
,FN(2U),DXE(6.),TPISON/1'130 
KODE,MODE,NBW,NBT  SON/1440 

SON/1450 
SON/1460 
Son/ 14 70 
SUN/146U 
S0N714.90 
S0N/150U 
S0N/161U 
SlIN/1520 
SON/153U 
SON71540> 
SON/1550 
SON71‘,60 
SON7167U 
SON  71 53 U 
SON/15,90 
SON7J60U 
SON  7 1,61  0 


00  (  MO  III- ,  J,  I  )=0  ./• 

HFAIK6,i2)  TM.'.FI 

M  f  MIIJF.,  I  )='  TM 

R  T-M  A  X  (  I)  =  MAXO(NIHAX(  I  ),NTH) 

Nf  <  MO Of ,3 )=• F! 

UFA |i(*,  11  )  (CO (MODE,  J,  I  ),  J  =  1,NTN) 
FORMA T( 6E12.8) 

FORMA  1(6112) 

RFTORN 


S(iN7l.o20 
SON71630 
SON/1640 
S0N7165U 
S0N71660 
S0N71670 
SON  71 080 
SON/1690 
SON 717 00 


r-Nf  S0N71  /i'O 

FOITKAN  oLSTiill,  (M.-Ch  SON71/20 

CODE  SON71740 

SliPPOII T  l=NF  l  Oli[.  SON71750 

COMMON/Cl/KI  Ox.(  2(HlO  ),  XF(5  ),  YF(3  ),  AR(3  ),  XI. ,  X2 ,  X3>,  X4 ,  Yl  ,„Y2 ,  BE  I  A ,  NOS  SON71760 
COMMON/C2/A%  0MACM,» MACH(6),NFREO,FhEO(10 ) , OMUDE, NSURF, L VPl C, LSSVPSON/ 177 0 
COMMON/ 03/C>  (  1 0,23, 2)',  Ml  (10,2),  NF(.1  0,2),  N ThAX  <  2  >  ,’FN  ( 28  > »  DXE  ( 6  ) ,  1  P 1  SON/178 0 
C0MM0N/C6/MiBU  nil  ),  NHL  (10 II  ),FO  ,  I F  R,  XL  ,  NS ,  NT  H,  K ,  J,  QR>,  U I ,  OAB ,  Cl  AN  S0N71790 
COMMON/ C6/X,  Y,  fi  x,,  l!Y ,  EM ,  1. K ,  EK B ,  FMI ,  NP  ,  MP ,  NO ,  NBOX ,  KODE ,  MOLlE ,  NB W ,  NB  I  SON71800 
COMMON/  C8/Rl>9  S  0  N  7 1 8 1 0 


HFTA 

=  EM 

SON71820 

XJ  = 

XF(3 )  -  X F ( J  ) 

SON/1830 

X2  = 

XF ( 3 )  -  <F(?) 

SON/184  0 

X  3  = 

XE(4)  -  X F ( 1 ) 

S0N71850 

It 

X 

XF ( 5 )  -  XF ( 4  j 

SON  71 861) 

XT 

II 

XF-  ( 5  )  -  X F  ( 1 ) 

SON/1870 

YJ  = 

YF (2)  -  Y F ( 1 ) 

SON/ lb80 

Y2  = 

YF(3  )  -  Y F ( 1 ) 

SON/ 1890 
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I F ( X2 . GT. XI. OR. X3 . GT . X3 . OK. X3 . GT. X5.0R.  X2<LT.;Q . 0  )  GO  to  50 
I F.(  Y1 . G,T.  Y2.0K.  Yl.LT  i  0 •  0  )  GO  TO  50. 

T  HI  =  0.0 

I F( Y2 . NF. Y1 )  TwL  =  (XI  -  X2)  /  (Y2  -  Y1 ) 

ARI1)  =  (  Y?«  (  X2  +  X3  3  -  Y1  *  (X2-X1M 

AR(2)  =  Y2«X4*?.II. 

AP(o)  =  -ABO  )  +  AR  ( 2 )  ,  - 

lit  OX  =  X1-/(FL*  A  TMBW)  -0.5) 

IF  (3’llG.O*  •<,  .  i>T .  x5 )  GO  TO  20 
15  NOW  =  NOW-1 
GO  IO  TO 
2r,  nr  =  iix/bota 
YOi  =  Yl/JlY 

YN2  =  Y2/0Y 

XNl  =  YH2  -  (cXl-X2)  /  nx 
XNT  =  YH2  +  XO/iJX 
XNI  I  =  X3/f)X 
Xim  =  X5/DX 
N»OX  =  XNTR  ♦  -.1  .b 
WHS  =  Y2/DY  *  1.0 

WOT  =  X4/IJX  ■*  II  .5 

.  nxRD  =  i.i* 

DXFI2)  =1.1* 

|lXF(3)  =  0.5 

DXF 14  )  =  A  I '  r  (  XhI.i*  ♦  1.5)  -  XNLfc 
nXF(5)  =  XNIF  -  riUAMNBGX-l) 

0*F(6)  =  II.!.  ■' 

X  =  U.5  *  Dx 
NO  =  (I 
sKUDE  =  1 

HO  ■'i  0  NP‘  =  1  ,  OB-JX 
XN  =  FLOAT  ( *  P  )  -il.5 
YW  =  YN? 

II  (TkL  -GT.  ii  . »» )  Y  W  =  AM  I  Nl\(  YNl  +  XK  /  (  ThL/BHT  A  )  ) 

Nil  =  f  F I  X  (  Y +  J 
28  MJiMNP)  =  N:- 

IFIMB.GJ.40 )  GO  fn  15 
If.  (NP  .:b‘U.  NOW)  KOl.F.  =3 
IF  (NSUKF  .10.1)  00  TO  29 
IT  (X  .GT.  xl  )  c.OiiE  -=6 
IF  (X  .GT.  X3  )  KuOfc  =4 
IF  IX  .GT.  X3H)*)  K(ini;  =  ? 

IJ  (NP  .FO.  NO  OX.-)  K  (Ml  B  =| 

2Q  I  F(,NR  +  NII.GT .  ?')  II  u  )  Gli  TO  15 
Mil  INP)  =  N- 
III)  30  MP  =  1,mR 
NH  =  NO  M 

J.!l  KHOX  (MB  )  =  M'.lF 
4(1  X  =  X+TTX 

OOPliO  =  0.5  M  AS*FM)«*2 
TO  =  -  Svifc*  I).X  *  ii  Y  »  0  ‘  i  R I  •  0/  F  h  *  R  H  li 
RETURN 

50  CAM.  EX'M 
RETURN 
FNH 

1  FORTRAN  NL 'I  Toll » lu-CK 

CPGI2H  Pol 2H 


SUN 719 00 

SUN71910 

SUN71920 

50N7193U 

SON/1940 

SON7.1950 

SON/ I960 

S0N71970 

SUIJ/198U 

SOTT/1990 

SUN7200U 

SON72010 

S0N72U2U 

SON/ 2030 

SUN72n^|) 

SON/2050 

SON  7  2  (1 60 

SON72070 

SON/2080 

S0N72090 

SUN/2111 0 

S0N721J  0 

SUN/2120 

SON7213U 

SON721^0 

SON 7 21 50 

S0N72160 

SON72170 

SON 723 80 

S0N72190 

SON/220  0 ’ 

SDN/2210 

S0N/2220 

SON/2230 

SON/ 224 U 

SON/2250 

SUN72260 

S  0  N  7  2  2  7  0 

SUN722bO 

S0N7229  d 

SON/230U 

SON/2310 

Sun/ 2320 

SON  7  2  33 II 
SON/234U 
SUN/235U 
SON/236U 
SON/ 23?  0 
S0N7236U 
SON7239D 
SON 7 24  0  IJ 
SON/241U 
S0N7242U 
SUN/2430 
SON72440 
S0N/2450 
S0N7247U 


SUBlfOUTlNfc  I  OT?  »  SUN/248 0 

COPPI  FX  VPI;(  ,|iS,l).j,u,PHIW,CA,CFX  SUN/24911 

C0MM0N7C4/VI  IC(  )0,3  5>>  nS(200l)  ),JlG(2H,2),0(3  .M3  0>3),PH1W(5U  ),CK(40  )SON7646  0 
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C0M^0N/C5/HiiH(]  "0  j,NBL(10u>,FU  ,  I FR,  XL,  NS>  NTH,  K»  J*  OK/ (If,  QAB,  QAN  S0N72510 
C08M0N/C6/X*  Y,  OX,  i)Y,  £H*  HK/EKB*  EKR, NP » HP*  NB » nBOX* KUDE , HOUE/NBK,  HBT  SON/2520 
COhuON/CB/RhQ.  S0N7253B 

K-?»HOR  C  NBOX  1  S0H/2540 

N=H  I  NO  (  NB()X  *  1 6  )  SON7255C 

flK  =  bKB  *  SON72560 


DK2=DK**2 
HI sB-l 
llKft  =  DK2/fl.O 
r)M  =  2.0*lJK8 

tiKi2=nK2/i?.n 

08=0.5 
pH=IiK*0  .5 


SON/2!>7Sl 

SON7258U 

S0N72590 

SON7-26DO 

SON72610 

S0N72620 

S0N72630 


DM  =  l>.5»!)H 

n»  =  2.0*l?K 

A  li  H  =*1)11 

ni=l..25*DK? 

H5  =  f.-K2/?4.ll 

no  A  J,=  1>M 
01=11.0 
H4=i  .  U/IiM 

B2=f!5/M-llH 

|i.i-iih#H4+H5 


SON/2640 
SUN7265U 
SDN72660 
S0NA267U 
S0N72680 
SUN/2690 
SON/2700 
S0N7271 0 
SON/2720 
S0N7273U 
S0N7274U 


n4  =  |«K8*fc4 


S0N7-27t>0 


1104=2.0*114 
0=1  .0 
n3=n.u 
C4  =  0.  i) 

C  /  =  li .  0 
CHrlr.  0 

n»  ?  j ■- 1 *  n 

ai=iiH/GM 
C1=(:n*  C0S(61) 

C?  =  -Oi*  SIN(AI) 

CAU  US.|N(Al.,ri>»C6) 
c‘j=ch#(;& 

r.6  =  r(;H*C6 

C9  =  r:j -C3 
cm=n2-C4 
r,  11  =  C  5  -  C  7 

n  i?=n6-r;b 

VHF  =  03*C;9-B4»,:i  0-'<5*C3-Hl*Cll-B2#C12 
V  l,M  =  R.4*(;9  +  B3*.:i  i'-'«5»C4  +  h2*C11-B1*C12 
vph:<  i,j)=opi  y ( v.wF/ v i m.) 

23  C3  =  (,l 
C4=c? 
r.7  =  t:6 

r.o=(;  6 

Hl  =  !  l  -111 
<H.S-s<3-riJ 
H4  =  M-I)4 
|).4  =  fM  +  UJM 
CN  =  (.N*2 .  u 
?  ROMIffUb 
Q,M  =  c  M  ♦  .1 . 0 

nH  =  HMHIUM 

.1'  nnM=iiiiM  *  mi 

no  *.)  j= i>*  n 

DO  4  1=1* Ml 
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SON 72760 
SON72770 
SON/ 27 on 
SON/2790 
SON/2HOO 
S0N728JJ.0 
S0N72H2U 
S0N/2830 
SUN/284 0 
S0N7285U 
S0N7286U 
SuN/2870 
S0N7288U 
S0N7289U 
SON/2601) 
S0N/291U 
SON72920 
SON72930 
S0N/294U 
SON72950 
S0N7296U 
SON72970 
S0N7298U 
S0N7299U 
SON/ 3ll  0  li 
S0N73II1U 
S0N/3D2U 
SON/303H 
SON 7 3040 
SUN  7  3  li  b  II 
SON7306U 
SON  73  U7 II 
SON/3U80 
SOH/3II90 
SON/310U 


■i 


*  K=H- I 

A  VPH.CK+1,  J)=¥PIC(K*1,  J)-VPIC(K,  J) 

5  vpicci. J).=2.a*viii«;(i, j> 

CM=fr.U 

nn=ii.o 
n»H=nK 
DM  12  I=1,H 
C7=0.U 
C8  =  fl.O 
C9  =  fJ.O 
cin=a.(j 
pi  =n.  i) 

P2=n . a 
r,H=i  .0 
H6  =  ii.5*flK12 
Djl  10  J  =  1,N 
ai=cm/cn 
A?  =  iiH/CN 

if  (A  l- 0.2)  7,7,8 
7  Pl  =  P.U-Al«'*2/3.'J 
P2  =  -nK/C6.0*Cri); 

oci  in  9 

0  03=  Sln(A.l)/Al 
R1  =  ?.II*H3 

H?r  (H3-  CUS(A1  )  )/'.2-nH/CN#B^ 

9  «3  =  P0S(A2)/CM 
B4  =  SIN(A2)/CM 
C3  =  l:l*fi3  +  B2  =  R<» 

(M-h2*R4"81>*34 

R5  =  l,  H*CH 

Cl  =  l'5*C4-2.t*#C3 

C?  =  -?‘.n*C4ri  5*C4 

C5-P3-C7 

C6=( 2-C8 

P3=P?-B6*CN 

P4~P3*2 . 0* DM  7*  ( C  i**1  ,  U ) 

VKI:  =  C5-Pl*Cr.r-P3*C,-P4#C9 

V  i  H  =  C6  +  P1*C5  +  P.3*C4-P4»C10 

V  P  J  C  (  I  ,  .J  )=VI'Ii;(  I  ,  I )  tCNPi. X  (  VHE  >  ¥3  H  ) 
PI  =  l’l*  l)H 

P  2= P  2*  bn*  DM 
cn=cn*?.u 
C7--C1 
CH-C  2 
09  =  03 
Cl  0  =  C4 
RA  =  lj£  +  DKl  2 
111  CON  I  I'NUP 
CMr.CM  +  Dh 
n«  =  l;M  *  r,  DM 

12  ni>M  =  |)i)M  +  l)D 

03  =  1  K /  (  2 . 0  * -i . 3  4  L 5  v 2 6 5  ) 

A 1  =  0 . 0 

DO  3  4  d  =  l«N 

CKX=IM*CMPLX(Sr  ' (All,  COS(Al)) 

DO  13. 1=1, H 

13  vp pc ( i ,  i)  =c»:X«  vp n;(  i ,  j ) 

14  A 1  =  A'l +DH 

TURN 

FNP 


S0N7311U 
SGN73120 
SON73130 
SON? 3 140 
SON/3150 
S0N7316U 
S0H73170 
S0N/3J6U 
Sun/ 3i 9 u 
S0N732CU 
S0N7321U 
SON/3^20 
S0N7323U 
S0N7324D 
S0N7325U 
S0N7326U 
S0N/327U 
SUN/3280 
SUN73290 
SUN7330U 
SON/333  u 
SUN7332U 
SUN73330 
SON/334 « 
S0N/3o5y 
S0N/336D 
SON/3370 
SON/ 338 1! 
S0N/339U 
SUN/34C0 
SUN/341 U 
SON73420 
SON73430 
SON73440 
SHN7345U 
S0N7346U 
SUN 73470 
SON/ 34 60 
SUN73490 
SON? 35 00 
S0N7351U 
S0N73  5  2U 
SON 735 30 
SON?”5'4  0 
SON7355U 
SUN/356U 
S0N7357  U 
SON7358U 
SON/ 3590 
SON/36  0  0 
SON73610 
S0N7362U 
SUN73630 
SON 736 40 
SON  73,6  50 
SON73A6U 
SON7367U 
SON 7 3680 
SON73690 
SON73700 
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1 

CCS  I  A 


rnkmAf!  NL‘;Toij#necK 

CSlN 

SUKF0IJT1NE  CSINtXJ  ,U,S) 

Sint-  and  COSINE  InTEGKAI  SUBROIJ  i  I NE 

C  AM.  s  ARE  THF  INTEGRALS  OVER  T  FROM  1  TO<  INFINITY  OF 

Ciistxn/f  anu  siN(xn/r 


s«=i.a 

x=yj 

IF  (X)  1,2,? 

1  S«=-SK 
X  =  -X 

2  X2=X«X 

IF  (X-l.fl )  vs,.),  4 

FOR  AhS(X)  LFSS  Thai.  1  A  SEN IHS  'EXPANSION  IS;OStD 

I  v=(  (  (  */V9H .  i>-il . 6  )  « .  !i5*X2  +  1 . 0,)*X2/I8. 0-1 . 6  ) *X<  3  ;>7079633 
L= I ( X2/45. 0-1 . 0  ) «  x2/24:;  Il  +  i,  6  )*X2/4  .  b-.57721*4o5-AL00 ( X ) 

GO  10  5 

FOP  AHS(X)  GREATER  THAN  1  APPROXIMATIONS  OF  HASTINGS  ARE  USER 

4  P=If(/X2*19.  ‘9  4  1  19  J*X2*47. 411538) *X2 *8. 493336)7  < ( UX2  +  21.361U55) 
1  *‘X2*7U.3  7M9A  J*  <2  +  3u.  1138227  )*X) 


% 

CP  ASP 


0=1  <  (X?*21..‘g  1724  i  *X? *49 . 719775 ) «X2  +  5 • 0895 04  )/(( { (X2  +27. 1779.50)  S.0N/398U 

I  «X2tll9.91  H93»'M  X2+76.707876)*X2)  SON/3990 

CO  =  COS  (X)  SON/4ll06 

Sla'.TN  (X)  SON? 4  010 

l=0»C0-P*Si  S&N/402U 

V  =  P«CQ4  0.*SI  S(JN74Ur30 

S=V*SG  SGN/4040 

W  F  T  0  R  N  SON74U50 

FNP  S0N/4U6U 

K*nTRAN  fJLsT Oil,  fll-CK  SON74070 

MASH  ,  SON74090 

SDRolUITTNE  oASH  S0N/4J  00 

COMPLEX  VP  |  G  <  OS, Oo, 0, PH  I W, CK  S0N74J10 

COMMON/ C3  /Ki  OX  ( 20  a  0  j ,  x£  <  5  ) ,  YE  <  3  ).-,  AR  (  3  ) ,  Xl>*  X2,  X3»  X4 ,  Yl,  Y2,  BETA,  NBS  SON/412U 
C0H(-0N/C3/Ci  (1  0>,2H,2),Nr(10,2),NF(10,2),NrHAX(?),FN(28),DXE(6),TPISON/413U 
COPI.ON/C4/V1  I  C  CIO,  15),  PS  (20  0  0  ) ,  DO  (  28 , 21 ,  Q  ( 1  u ,  1  0 , 3  ) ,  PH  I  H  ( 5U  ) ,  CK  (  4  0  )  SON/6460 
rOMMM/C5/MoO(l<iO  ),NBJL<10U),FQ  ,  IT  R,  XI. ,  NS ,  N  I M,  K ,  J,0R ,  01 , 0A8, 0  AN  SON  74 150 
C'lMh0N/i;6/  X'»  Y,  0  <>  •» Y ,  FM,»  f  K>  EM5,  EKR,  NP,  PP,  Nfl ,'NBOX»  KpDE,  MODE,  NBW,  N8T  S(iN74i6U 


S0N73710 
SON73730 
S0N7374U 
S0N73750 
SON/3760 
S0N7377D 
SDN/378U 
SON73790 
SON73O0O; 
SON7381U 
S0N7382O 
S0N7383U 
SON 73 840 
SON/ 3850 
S0N7386U 
SO.N/38  7  o 
SON/3880 
S0N73890 
S0N739HU 
SON/3910 
S0N7  392ji 
SUN/393U 
SON73940 
S0N/395U 
SON73960 
SON/397U 
S.0N/398U 
SON/3990 
SON  74 II 0  0 
SON74O10 
SON/402D 
SON  74  0,30 
SGN74040. 
SON74U50 
S0N/4U6U 
SON74070 
SON74090 
S0N/4J  00 
S0N74J10 


COMkON/(;«-/RoO 
XP  "  X 

II  (MS.F0.2)  X,P  -  X  -  X3 
NTH  =  N4  (MOi-F/NS). 

7  =  C<IJ(MUUE,1  ,Ns:) 

nz  =  O.ii 

IFUilM.EU.I)  on  111  70 

K  =  1 

'IF(XP.EO.O.O)  G’<  lO  50 

PX  =  XP 

YX  =  Y/XP 

no  40  N  =  2 , 7 

PXY  =  PX 

no  30  h  =  1,N 

K  =  K  +  3 

IFIK.GI.NIM)  i;0  To  2  0 


SUN/4170 
SON7418U 
SON/419’0 
SON  7420  0 
SON/421 U 
SUN/4220 
S0N7423U 
SON/424  U 
SON/4250 
S0N/426U 
S.ON/-127  0 
S0N/428U 
S0N/4290 
SON74300 
S0N74.51  0 
S0N7432U 
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?P=i*XY*CO(Hi»or^of<S) 

S0N/433U 

? =7+,2P 

S0N7434U 

-- 

OZ=l)Z  +  FJHX  J«ZP 

SUN74350 

Jn 

FXY=PXY*rx 

SON74360 

4 11 

PXiKX#XP 

SUN7-4370 

2d 

n/'iiz/xp 

S0N74380 

no  7  0  7  &  , 

S0N7439U 

6n 

PXiJ ,0 

SUN74400 

no  60  N=2,7 

SGN/441 U 

K  =  K  +  N 

SGN/4420 

1 

X. 

II 

& 

SUN/ 441(1 

irok.l.£.NT».  )  nZ  =  liZ*PX«CO(HODfc,KK,NS) 

SUN/4440 

P  X  =  P  X  *  Y 

S0N7145U 

Jf(K.GI.NTM)  nil  TO  70 

S0N/44O0 

on 

7  =7  *  P  X  *  C  0  ( H  < «  H  I- ,  <  ,  •'<  S  ) 

SON/4470 

7  it 

GO  JO  (HO, 90), H$ 

S0N/44.H0 

on 

nS(hB).  =  fcMH  XP'JZiEK*Z7 

SUN  7  44  9.0 

RFTuPN 

Sun 7 46 u  u 

9  n 

OS  (PR)  =  GMPI  X(-iZ,EMZ) 

SUN/451U 

RFil URN 

SUN/4620 

r-Ni 

Sun/ 4i>3u 

.1 

f  hi.  1  KAN  M.S-TuU«'PcCK 

SDNA464U 

CPhin 

PHI  B 

SUN74660 

SUPkuiJT  I  Nf  i-HfS 

Son/ 4670 

GUPPIPX  VPfr;  i)S,  l)u,b,  PH  l  W>  CKyCZI  RO,  PH  j  ,  PHITE,  UPHI  ,SPHI 

SON/468U 

r;iiHpfiN/i;4/V»r !.:('!« >.1S  y,  DS(20UO'),iJ<M28,2),G(;l-0,10>3j,PH:iW(5U  ),CK(40  )SUN/646U 

COHMOM/OS/H'-HUr'O  )>  W.<l.(100),F(l  ,  1FR,XL,NS,NFH,»K,J,QR,  oi>  oab>  oan 

SON74600 

COMPON/r;6/X,  Y/pA,.iY,‘:-K,bK>EKB<EKR>Nf;,PP,Nb,H8LiX,KODfc,HODE»NBW,NBi 

SGN74611i 

C0MM0=\/f;7/c/PRn,  pH,,  t,ii’UE>  DPHP,  SPHI 

SON7462U 

cohmon/gh/rho 

SGN74&3U 

MJi-MIMn  (NP>  15) 

S  U  N  /  4  6  4  U 

no>  pn  i  =2^=K'<j 

SUN  74  65  il 

NO-  i'lP- I  +1= 

SUN74660 

JP=MI)0(WU) 

SUN/ 4670 

fv  J-NRL  (  NU ;  +)' 

SUN/ 4^80 

no  po  .ui,ji 

SliN7  46911 

K-1  +  1  AHS(llP-J) 

SUN/ 47 00 

fjPH  1  SVP  I  C  (  K»  F) 

SON747I.U 

IP  (.1.10.1)  GO  ro  ID 

SUN7172U 

K  =  MI>4  ,j-j 

SUN/4730 

uph  i*i)Piii+vPii’(  <-*  i ) 

SUN7 4/^0 

J II 

SPHI  =SPi*j+DKHJ«  'S(N.i) 

S  0  N  7  4  7  5  0 

2n 

N J-uJ+1 

SUN?  476 1) 

P  b  1  1  U'ei 

SON,/ 4  7,7  it 

P  NU 

SON/478  ll 

(OUTRAN  NL  jT'OIJ,  DF  CK 

sun/  4:79U' 

CJ'STi'f-C  MSI  HRC 

SuN/481 U 

FUNGI  ION  MS  l-MFC  (Mi  N,  |.,  A.,  B  ) 

SUN/ 4 820 

r.oHPi  ex  a,  bug 

SUN  /  4,83  !k 

[J  I  M 1  N  S  1  l)N  AJ  M,  1  ),  ‘Hr.,  i  ) 

SON/484U 

no  3  0-  I  -  U,N 

SUN/48biOi 

c  =  o,.o 

SUN/4860 

no  m  j  =  i,n 

SON/487 u 

lu 

C=AMAX1  (C,  AKS(  HI: At  (  A  (>J  ,  J  )  )  ) ,  AB S  <  A  J  M A G.(  A  (fI ,  J  )  )  )  ) 

SON/48811 

FT  (C.hO.O.O  )  GO  TO  7  011,0 

SON /4 890 

no  ? 0  J  =  1  , N 

SUN/ 49  0  0 

2(1 

a  ( [j>  J )  =  a  ( i ,  ,i )  /n 

SON74  9iO- 

II 0  ■)  0  J  =1.1 

S0N7492U 

v5  0 

8(  i,  j)  =  Ji(  i ,  j ) / r: 

S0N/493U 

FF(N.J=G.l)  no  TO  20b 

SON/494  0 
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N«  -  N  -  1 

S0K/4950 

no  20.11  .!  =  1  ,  fJM 

S0N/4960 

C  =  0.0 

5GN7497  U 

K  =  0 

S0N7498U 

no  4  0  I  =  J,N. 

;S0Vi7499U 

n  =  A«S»R£Al  (A(  |  . J)  )  )-*  AUSUIMAGCAi I, J}) ) 

SON75000 

i  f  u  .  t;E.D)  i=n  i  i  io 

SON/ 5 (i*3  iJ 

K  =  i 

SON75020 

C  =  ji 

SON/5030 

4  II 

COM  l  N 1 1 1* 

SON  7  5  MU 

If  (f  .tO.U.OF  .-C.U  .1  .F-7,)  GO  TO  1000 

S  U  N  7  5  U.5  U 

if  (r.fo. j>  <;n  t  i  ;o 

SON/50V6U 

III!  .  0  JJ  =  .!„f. 

SON  75  07.0 

r.=A(..i,  jj) 

S0N75ii8u 

A  ( .1.  .1.11  =  ,A(K  ,J  J) 

SON  7511 90 

5*1 

AIM  .1.1  )  =  G  - 

-.SON751  Oil 

n f»  M«  j.jj  =  i,i. 

SON:/ 5 1.1  u 

g=h i j . j j ) 

SUN  7*3.20 

H(J,.JJ)  =  H'K.JJ) 

SON/5livO 

R  (  K  ,  J  J  )  =  G 

SON  7  51  4  il 

/••l 

G=1  .11/ At  J,  .11 

SON/5150 

IP  =  .1  +  3 

S0N/5li6  0 

nn  Ml  JJ  =  »P,  N 

SUN  7  5 17,0 

Ml 

A  IJ)  =  A(  J,  1  D*R 

S.0N/518U 

9  <1 

n li  i nn  jj  =  i.i 

SUN  751 911 

l  o  n 

H(  I .  J  JJ  =  W.(  J  ,  I.I )  *14 

S0N752P  II 

nil  2(10  I  =  3.'! 

SON  75.2  111 

ii  ( i  .to.  J)  on  Ti<  2(in 

IS  ON  75.22H 

G  =  MI,J) 

S  ON /.523D 

nn  ii>o  jj  =  jp.  i 

SON/ 5.24  0 

11  n 

4(1  ,  JJ)=A(  I  ,  JJ  )  —  f» «  A  (J,  JJ  ) 

SON75250 

f)!l  1211  JJ  =  3,1 

S0N752&U 

1  2n 

B(  1  ,  J  J  ) -  B  (  I  ,  J  J  )-R«  B  (  J  ,.J  J  ) 

.SON  7  5  27  0 

2  (i  ii 

COr-UMif 

S0N7528  U 

,  2(10 

R=A(M.M 

:S  UN  7  5  29,0 

if  ,(  A  MS  (  KF  A 1  Ci!)  >  ♦  ABS't  A.!  HAGXli ) :),.  1 T  .  1  «.E-7  )  GO  TU  100:0 

.SUN /53  On 

no  /111  J  =  I  ,4. 

SUN  7*  3 1,0 

2  3  n 

H(M,.j)=H(N,.n/n 

SON./*  32.li 

I  r  ( ■•.  FO.l  >  1-0  10  23 II 

sS-.0N7fS.A3iU 

rj  0  .’2  il  1  =  3.  f|p 

S  0  N  75, 3  4.0 

no  220  JJ  =  3,1 

SO  N7*. 3  5,0 

22  n 

H((  1  ,  J  J  )  =  B  (  i  ,  J  1 )  -  A  (  I  ,  N)*h(N,  JJ) 

ISON./ 5  3  6,0 

2  01' 

PSrriFC  =  l 

SO  N  715.3  7  0 

RETURN 

SON/** 80 

3  II  li  ii 

PSI  MEC  =  2' 

SON  753911 

PP1 URN 

SoN./*4i0iU 

MI» 

SON/*  41.11 

i 

fOl’TK/Ui  NU  T'ilMltCr. 

IS  ON  7*  4,2.0 

n-iic  1  .j 

no  i  j 

S0N7544.il 

SOWMIinilfh  1  Q  I  J 

S,0N,/*4‘5iU 

COM  PI  ex  VP  i  ,|.S,  II-..  U.,  PM  I  w,  CK  ,  CZfcRU*  P.H1  i- PHI  TE.UPhL.S  phi 

S:0  N,7*  4  6:0 

TOP'  Oo/Cl/Ki  oy(  20u0  ),  XF<5  ),  YE  (3.)  ;  AR(*.),.Xl,X2>  X3,  X4.»iVl,  .Y,2,!BE1.A;,  NI»:S 

S.0N/*47‘U 

C  OHM  (IN  / 03/  C;  (  1  n  ,  2'.,  2  ) ,  N  1  ( 10 , 2.) ,  NF  W  (1 , 2)  .NT  MAX  (.2.)  ,!f'lN.(;28:)  >|DXE\<'6  )  ,|I!P,IS0N/S4:8U 

COM!  0:1/ 04  Ay1  |i.  (  nil ,  3  6  )  ,.i)S(20U0  ) ,  1>Q(28,;2.) ,  Q  (l.u ,  10 ,3)  »iP,Hl.W((*!U.})  ..tCK14iO)SON/,6  4,6iO 

COMM)N/t*/Mi»nQ  10  ).,  l*HL  ( 3  Oil )  ,FU  ,1PR,  XL  ,  ,N  S .,  :N  i|  M » ,K  >  J.»  ,Q  R ,  U  ii .0  A’B  >«0  AN 

SON,/*  6:0,11 

0DMP0i'l/C6/ X  ,  Y,Ji)X,  ii  Y ,  t  rt ,  fcK ,  EK8 .  FKR ,  NP  ,,MP ,  NB  #  NB  0  X..:K  ODt  »:M0  U,E.»!NB.W.,;NbiT 

SON,/*  6 1.0 

C!' )*•  I- C N / C ?'/ C / F M) ,  P t ■  I ,  P H 1  I  E,  I) PH  1  ,  SPHl 

:S,on  /*  62ii 

COPi-iON/OB/R.-n 

S‘ON.7*  6.3:0 

OPH  1  =  PH  1  *•  Xf  (  SO  IF  ) 

SON. 7**4. U 

IF(K-P.FO.I)  DPMI  =  II  .  V  *  llPil  I 

SON  7  5  65‘11 

s 


'X 
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XP  =  X 

I RfN5.pl*. 2)  XP  =  X  -  X3 
PX  =  XP 
K  =  1 

NT*  =  NTHAX(NS) 

nO(l*NS)  =  non  .  N  ;S)  +  DPH  I  *CHPL>:<  Q  ;  0  ,~EK  ) 
JMNTH.BQ.14  GO  ro  50 
BiJ.O.i' )  0*i  10  3n 
Vx  ■=  y/xP 
IliJ  20  N  -  2,7 
i 7  =  PX 

no  in  n  -  i, n 

k  =  k  -+  3. 

;MMi.;OT..NhM;  i;Q  Tu  b  ii 

nUTK/NS)  =  l.'HKVNS)  +  •.PHH.1*Z*CHPLX  ( F  N  (  K  3  /XP  ,~BK  ) 

1(1  ?  =  2  *  y.x 

2il  PX  -  PX*XP 
oil  10  5(1 
3n  px  =  i.,,n 

-DU  4,0  m  =  2, 7 

;K<  •=  K  +  N 

;kk  ==  k  -  1 

TK(MK,.'I  B...IJT’  =)  n.»(\K,JNS)  -  HU ( KK, NS  )  +  DP  HI  »PX 
PX  =  P X  *  Y 

TP  ( t-:i  ,41..'NTM  )  00  Vi  t  9(3 

4(1  6  0  <  r.  > '  5  )  -=  "»)PK  ,  NS  3  •+  TlpHl*OMPLX(U.O,-.EK*PX')  . 
t>'i  lFXTOT(-,.NE,.;v)  Bo  Id  9IJ 
ii ph  j  *  hPur/nxF  (K-iiip  ) 

K  -  a, 

yx  -  Y./XP 
PX  -  "...O/IJX 

nn  /n  n  =  1,7 

7  =  (P  X 

l|MI  Ml  M  -  1,  N 

nn(ip,NS)  ■=  mitk.nv)  -  n PHI *./ 
k  -  k  +  l 

;|;f'fi',.oT,.NilM)  mi  if.  Mi- 
<Mi  7  =  7  *‘iY  X 
7  II-  P  X.  -  P  X  o  X  P 
hu  AKG  •=  RK«,CX*'  -  Xtl  * 

ill  PHlI  =  -jD:P>Hil  »  C  M  I’d  xi(;OS;<vARO  ).,  -.S  I:N  (  ARC  )  3 
XI  =  0..I0 
iO  AI  !l.  tOE.OF 

,90  iliFdihiO n,F..tNB..‘3'.)  iRfc1T*JR»v 
X  I.  ■=  'X-<1 

i!!F'(«iPi.>PQ..l  )  'PmiIsPF  =  <)..'•¥ *PH!Ii1:E 
ill  PHI  -  iPHiliK 
CAl  I.  OF  OH 
!REir">RN 

lE  H  P 

1  lF<0 RiTiRAiN  iNL'‘'-!Pnil.,:UH(M< 


S0N75560 

.S0N7557  0 

SQN755813 

S0N/559H 

S0N7560U 

S0N7561U 

S0N7562U 

.S0N75630 

S0N75M  0 

S0N75O50 

SON/56.60 

S0N75o70 

SON75M3U 

SUN,/5ft9  U 

SUN  75.70  II 

,S(jN75/1iJ  . 

S0N75720 

S  ON  7  5731} 

SHN75/4  il 

.S0N7575U 

;S0N75/6U 

.S0H7577U 

S0N7578U 

S 131797  91k. 

,iON75hflU' 

S.UN/581U 

-SON  75,820 

SUN7583U 

5 ON 75 8 4  0 

S  0  N75  «  5 1) 

SUN/J»860 

5007987 u 

SON  7/58  8  U 

SUN/589U 

S  0  N:7  9  9.0  II 

;SON /591U 

SON/5920 

IS  UN  /  5  9,3  U 

SON/59411 

S0N/595U 

:S  UN  7*59  6  II 

S0H/59.7.U 

.SON759HO 

.SON  759  9  U 

5UN76  0  Oil) 

:SON7.6|iaO 

:S  ON  76020 

(SON  A6.1J.30 

5  0  N  A6  0  {f  l) 

:S  UN/ 6  6,5,0 
lS  ON  7611 611 


if'illEiGr-  (Cl  FiG  h  :S0N/»6,U;8l) 

5‘UHROimiiiE  I'.I-OB  :SON7'6ii9U 

tC.oHhl  .EX  'V:PiJC„.|iS.,i»..i,.n  ,ii’Hil!W,;CK*i!C.2tfRtO>(P.HlI.»iP,Hil  (IiL»!l).PH>l,*;SP,H!l  :S.(iN/6HJ  u 

'CiOPr,no/iC4,/..Vl’iIa'.t(  oil •»  .1*5 !)  >» if).S((,20llnp ') ,» iD.QTP.O  >2  >) •» -Qli  *l  >  J<0„.3])  .iP.HIlWTSi'1)  >iCKi(-4  0  )50N/  64  611 
iC'O H fi.O N /(05,/'M i  'ltd  iiO')  ,.NHUliOi01)  ,;FtQ  liHR„.‘Xi ,iN.S„!NlliM.,iK.*.J.»iQR,,iUiI  ,.QAfci  ,.QANi  'SuN'A6120 
.CO Mt'.O si /tC6  7.X.,  T  /n .<,Y.,iB Mv.r;,K  ;.BM3  ,lENiRo'N.H  ,iP;P  „N (),,  v, H U.X.„Kn).D  t  ,  M 0 fl E  *  N,B  w  „N til  :S ON  76  I  3  0 
,C.OM,  0?l 7(0.7  Aciil.-t/'O  AP  III ,2, Ml  TE  aDPHII  AS’P.hU  SON 76 1-4.0 

(COMiin.'l/iOS./iRt  .0’  ,S ON 761,5 U 

!p:X-a..'Or/!l''X  S  0  N  7616  U 
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j* 


A 

o‘ 


c 

5 


<5  * 


V 

*•*  Jfi 


il 


■a 

i 


!F{  XL.FU.O.H )  fi'J  10  30 
Y  X  =  Y  /  X  l_ 

K  =  1 

no  20  n=1  ,7 
PY  =  fX 

nil  JO  M  =  l,N 

1)0(1',  2)  =  ])u(K,2)  -  OPH !  *PY 
lFtK.EO.NrM)  KFTUkN 
K  =  Ki  3 

111  PY:=PY»YX 
2ii  PX=I’X*XL 
RETURN 
<>  (I  K  =  ft 

1)0  4  0  M  =  1.7 
<  ~  K  N 

IF  <K  .I.I.NlH)  REHIrfN 


SON76I70 

S0N7618U 

S0N7619U 

SON7620U 

S0N/621U 

SON/6220 

SON76230 

SdN76240= 

S0&7625H 

SUN7626U 

S0N7627U 

SON/ 628 1) 

S0N7629U 

SON7630U 

SQN/6.liu 

S0N76,i?U 


4  il 


>} 

OHfl'T 


hlHK,?)  =  If  (K,7)  -  OPHl  *px 
p  X  =  PX  «-Y 

rftiph 

FNn 

l  ni'  i  ‘i a fi  ns  M.iij.  n>  On 
-P  iiii=I 

SORPOllI  INI-  I  mil  t  I  id) 

COMPLEX  VPIi  ,  liS,  I)<J,  ij»  PH  I  W  ,CK , CZERO-,  PH  I  »  PH  I TE, OPHFi SPH 1 
I)  I  R|;NS  ION  $!■{*».  >)s  Sl'RF ( 2» .i ) » COd.O.7  )»C(-50 ) 
niMINSloN  CAP'tS(?->»*jO  ) 


S0N7633U 
S0N/634U 
SUN7$35?I1 
S0N/636M 
SON/ 6370 
SON 76390 
SON/ 6400 
SON  76  4  111 
SON/6420 
-**i'#**** 


CIIMPUN/OJ /K.^nX(MllilO) .  XE(5)i<Y£(3>,  AR(3),X1,XH>X3,X4,  Y=l,  Y2,BE1  A,NfiS  SON/6430 
H 0 M M 0 N / C 2V A s >  iH 4C-i'»  f  hACH( 6  ) ,  NFRbQ,  FREC  (,10  ) »  NHODE,  NSURF, LVP I C >  LSSVPSUN7644U 
n')Mf*0N/i;3/ci‘(  \0,2<S,2\,  NKiO,2),NF(l>,2 > » N  TKAX  ( 2  )  >  FN  ( 28 ) DXE  ( 6 T,  TP  1  SGN/64&U 
C0HM0N/C4/VI  I0(  HO,  15).,  OS  (20  00  b,  l)Q(28>2>,Q(lli,lO,3),PHIW(5U),Ck(4U)SON7646U 
C0MHUi|/C5/MiiR{JO0)/NOL(  JOU),FO  .  I F R,  XL , NS , N TH» K, J, OR, U I , OAB, QAN  S0N7647H 
nOMpOf,/ 06  /  X  ,  Y ,  Dx ,  >•  Y ,  tM,(E.K » EKB,  EKR, NP >  f?P*  N8 ,  MB  OX  ,  KUDfc,  MODE,  NOR,  NOT  SUN/6480 


>1 

2 

o 

4' 

r. 


r.0Mhi)i|/C7/C/.PKn,P*l! ,  Pill  TE,DEHI  ,,SPHI 
CiihmOn/cb/R'-o 

DATA  (  Swtl,  I  )  >  !~1 .6  )/26||MAP  OF  SONIC  BOX  OVERLAY 

26llONt  RING,  TAIL  AND  R4KE 
26H  (S')  -  RING 

26 H  ($>  -  TAIL 

26'H  (,)  -  RAKF 

26  Fi 

data  FSDRF  (  i  ,  i  ) ,  1  -1 , 3  }/HHW  I NG  ,8HTAIL 
DATA  cnii/i<H,:,4  H  F,  I  HS,  lHli ,  1H3>,  1H, ,  1H .  / 
fit)  50,  (10  ,211  ,  <0  ,40  ,dO  ,60  ,  70  >80  ),INl) 

10  R W  ‘  I  b  ('6 » 1 1’ )  EM,  AS,  PUD.,  x£  ( 1 ) ,  XE  (  4  ) ,  XI ,  X4,  Y1,T2,  T2,  Y2,X2,X4, 

I  A R  ( J  ) ,  Al<(  2  ) ,  NOR »  N;i  I ,  NDS,  NIIS 

11  FORMA!  (1HJ//////3 *X,43HHAC/NAA  MISSILE  IRANSONJC  AIRLOADS  PROGRAM  SUN7662U 
1  ///37X>  3mHI  L  I  G»<T  Cn.il)  I  I  I  dNS  A N I»  GECMfc  TRY/1  H0.//15X,  1.3HMACH  NUHBERS0N7663'.f' 
7  =,l  R.S.4X,  l6-iSiJE!  R  01-  SOUND  =F10.3,4H  L/T ,  4  X ,  4rlRH0  = ,  E14 . 8//  1H0/  SON7664U 


/ 

llriWING  +  TAIL.  / 


S0N7649U 

SOI97  6&OU 

S0N./651Q 

S0N/6&2U 

SO,N76d3H 

S0N76d4U- 

S0N/655D. 

SON/6560 

SON/657U 

S0N7658U 

SON/,6590 

S  0  N  7  6  6  0 II 

SON76610 


X54$,4HW1NG,  IHX, 

.5  4 M  |  A  I  i. / /  / ? 7  Xj,  1  6Hi  .1 
2F  72.3//  2/X ,  1  Mil  .! 
2F72.3//  2/"  X  ,1  *>H  •  i  P  CHORD 


S  I  A  I  ION  (L), 21-22. 3//22X,  16HR00T  CHORD 
SPAN  (  L  ) ,  2 F  2 2 . 3 / / ?,2 X ,  1 6 H t ,  t .  SPAN 


S0N/665U 
(  L  ) ,  S0N7666IJ 
(  L  ) ,  SUN76<>7 1) 


(L),?F,22.3//27X,16HT0!AL  AREA  (  L*L  ) ,  SuN/6t>8U 


2.F72.3//  27X,1-..Hi:miRllWISE  BOXES.  ,  1 19,  1 22//  S0N76I.91I 

? 7 X , 1 6 H S P A N W I S E  BOXES  ,.  119,  122)  SDN76700 

WRilF(6,12).  ROX,IU,ilY  SON/6/10 

1?  fOPMAT  (JHU/,1 1  X ,  2oH  I  0  f  Al  CHORDWiSfc  BOXFS  =,!3,  liX,llHB0X  CHORD  =,  SUN/672U 
.1  1>P  IF12 .  i> ,  21  I,  5  < ,  1  'I  hHIjX  SPAN  =  ,4P1F12 .5 ,  ?H  U  )  SUN76730 

RR I  I  E  (  6 , 9 1 )  S 0 N 7 6 74 1) 

Nil  =  1  SUN76750 

no  17  NP  =  1  .  IR  IX  SUN76760 
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II 


MB  =  HG8  ( NP } 

HO  13  HP  =  1  , MB 
K  =  KHOX(fcR) 

C  ( HP  5  =  C'OD  < K  ) 

13  NB  =  MB  +  1 

IF  {Np.GT.6)  iHl  10  15 

fcftME(6>i4)(SH(  I,  HP),  I=1,5),TC(HP),,HP=TL,MB) 

14  F  OPhAT(10X,5A6,<>OAl  > 

GO  10  17 

15  WRITE (6, 16)  (  i;  C  -(PJ.P.Pd.HB) 

16  FOPMAF(4  0X,5nA1 ) 

17  CONlJNllf- 
GO  10  logo 

20  klMIF(  6,.5l)FRt:0(  If R ) , NROXVEKR  ,  EH 
WR  I  I F(6»21 )  t-H.fcKB 


SON 7 6 770 
SON76780 
SON/6790 
SON76800 
S0N7681U 
S0N76820 
SON76H30 
S0N7684U 
S0N76850 
S0N7686U 
SON76870 
SUN76H8G 
S0N7689U 
SON/6900 
SON/691 U 


21  F.lRf'AKlfl  »  28X » 48 -iPi  4NAR  VELOCIlY  POTENTIAL  INFLUENCE  CUEFFIC1ENTSS0N/692U 

1  /1n0,30X,5i*H»iACH  Nii.  =  ,  F8..5ilQX‘,  6HKBAR  =  ,  F9.5/1HU ,:L3X„lHf  5X,  1HJSON76930 

2  inx,5HNUBAI'<HX/5l:HU8AR,llX,i4HkEAL  VF IC  (  J>  I )  ,  8X,  14H I  HAG  VPIC(  J,  I  )SON76940 

3/1 H  )  SUN/6950 

J»Oy  =  2*HOl'(plBiiX)-2  S0N7696Q 

I  BOX  =  H  i  NO  { NiiOX,  1  5  )'•  SGN/6970 

PO  22  I  =  l.HOx  SON76980 

RO  22  J  =  1,.'J.«0*  S0N7699li 

HAPnil  —  I  r-  1  SON/7O0O 

BAPHll  =  o-  -  t  SON77U1U 

22  W  R I  T  F  { 6 , 2  3  )  1 , J.BaRNU. BARMU. VPI C ( J, J )  SUN/7H2(J 

21  FnRHAr(9Xf2I6,4x,.;Fl:(l.l,2X,lP2E22.5)  SUN/7ii30 

GO  10  1000  -  SUN/ 7.04  0 

Jil  Uin  i[;(6#&l)  FKF  J(  IFld.NbriX.FXR  .EH  S0N/7U5U 

kRMF(6,-3l)i*OHF  SON7706U 

3.1  F Ul/.MA  1  ( 1H  .21  X.o9mUPPER  VELOCITY  POIEMULS  AND  SURFACE  UPWASHES  FSOn/7070 

(COR  M0HF  NO..  13/  lilO,  9X,  1HN,  6X  .1 IIH.  5X  »  2HNB.  7X,  10  HR  PH  I  ( N,  M 1 .  7X.10H  1  SON/ 7 080 
2  PH  I  ( N.  M ) » 10X .  9-|R  k(N.H),  «X,9hj  W(N,M)/1H  )  SON/7II90 

GO  10  1000  SUN/710  U 

40  SPH 1  =  DSCNH1/E1  SON/7lio 

WRl  II  (.6,11  )?.P.H-'.NB,SPHI  .CKlHPl  SON//120 

41  PORHAJ  ( 4X,  3  I  7 , 1  P2irl  7 . 5.  3X,  1P2L1/.5)  .  S0N/713G 

GO  1  0  10  0  0  SON  77,1 4  0 

6ii  I F I  OS.FQ.3 )  GO  10  100  SON/7150 

5.1  hi?  I  i F( 6, 51 ) I  RI-f/('I  FR } , NHuX , EKR  »  EH  SON/7J60 

51  F  OPMAT'd-Hi,  33x,  -43mm  I  SSII  E  TRANSONIC  AIRLOADS  PRUGRAH  (CONT-D  )//lH0SUN/7  j.7u 

1  8X.27H0SCIM  AT*'RY  FREQUENCY  ( CPS  ) , F 12 . 5 , 1 4 X, 1 2 , 25H  BOXES  IN  CHORUSON7/18U 

2  IIIHFCTION  /lM0,8X,3flHREDUCFD  FREQUENtY  (SEMI  CHORD  ) ,  F9 . 5, 14X ,  SON77190 


3  2.3PFREE  SThFAH  MACH  NUMBER ,  F9. 3  ,./lH  ) 

HRT  l\F(6, lilt)  SOkFd  ) 

101  F  UP'-iA  r  (  28X ,  i5»l  I  iPoT  MODE  SHAPE  POLY  EOF  I  A'L  COEFFICIENTS  FOR 


SON 7 72 00 
SON77210 
SON/7220 


J  22XV62HREFFRI  NCE.i  L<1  THE  SURFACE  LEADING  ElJGb-CENTLRL  I NF  T  NT  ERSECS  ON/723  0; 
21  ION  //2X,4iiMUl!fe,20X,  7HCUEFFS.)  S0N772.4U 

HO  102  l=l,>Mfi[ifc  SON77250 

N  1  M  =  NT ( I , NS )  SON/726'0 

02  W  R  PI  F  (  6 , 1 0  3  )  .1  ,  (CO(  1,  J,  NS  ) ,  ;J  =  1. ,  NTH  )  S0N/7271I 

Id  FORMAT (1H0, i4,4X, lP/Pl3i4/(9X,lP7E13. 4) )  SUN/7280 

Oil  HR  I  I  E  (  6 , 51)  P  RE  (J  f  IKR  ) ,  NBOX,  EKR  ,  EH  SUN/7290 

KIM  1 F  C6 » 61 )  (SURF  ( I »  NS  ) » I  =1, 2  )  SON'/ 730  0 

61  FORMA  F  ( 1H  ,  35X,  23HGFNERAL  1  ZED  FORCES  F  OR  »  2  A6/1H0  i.6X ,  4HDEFL  >  3X »  SON/,7310 

1  4 H I  GAD.l  OX,  9 H REAL  PART ,  ll)X>.9H  I  hAG  PART ,  10 X ,  9IIAUS  VALUE.1UX,  S0N/732U 

2  1-1HPHASE  Ai.GLF  //)  S0N/733U 

0«  10  1000  ,  S0N/734U 

7(1  WJ?!lF(6,7l  KJ.K,  JR, 01  .QAH.UAN  SON/7350 

KKK  =  ?*NMQ|)E  **t*«»«» 
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NNN=?*K  ******** 

NNHN=N!1M-1  ******** 


CAPOSt  J,NNNM=QH 
CAPOS (  J,  NNN )  =il  1 
IF  INSURF  .1-0.  i)  BO  TO  667 
IF  CNS.NE.3)  lin  TO  3  00  0 
66 7  IF  (J.NE. RHODE  >  GO  in  inOU 
If  (K.NE. RHODE)  GIJ  10  1U0U 
KKK=2*NH0r)F 

PUNCH  666 ,  i  ( CAR'iS (  I  I ,  J J ) »  J  J=1 » KKK ) •  I  I=l»NHODE) 
666  FORMA  T ( 6E12 . b ) 

/I  FOPHATUHO#  19.  |/,/X,l<P3El9.b,0PPi6.3,4H  DEG) 
mi  io  loon 
bO  RR I  IF (6.913 


******** 


*#•*••** 

******** 

*******v* 

******** 

******** 

******** 

S0N773&U 

S0N773711 

S0N7738U 


kR1 1 F ( 6> 81 )  FWFj(IFR)  *  S0N7739U 

hi  FUPhAf(lfUi//////?4X,b6HUNA8LE  TO  OBTAIN  COMPLEX  SIMULTANEOUS  E0UATS0H/74UU 
HOPS  S0LUTIn?I//.m,42llC0HPUTAT!0t(S  I  OR  THIS  FREQUENCY  TERMINATED, /S0N/74iu 
?/49x,  6HFKE-.  =  #  I  9 . 3  )  SON77420 


9t  FORMA  I  (lHl,3ftx»  »3  iMISSI.t  E  TRANSONIC  AIRLOADS  PROGRAM  (CUNT-D)//)  Sun7?430 

3:0 (in  RL  T URN  ‘  SON/744'I 


FUii 


S0N7745U 


5,0  SUPERSONIC  UNSTEADY  AERODYNAMICS 
5,1  Theoretical  Derivation 

For  Mach  numbers  greater  than  1.6,  the  linearized  flow 
equation  becomes  a  hyperbolic  differential-  equation: 

(M2-!)  4>-  -$  ,  -$  i-M2  2ilc  $  -  k2  £  .  . 

■  3cx  yy  Yzz  |_  Yx  YJ.  (5.1.1) 


The  equation  Only  has  solutions  within  characteristic  regions  called 
Mach  cones.  Linearized  supersonic  flow  theory  has  led  to  closed-form 
solutions  for  many  types  of  lifting  surfaces  in  steady  flow.  These 
solutions  are  easily  derived  because  the  influence  of  a  small  perturbation 
is  confined  to  its  downstream  or  aft  Mach  cone.  Conversely,  the  only 
disturbances  that  can  influence  a  particular  point  are  confined  to  its 
upstream  or  fore  Mach  cone.. 

Equation  (5.1.1)  is  satisfied  by  a  pulsating  source.  The 
source,  placed  at  (C,h,£)>  eumiiates  spherical  disturbances  and  has  a 
velocity  potential  induced  at  (x,  y,  z)  given  by 


$g  =  A  (5,  h,  ?)  G  ,(x-C„  y-n,  z-£)  -  (5.1.2) 

» 

where 

G  (x-5,  y~n,  z-s)  =  -  exp  [-  ifc  (x  -  5)1  cos (5.1. 
R  -V(x  -  5)2-(M2-l)[(y  -  n)2  +  (z  -  cl2]  (5.1.4) 

k  =  k  M2  /  (M2-l)  (5. 1.-5) 

arid  A  (5,  n,  5)  =  source  strength 


This  type  of  disturbance  has  no  influence  outside  the  down¬ 
stream  Mach  cone  and  is  discontinuous  at  the  point  (5,  n,  £}.  To  provide 
the  necessary  antisymmetry  of  disturbances  associated  with  lifting 
surfaces,  we  place  a  pair  of  sources  on  either  side  of  the  z  =  0  plane  and 
require  the  lower  source  strength  to  be  equal  in  magnitude  but  opposite 
in  sense  to  the  upper  source  strength. 
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A:>pl\  lag  this  source  superposition  technique  to  the  wing. 

Wake  auc  control  surface,  the  entire  z  =  0  plane  is  covered  above 
and  below  with  source  sheets.  The  strength  of  the  ^pposing  sheets  are 
equal  in  magnitude  but  opposite  in  sense;  -A(£,  n»  0  )  *  -A(£,»  n,  0  ). 
The.strength  distribution  has  been  shown  by  Garrick  and  Rubinow  (Reference 
8)  to  be  equal  everywhere  to  the  lofcal  downwash.  Using  this  condition, 
we  have 

A  (5,  n,  o+)  -  w  (5,  n,  o+)  (3.1,6) 

+  4  > 

The  velocity  potential  at  (x,  y,  0  )  can  be  written  as 

$  -fj w(5,  n)  G  (x  -  5,  y  -  n)  d^dn  (5.1.7) 


The  range  of  integration  extends  over  the  region  contained 
in  the  fore  Mach  cone  of  the  point . 


In  the  wake  region  between  the  wing  and  control  surface,  the 
downwash  must  be  determined  to  satisfy  the  continuous  pressure  condition. 
This  is  the  same  relationship  given  by  Equation  (4,1.7)  in  the  transonic 
unsteady  aerodynamics  discussion.  Substituting  Equation  (5,1.7)  iii  the 
wake  condition  gives  ‘ 


fj" 


(5,  n)  G  (x  -  5,  y  -  n)  d£dn  =  4* 


‘T  • 

exp-ikfx~x  .  1 

t.e.  L  1  Wing  t.e.'J 


(5.1.8) 


This  relationship  which  requires  knowledge  of  the  upstream  down- 
wash  within  the  fore  Mach  cone  to  solve  for  the  local  downwash. 

The  preceding  discussion  applies  only  to  lifting  surfaces 
with  supersonic  leading  edges.  When  the  leading  edges  are  supersonic, 
the  fluid  flow  over  the  top  and  bottom  of  the  lifting  surfaces  are 
completely  isolated  from  one  another.  To  treat  the  problem  of  subsonic 
leading  edges,  the  concept  of  a  diaphragm  of  unknown  downwash  is  added  to 
the  planform.  The  diaphragm  region  is  between  the  Mach  envelope  and 
subsonic  edge.  The  diaphragm  downwash  is  found  by  applying  the  condition 
of  zero  pressure: 

0#  (£,  n)  G  (x  -  5,  y  -  n)  d?dn  (5.1;9) 

The  downwash  of  the  diaphragm  creates  an  additional  con¬ 
tribution  to  the  pressure  profile  over  the  planform. 
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the  velocity- potential  distribution  is-  found  by  the  so 
called  Mach  box  method.  The  wing,  control  surface,  wake  and  diaphragm 
regions  are  covered  with  a  grid  of  rectangular  boxes  of  length  A*  The 
width  of  the  boxes  is  set  equal  to  <A  /  V®-1  thus  making  the-  box  diagonals 
parallel  to  the  Mach  lines.  The  source  strength  (downwash)  is  assumed 
constant  over  each  box  and  equal  to  the  value  ^at  the  box  center.  A 
typical  Mach  box  overlay  is  shown  in  Figure  (5.1.1). 

From  the  conditions .imposed  on  the  pianform,  the  downwash 
can  be  numerically  evaluated  for  the  diaphragm  and  -wake  boxes  and  the 
velocity  potential  profile  can  be  established  for  any  deformation  mode. 

The  pressure  can  then  be  calculated  from  Equation  4.1.8  arid  generalized 
forces  can  be  evaluated  with  Equation  2.0.1. 
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5.2  PROG&AIi  DESCRIPTION 

The  Supersonic  Unsteady  Aerodynamics  Program  calculates  generalised 
forces  for  up  to  10  deformation  modes .  The  computer  solution  is-  based 
upon  the  -Mach  box  technique.  The  various  configurations  which  can  be 
analyzed  are  shown  in  Figure  5.2.1  and  Table  5.2.1.  The  analysis  includes 
interaction  effect,  between  tandem  surfaces  and  wake  effects  on  the  trailing 
surface.  Single  surfaces  may  be  analyzed  by  inputing  a  second  Surface,  with 
a  zero  chordlength. 


The  supersonic  Mach  box'  method  uses  source  superposition  method  to 
approximate  the  aerodynamic  forces  on  an  oscillating  thin  planar  surface. 

The  diaphragm  concept  was  employed  to  handle  subsonic  leading  edges  and 
gaps  between  the  aerodynamic  surfaces .  For  purposes  of  calculating  pres¬ 
sures,-  it  was  assumed;  that  the  source  strength  over  the  area  of  each  box 
is  a  constant  value  which  satisfies  the  condition  of  tangential  flow  at 
the  center  of  the  box.  The  Mach  box  procedure  is  basically  the  same  as 
the  method  of  Pines,  et  al.  Reference  21 ,  differing  only  in  that,  the  surface 
and"  diaphragm  is  overlaid  with  a  grid  of  rectangular  boxes,,  the  diagonals  of 
which  are  parallel'  to  Mach  lines.  As  in  the  subsonic  and  transonic  cases,  the 
potential  equation  is  written  as  ah  integral  equation  (this  time  relating  the 
downwash  to  the  source  strengths)  and  approximated  by  a  matrix  equation.  The 
numerical  difficulties  are  primarily  ones  of  computer  logic  since  the  zones  of 
influence  of  a  given.  Mach  box  is  limited  to  the  region  within  the  aft  Mach 
lines-.  The  matrix  formulation,  leads  to  a  partitioned  form  since  there  are  two 
boundary  conditions  to  be  matched.  The  first  boundary  condition  is  the  downwash 
on  the  surface,  and  the  second  is  that  there  be  ho  pressure  difference  off  of 
the  surface  in.  the  diaphragm  regions.  The  zero  pressure  conditions  leads  to 
a  relationship  between,  the  diaphragm  potentials  and  the,  surface  potentials , 
and  the  surface  downwash  condition  then,  leads  to  the  surface  potentials . 

The  surface  pressures  then  follow  from  the  surface  potentials. 


The  solution  for  the  generalized  aerodynamic  forces  requires  the  input 
of  the  deformation  inodes  due  to  vibration.  The  program  considers  the  modes 
to  be  expressed  as  analytic  functions  of  the  form: 


w(x,  y) 


'(n-m)i 


(n-m).  m 

x  y 

m  J 


To  meet  this  requirement  only  the  coefficients  "c"  are  required  as  input 
into  the  program.  These  coefficients  can  be  obtained  in  several  ways, 
the  most  common  way  is  to  surface  fit  the  modes  by  the  least-square 
technique. 


DELTA 


T/?AP£ZO/0AL 


T/?AP£ZO!DAL 

CC&OPP£D) 


<t 


DELTA  PECTANGULAL? 

(.CHOPPED) 


E7GOPE  S.ZJ 

TAPDEM  COPLAJ'/Af?  COA/P/EOPAT/OHS  AT  SOPPESOaP/C 

MACH  A/<PMBE/P 


TA-iLli  5i2.1-  OPTIONAL  00I42CURA'JLT0N3 


CONFIGURATION 


RECTANGULAR 


DELTA 


TRAPEZOIDAL 


CHORBWISE 

COORDINATE 


TRAPEZOIDAL  (CROPPED) 


DELTA  (CROPPED) 


x(i)  «*  0.0 

X(2)  -  0.0 

x(3)  >  0.0 
x(:4)  >  x(3) 
-X(5)  £  X(4) 


X(l)  =  p.o 
X(’2)  >  0.0 
x(3)  =  X(2) 
X(4)  >  X(3) 
X(5)  *X(4) 


X(l)  =  0.0 
X(2)  >0.0 
X(3)  =  X(2) 
X(4)  ^  X(3) 
X(5)  ^  X(4) 


X(L)  -  0.0 
X(2)  >  X(l) 
X(3)  >  X(2) 
X(lf)  *  x(3) 
x(5)  2  x(4) 


x(i)  =  6.6; 
x(2)  >  0.0 
x(3,)  >  x(2) 
X<4),g  x(3) 
X(5)  §  X(4) 


SPANWISE 

COORDINATE 


y(i)  ■  o.o 


y(3)  >  o.o 


y(a)  =  o.o 

Y(2)  =  o.o 
Y(3)  >  O.o 


Y(I)  =  0.0 
Y(2)  >  Q.0 
Y(3)  >  Y(2) 


Y(l)/=  0.0 
Y(2)  >  0.0 
Y(3)  >  Y(2) 


Y(l)  =  0.0 
Y(2)  =  0.0 
Y(3)‘>  Y(2) 


5-3  INPUT  INSTRUCTIONS 


t-S 


Instructions  for  preparing  input  data  for  the  .'supersonic-  doippuiar  „ 
program  are  presented  here.  The  field  location  and  format  for  each  quantity 
is  specified.  Any  set  of  units  may  Ve  used  for  geometric  .dimensions  .and  acoustic 
velocity  as  long  as  they  are  consistent,  e.g.,  if  inches  is  used  for  length,  then 
the  acoustic  velocity  must  have  dimensions  qf  inches  per  second.  The  required 
data  and  the  sequence  in  .which  the  information  is  entered  is  as  fellow ss 
1.  Streamwise  Coordinates  (6E12.5 ,fd;mat) 


Column 

H2,  ... 

13-24 

25-36 

"  37-48 

49-60 

61-72 

Name 

Xd)  . 

X(2) 

X(3)  1 

..  X(4)_ 

X(5) 

Item 

(1) 

(2) 

-131—, 

_ 

_ (5)._ . 

(1)  X(l) 

(2)  X(2) 

(3)  X(3) 

(4)  X(4) 

(5)  i(5) 


Wing  root  leading  edge  coordinate  (see  "Figure  5.3.1) 

Wing  tip  leading  edge  coordinate 

Wing  trailing  edge  coordinate 

Control  surface  leading  edge  coordinate 

Control  surface  trailing  edge  -coordinate 


A  single  surface,  the  wing,  may  "be  analyzed  by  .setting  X(4)  and  1C (5) 
equal  to  X(3).  The  various  configurations  are  generated  as  .shown  in 
Table  5-2-1*  The  origin  for  the  planform  and  A1C  station  coordinates 
must  be  at  the  leading  edge  root  of  the  wing,  therefore  X(l)  and  T(l) 
described:  below,  must  always  be  zero. 


2.  Spanwise  Coordinates  and  Acoustic  Velocity  (6E12.5  format) 


Column 

1-12 

13-24 

25-36  , 

37-48 

49-60 

61-72 , 

Name 

Y(l) 

Y(2).. 

-  Y(3) 

.  SOUND 

RH0 

Item 

_1D _ 

Q  ‘  - 

(2) 

(3) 

-14Y_ 

(5)  : 

. 

,(1)  Y(l) 

(2)  Y (2) 

(3)  Y(3) 

(4)  SOUND 

(5)  RH0 


Wing  root  spanwise  coordinate  ■ 

Wing  leading  edge  spanwise  coordinate 

Wing  (and  control  surface)  tip  spanwise  icoordinate 

Speed  of  sound  at  altitude  for  which  analysis  is. performed 

Density  of  fluid  *  1000.0  ( ft. /I?) 
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3.,  General  Information  (6112  format) 


Column 

Name 

Item 


1-12.  . 

13r24  ' 

25-36= 

,  37-48 

49-60  _ 

61-72 

NMACB. 

NFREQ 

„ NMODES 

NBW  c. 

LYPIC 

„LSSVP 

■  (D  -  ■ 

(2) 

__I3j _ _ 

w 

(5),,. 

_!§)- 

CD 

NMACH 

(2) 

NFREQ 

(3) 

NMODES 

(4) 

NBW 

(5) 

LVPIC 

(6) 

LSSVP 

Number  of  Mach  numbers  (max  6) 

Number  of  Input  frequencies  (max  10) 

Number  of  input  modes  (max  10) 

Number  of  chordwise  wing  boxes  (max  10) 

Print  velocity  potential  influence  coefficients;  O'-'-'NOj  i-^Yes 
Print  upwashes;  0  ''■'No,  1  ^Yes 


4.  Mach  Numbers  (6E12.5  format) 

Column 
Name 
Item 

(1)  FMACH(i)  Mach  numbers  for  which  the  analysis  is  to  be  performed. 


1-12 

13-24  .  * 

25-36 

37-48 

49-60 

61-72 

FMACH(l) 

FMACH(2) 

FMACHC3) 

:  FMACHC4) 

FMACH(5) 

FMACH(6) 

(D 

(2) 

—13)- _ _  .. 

_I4) 

-  (5)  ~ 

(6) 

5.  Frequency  (6E12. 5  format) 

Column 
Name 
Item 

(1)  FREQ  Frequencies  for  which  the  analysis  is,  to  be  performed. 
Continue  on  next  card  for  FREQ(1)>  6 

6.  Deformation  Modes.  Repeat  the  following  Cards  NMODES  Times 

Column 
Name 
Item 

(1)  NfMl(i)  Number  of  deformation  mode  coefficients  for  the  wing,  mode  i 

(2)  NFI  Compute  generalize  forces;  0<~'No,  l^  Yes 

If  NFI  =  0,  the  program  will  compute  the  VPIC's  arid;  stop. 


(6112)  Format. 


1-12 

13-24 

25-36  . 

37-48 

49-:  60  , 

.  61-72 

NTMl(i)j 

NFI 

(D  i  (2) 

' 

1-12 

13-24 

25-36 

37-48 

49-60. 

61-72 

FREQ(l) 

FREQ (2) 

FREQ(3) 

FREQ (4)/ 

FREQ (5) 

FREQ (6) 

73’ 


Column 

Name 

Item 

(1)  C(5(i) 


Column 
Name 
Item 

(1)  NTM2(i)  Number  of  deformation  mode  coefficients  for  the  control 

surface,  mode  (i) 

(2)  NFI  Compute  generalized  forces;  0  »-»No;  l^*Yes.  If  NFI  =  0, 

the  program  will  compute  the  VPIC's  and  stop. 


.  (6E12/5)  Format 


Column 

1-12 

13-24  ■ 

.  25-36 

37-48 

49-60 

61-72 

Name 

C0(l)  „ 

C<K2) 

C0(3) 

C0(4)  . 

C0(5) 

C0(6) 

Item 

\  CD 

(2) 

(3) 

<4)  , 

(5) 

:  .  (6) 

C0(i)  i  -  1,  NTM  deformation  polynomial  coefficients  to  be  input 
in  the  following  order;  0,6;  1,0;  0,1;  2,0;  1,1;  0,2;  3,0; 
etc.  where  the  first  integer  is  the  power  of  "x"  and  the 
second  is  the  power  of  "y" .  Continue  on  successive  cards 
until  all  polynomial  coefficients  are  input. 


i  =  1,  NTM  deformat'ion  polynomial  coefficients  to  be  input 
in  the  following  order:  0,6;  1,6;  0,1;  2,0;  1,1;  0,2;  3,0; 
where  the  first  integer  is  the  power  of  x  and  the  second'  is 
the  power  of  y.  Continue  on  successive  card  until  all  poly 
nomial  coefficients  are  input. 


5.4  SAMPLE  PROBLEM 

The  generalized  forces  are  calculated  for  the  configuration  below. 
The  flight  parameters  and;  pertinent  (input  data  are  on  the  first  page  of 
the  computer  print  out. 

The  coefficients  of  the  deformation  modes  for  the  forward  surface 
are  shown  oh  the  third  page  of  the  computer  print  out,  and  for  the  aft 
surface  on  the  fifth  page  of  the  computer  print  out. 
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63.191  DEO 


PROGRAM  LISTING 


FORTRAN  NLSimi.DtGK  SUP/U210 

lm.nuF.  lR.-r  slip  /  o  2211 

CliK  IV  DHl  V  *  SUP  7 1)2  3U’ 

COMPLEX  CZF:  0.  VPIC.SS.nii.U.PHlW.SPKJyPHl.PJHTfc.GPHliFXF  SUPVUVri.li 

GHHHO^ZCl/Knnx(vnuO),XF(5),1fE(3>r;*{:i3),lXl,X2^X3^X3»»YA,Y2,BETA,HBS  SUP7  025U: 

rUHH0N/C2/ AS.  ‘jRaHh,  FHaCH(6> » NFR£Q»  FRECi  Mil ) . NHOBE. NSURF *LjfPlC»L-SS¥PSUP7 02.60 
GUHMON/C3/CI  MU  »2d. 2  ) » NT  M0.2 ) » NF (1,0. 2  ?» NT MAXI  2 ) . PXY (28 ) if  MC2U)  SUP/ 02? U 

rfiMH0rl/C4AV»'H.(  l275)»SS(2U0U).I)U{28.21.sU(10»:lU»3).PhIM(S#).SPMI  SUP/b(i5u 

nilHhON/C5/H..R{5  i )  . NHL (50 >. KC (50 >,KL (2b T.BSL (20). DXE< 7 1.TP1.U  SUP7029U 

CUMKON/C6/X. Y, h  <,  iiY . EH. EK.FKB.EKR.NP. HP. NB. NBUX. XODh. H0UE. NbW. NBT  SUP/0  3UII 
CUMMON/C7/C7EPO,  PH  I ,  PH  1  TE.  DPH 1..  XL.  NS.  NTH.  K.  J.  UR.U  i.  UAj.  UAN.  1  Frt,*  T  ULSUP/031  ti 


Poll  I  (2) 


•SlN(ARG)) 


CUMMOn/CB/Ri'O 
CAM.  DA  IN 

111!  1  n  It  0  HAC«  =1  .  IMACII 

FM-I<MACH1  HAI  II I 

IF  (l-M.LT. 1 .1  )  on  To  1000 

CAt i  COoL 

J 11 H.  -  T hL/Hp  T  A 

ca!i  iMn*ra) 

0  =  AS*FM 
TIM<  =  T  |»1/  II 
RFF=f>X»(Fh/t  F  1 A  ) *  =  2 
lili  VOG,  iM<  -Ti.NMtrU 
FK  =  l  Rpbl  I  FR  )  *  IPO 
FKH  =  (FK*KF* 

FKP  =  FK*Xl/?.li 
CAI  I  CAf  1 

I  F  <  |  VPjC.NE.ll  J  U/il  I  Poll  1  €  2  ) 
•A'Mi:=EK*  «* 

FXF  =  CMP'LX{C«‘S(  A’((}),-S|N(ARG) ) 
no  ».nn  mliiE-1  .NiOuF 
no  in  j=ii5f. 

-io(  j.-i)sczE>«n 

II-  (NF(MOllF.l).F...  |)  go  TO  210 

x=n.5*nx 

Nll-1 

JF(4  SSVP.NE.O  )  GALL  nuUT<(3) 

no  c ii u  hP=i>ni«o*. 

K II  =  |(OOX(NI  )' 

NS  -  l 

GO  10  (  7 0  *  6i’,.  7 0  >  6 «!.  on  ,70 »  70  ) »  KD 
NS  -  2 

MM  =  MOH(NP) 

Y  =  0  ;  IT 

1)0  Kill  MP=1,M«) 

HP  =  MP 
KOHI  sKUUXlNI*  ) 

SPH1=C/C«0 

IF  (OP. ‘GT. DfAi  l  P-llo 
SPH  K=SPII  I  * li > 

P  H  1  =  C  Z  F  H  0 

GO:  i  n  <  4  U  ,  4  n  ,  4  TP,  4  >l  ,1 0  .  ?  0 , 3  0  )  ,K 01 
FPH  l  =SPHl -PM-ii  (  -IP) 
ph I =PH I W( OP  ) 

PII]l(lMP)  =  l,HIW(MP)»FXF 
GO  TO  50 

IF  (Kll.l  1 .6)  GO  10  40 
SS(=NB)=-SPH'l/VPin/n-Y 
fill  10  911 
CAl.l  WASH 


P|JHT(3) 


4  <1 , 4  0  .  ?  0  »  3  0  )  ,‘R  Oil  E 


SUP/032U 
SUP /ill  33/1 
SUP/ 031  ll 
S  (JP7  6  3  50 
SOP/ 0  36  0 
SOP /U 370 
SUP/aiSrtO 
-SliP/ii.iyii 
SUP/U4UI) 
Si.P/lMl.U 
SnP7cii4.?n 
SOP/ 11430 
SOP / 04 4 U 
SUP/ir4«>o 
SUP /O 4 60 
SUP,/ 0  470 
SUP/048U 
SUP/ 0490 
SOF/OSOU 
SUP/  (IH  2  II 
SUP/nS2li 

Sup/ 0530 

SUP/11&40 
Sup/nb5ii 
SUP/ 0560 
SUP 7057 (i  ! 
SUP7II580  - 
SiJP7  059li 
SUP  / 0  6U(j 
SUP7063  d 
SUP/ 0  6?.,d  ' 
SUP /l)  63  d 
SUP70640 
SUP / 0650 
SUP/ 0^60 
SUP 70 670 
SUP/ 0630 
SUP/ll69r<j 
SlJP/0  70  u 
SUP/0,7111 
SUP 7 0/20 
SI :  P  ?■  II 73  ll 
SUP70740 
SUP  7 0  7b U 
SUP7U/60 
SUP 7 077 n : 
SUP 70 /HO ; 
SUP  7  0/19  0 
SUP  7  0  rill  u 


88 


;  ft 


to 

Ho 


J.F(KD<»L1<»6.}  SS.fNfi)  =  SS(Na)-.ARLKfTqB5*(5S(KH)+SPHI/VplC/Dr) 
lF;(*OD£.Ufci6)  R.u  10  90 
P U4-.SPH  1;+  SS  (  NH  )*  V  PI  O/J  Y  ' 

I F  ( KOMF.fc  0 t  A  )  pH!  H  ( MPj=PHl*bXF 
jrCkNP.EO.NRlvx-j  J  ZH1X.(!HP)=PM 

J.y  ( MPVJEO.NBIIX )  °  HI  Tit -P, HI  +:(PHI-PKT  k  (  HP )  JtDXJE  (5 ) 

CM  b  DOI J  i 

•1FCI  SSVP.NE.T)  )  CAl-t  POUT (4) 

NB  =  NB  +  1 

ki>  »  xoue 

Y-Y+DY 

x=x+nx 

D  (I  4  0  0  HO  -n  ,  NH-lML 
DO  ,30  0  NS=1  i  NSUTF 
f)(‘MUPF.,M0,N5:>=r.7ER0 
HTM=  NH-MO,  NM 
no  3oo  u-i.’-ia 

0  ( none,  HOj.NSI-O ( MuDj-.,  MO,  NS  )  +  CO;(HQ,  N.»  NS  )*DQ( N, NS> 
0(MUDF,'H0,3)  =  C2-R0 
Q  { MODE,  MO,,  3  J=U  (  HOilfc,  MOi  1  J  +  Q  ( MODE.-,  MO, .2 ) 

C  QNI INIIR 

OF  =  B.OfllXMlY^HHO 
D(J  HOfl  NS=1,3 
CAI.I.  pnilT(6) 

DO  7  0  0'  J=l,UMi5W: 

DO  70  0  K=i,  f.MfiDt 

OR  -  uMRBAl  (•!(*,  J,  US) > 

oi  =  'jF*ATMAO(n(K,  J,  NS,),} 

GAP=S'JR.l  (  U  ]  t-n  [  +UR*>IJP  j 

OAN=0^n 

IF  ( (jAM.'NE.  0 .  tM  »A  J=57*,?9578»A  I  A'N2(QJ  ,CRJ 


SUP70  61(i 
S  UP  70  M2, 0! 
SUP/D630 
SUP  7:0  8  4  0? 
SUP7  085,0 
S UP 7  0  86;S  • 
SUP7'0  870  ; 
SUP /  O  ft  8(0  = 
SUP/0890 = 
SUP7090I0  | 
-SUP 7jl '910  ' 
SUP  7  09.20  I 
SUP70  9“3i)i 
SUP 7 094 O' 
SUP  7.0  5*5,0^ 
SUP/ 0960 
SUP 70V? U 
SUP 7 098 6  ; 
SUP7II9V0 
SUP  71:0,0  0 

sup  )u\  ro  i 

SUP  7H)  ?  0  ; 
SUP 71 030  < 
SUP 71040  < 
SUP 710 5(0  l 
SUP71-0  60  " 
SUP/1.070 
SUP 71080  i 
SUP 71 090 
SUP  7110  0 


'40  Tn 


8  0  0 
.  9Utfl 
10  on 


% 

$ 

CUAiJ  N 


IF  ( OAM .t/E.  0 .  .<A  1=57  *,29578*  A  I  A‘N2(0i  ,45  R.)  *UP71-110 

TALI.  POUT  17)  .5 IIP/ 1.1 20, 

,IF (  NSiJRF TO  Tn  900  SllP'?il3U 

CONTINUE  SOP71140 

C0NIl4uF  .  SUP  71150, 

Ci.mil  N II F  SUP  7116  0 

GO'  fO  1  SllP/li.70 

SND  SUP71 ISO 

FORTRAN  NLSTiil),  Dl.-CK  , SUP/1190 

INCOME  J  H^F  SUP71200 

MAIN  SUP712lli 

SUBROUTINE  ■  ■  AT  N  SUP71220 

COMMON/Cl/iK<*OX (-7  0 <!.():)  >  XE  <5  )„  Y,E  <3  ) »  AR (3  > ,  XI, X2»X3»  X4,  Y(l  ,  Y2»  UETA, NBS  SUP7123U 
G0MM0N/02/.AS.  TMACn,  F  MACIi(  6)  >  NF  RF  0 »  FREtTlO  ) ,  NMODfc,  NSURF,  LVPI C,  LSSVPSUP/1240 


NlSTiill,  Dl.-CK 
J  H^F 


SUP-/ 1.12.0, 
SllP?il30, 
SUP71140 
SUP  71150, 
SUP 71160 
.SUP/li70 
SUP71180 
, SUP/1190 
SUP71200 
SUP71214 
SUP7122U 


C 0 M U 0 H / C? 3 / C 0  ( 1  0, 2 d , .2 ) IN T  ( 1 0, 2 .) ,  N F  ( 10, 2  ) „  N T M  A X  (  2 .)  >  PXY  (  28.) ,  F,N  ('2.8')  S 1 1 P  7 1 25  0 


COMM0N/C6/XUY,nx,  1) Y,  FM,.EK ,EKB, FKR,  NP , HP ,  NB,vNB4X,  KODE ,  MODE*  NSFF,  NBT  SUP7.1260 


COMM0N/C8/RH0 
R  E  A  l.i  (  5  v  11 )  <Xh'(  I  h,  1=1,5  ) 

RE  A  U.(5 ,  il )  (  YE  (I .)  •  .1=1  ,3),AS,RHO 
RHO  =  RHO/IUOU.O 

REAM  (  5 „12.)  * 'M AC > f>  N FR FO,  NMODfc ,  N HR.,  L V P  I C  ,  I. S SVP 
R.EAU4  5,11.)  (  FMAUH  (  f.)lsl=l ,  NMACH  ) 

REA:m  (  5  .ill.)  (  F  RF  JU  ),  l.=l  ,/N)  REO  ) 

NSIIPF-2 

I  F  (  XE  (4  )>.  L  i‘.  XF.  ( 5  ) )  i;0  T(R  10 


I  F  (  X  E  ( 4  )>.LI‘, 
NSURF=i 
XE(4:)=Xh(  3  ) 
X'F  ( 5,  )='XE  (3  ) 
NTMAXCl  )-fJ 
K'TMAX(.2)=0 


SUP 71270 
SUP71280 
S  UP7129IR 
SUP  7130-0 
SUP/1310- 
SUP n 320 
SUP71330 
SUP? 1340 
SUP  / 13  5.0 
SUP7136U 
SUP  71 37  0 
SUP  7 13 8  0 
:S  UP  7 139  0 
sup/ J  4  0  0 


c  •  1 


89 


nn  30  HllUE=1,MMGDfc  SUP71410 

no  on  I  =1 , N'-llRF  SUP7142U 

DG  20  J=1  / 2>-  -  SUP71430 

2n  COGNODF.J,  1  )  =  «..)  SUP7144U 

RFAIi ( *i » 12  )  Tm.  -iF  I  SUP714&U 

NT (MODE. I )=•  Tm  '  '  SUP71^6U 

MHAX(  I  )  =  MAxn(NTHAX(  IV),NTH)  SUP714J0 

NFfhOiiF.  I  )  =  i  F I  f  SUP71480 

30  RE/Al:(3,ll  >  (C»V<dQ|)F,  J,  I,),  J=!,NTM)  SUP7149U 

i  i  FORMA  F  (  6fcl2 . 8  )  ’  SUP71J>h  U 

1?  FORMAT (6 112)  .  SUP7151U 

RETURN  SUP71520 

END  SUP73 530 

•5  F OpTiRAN  NL^roiJiPcnK  SUP70030 

*  INCOME  IRnF  SUP7004O 

CCCN?  CONS  SUP7005U 

BLOCK  PaTA  SUP/ 0060 

COMPLEX  CZE^O,  VPfi;,SS,mj,(),PH!W,SPHl,PHi,PHlTt  ,  UPHI  SUP7007O 

COMMON/C3/C«i(in/2«»2  j, NT (10,2)  s>NF (10,2) i NTMAX < 2 ) , PXY < 28) , FN < 28 )  FUP/U0SIJ 

COMMON/ ij4  /  VI*  I C  ( 12  75  ) »  SS  ( 20  0  0  )  >  DC  ( 28. 2  )  ,  U  ( 1 0  >  10 , 3  )  »  PH  j  W  ( b'Q  )  ,.SPH  I  SUP/5u5o 

COMMON/ C5/MUB  { 8»i=):.  NHL (50  ) ,  KC ( 50  )  , KL ( 2b  ) , BSL (20 ) , DXE (7 )  ,1P>1,  U  SUP7010U 

COMMnN/C7/CZFpri.PJlI,PHl  IE,DPHl,XL,NS,NTK,K,  J,UR,QI,.UAB,UAN,  IFR,  THLSUPZllilU 
COMHON/CH/RcO  SUP70120 

OAT/i  KC/1,2»4, 7.lL,»16,22,29,37,46,56,6?,79,92,lU6,12t,l37,154,l'72,SUP7U13u 
3 191, 211^232,254, 2/7, 301,326, 352, 379, 4 07, 436, 466, 497, 529, 562, 596,  SUP7G14U 
2631 ,667, 7 1)4, 74?,  701 ,821, 862, 904, 947, .991, 1036, 1062, 1129, 1177, 1226/,  SUP7  0150 
3FN/0  .;  1 .-» I)  i  .  2 .  ,1  ..0«.3..2.,1  • 9  •  ,4  3  »,  2  ,  ,  1  • ,  [i 5  •  ,  4  •  ,  3  *  ,  Z  •  ,  1,  .l1 «  ,  SUP7  0160 


4  6 .  ,  5  . ,  4  . 


?.,3  .  >  0<./-»  TP  1/6. 2831853/ »C2ERO/(  C.U,  0,.0:)/ 


SUP 7 01 70 


1 

f 

pcopf 


DATA  KL/1,1  ,<Vi2.3,l,4,5,6,l,7,8,9,in,l,li,12>13,14,15,l,16,l7,18i  SUP7M8U 

1  19,20,21,1/  SUP 7 0] 90 

PNU  SUP7U20U 

FORTRAN  NLI.ToUi  IK  Ck  *  SUP71540 

.PNC  OPE  1FKT  SUP71550 

CODE  SUP7-1560 

SORPOtJTINE  CODE  SUP7157U 

COMPLEX  CZE'-  0  j’VP  J  C,  SS,  PH  I H,  SPH 1 ,  PH  I ,  PH  1  tl ,  UPHl  SUP71580. 

C0MHON/Ci/K( nx(20u0),XF(5),YE(3),AR(3),XliX2,X3,X4,Yl ,Y2,BETA,NbS  SUP/1590 
C0MM0N/C2/ A '■ » .JMAC'f ,  f  rt ACM  ( 67  ,  KF  RfcQ,  FHEC  ( 10  ) .  NMODt,  NSURF,  LVP I C,  LSSVPSUP7160  0 
C0MM0N/C3/CU { 1 II ,  28 ,  ?->) ,  NT ( 1  0  » 2  ) »  NF  ('10 , 2  ) ,  N  TM  AX  ( 2  ) ,  PX  Y  ( 28  ) ,  F  N  (  28 )  SUP/ 1610 

C0MMnN/C4/VP!C<  1275  ),SS(20  00),  00(28,2),  0(10,10  i3),  PHI  H(  50  ),SPHi  S0P75(l5ij, 

CUMMGN/L5/Mi.B(5  L) /NHL  (50  ) .  KC  ( 50  ) ,  KL  <  28  ) ,  BSL  ( 2>0  ) ,  DXE (  7  ) ,  T-P I ,  U  SUP/lft3u 

QOMK0W/C6/X.  Y,  fix,  )Y,  EM,EK,EkB,EKR,NP,MP,NIT,  NBUX.  KUDE,  MODE,  NBW,  NUT  SUP7,lf>4  0 
CUMi-  0N/C7/CJF  nO ,  P.i  I ,  Ph  IT  E ,  I1PH  i  ,  XL ,  NS ,  NTM,  K ,  J ,  uR ,  0 1 ,  OABv  UAN,  I FR,  IWI.SUP71650 


COMHON/Cb/Rl.O 

BETA  =  SGRTUj'FM  #  EM) -1.0) 
XI  =  XF  <  3 )  -  iXF(1  ) 

X2  =  X E <  3  )  -  X  E  ( 2  > 


X  E  <  3  ) 
XE  ( 4  ) 
XF  (5  ) 
XF  ( 5  ) 
YE  ( 2  ) 
YE<3) 


XF  (1  ) 
XF  ( 4  ) 
X  f:  ( 1  ) 
XF(1  ) 
Yf  (1  ) 


IF (X2.GT .XI. OH. XI  .Gl .X3. UR. X3.GL.X5. OF. X?.L T.O.U)  GU’  TO  5U 
if  (Y1  .GT.  Y2.0R.  n.LI  .0.0  )  GO  TO  50 
TUI  =  0.0 

IF(YP.Nt.Yl)  T  41  =  (XI  -  X2 )  /  <(Y2  -  Yl) 

A R  ( .1  )  =  (f?a(X2  +  XJ)  -  Yl  *  (  X2**.X1 )  )’ 

A  R ( 2 )  =  Y2»x4s?.0 
A R ( 3)  =  AK(1  )  ♦  AK ( 2 ) 

1(1  DX  =  XI/ ( F  L »  VH  <FH)  -0.5) 


SUP/l66lf 
SUP/ J  67  0 
SUP74ftfi(l 
SUP7169U 
SUP/ 1700 
SUP71 71U 
SUP7i72li 
SUP/ 1730 
SUP 7 1740 
SUP/ 1750 
S.UP/1760 
SUP7il770 
SUP  7.178  0 
.SUP/ 17  9  0 
SuP'7100  ll 
SUPZlHlu 
S.U  P'7 1 8  2  0 


90 


IF  150.0  *  HX  <r> G T #  ;X5 )  GO  TO  20 
15  NB7J  =  NBW.-l 
GO  70  10 
20  DY  =  DX/8ETA 
YNJ  =  Yl/UY 
YN2  -  Y2/DY? 

XNI.  =  YN2  -  (X1-X2)  /  DX 

XNT  =  YN2  ♦  X5/0X 

XNU  =  X3/DX 

XNTF=X5/DX 

NBOX=XNTE+0 .5 

NBS  =  Y2/DY  +  1.6 

NBT  =  X4/DX  1-0.5 

OXF(l)  =  l.n 

RXF(2)  =  1.0 

0XF13)  =0.5 

DXF(4)  =  AlnTtXMLt  ♦  1.5)  -  XNLfc 
fi  X  F  ( 5  )  =  X  N  Ui  -  Fi:  0  A  T  ( N  H  0  X  -1 ) 
nXF(fi)  =  0.:< 


l)XF(7)  =  0.0 
X  =  0.5  *  DX 


l]IJ  40  NP  =  1 ,  nbux 
XN  =  FLOAT (KP)  -  U.5 
YW  =  YN2 

IF  (THL  .07.  p.«ij  Y:W=A«INi(YHyYNl?*Xh/(THt/BFTA>) 

I F  (  X  .  U  J  .  Xb  \'V  )  )  GO  TO  24 

KH  =  HINK  AhAXI  (Yw,XN+YNl)yXNt-XN)  +  l- 

GO'  10  ?«  • 

KB  =  MINKAnAXI  (  XML+XN,  XNi+ YN1  >  ,  XNT'-XN ) +1 
MOB(NP)  =  MB 
KODF  =  1 

IF  CNP  .El).  NOW)  KOBE  =3 
IF  ( M SURF  .EQ.4)  00  TO  29 
IF  (X  .IFT.  XI)  KOBE  -6 
IF  (X  . G7 .  X.l  )  KODE  =4 
IF  (X  . G.T.,  X3  +  DX)  KODE-2 


IF  ( M SURF  . F Q . 1 )  00  TO  29 
IF  (X  .IFT.  XI)  KOBE  -6 
IF  (X  . G7 .  X.l  )  KODE  =4 
IF  (X  . G.T.,  X3  +  DX)  KODE-2 
IF  <NP  .ElJ.  NROX )  KODE  =5 
29  lF(NB+MB.GT.?OnU)  GO  TO  15 
'NBI.(NP)  =  NP 
1)0.  30  MP  =  3  .  MR 
YN  =  MP-1 
NIL  =  NB  +  1 

IF  X  YN  .GT.  Yw)  KODE  =7 
3 1)  KHOX  (kill  )  =  KOUE 
40  X  F  XfDX 
RFTliRN 

5l)>  CAT  1  EX  I  T 
RETURN 
FNI) 

MikTKAN  NLf.TfMJ.  DECK 
I MCOOE  IR.’  F 
f  I  CAFrl 

SUBROUTINE  OAF) 

COMPLEX  OZEPO,  VPIC,SS,no,«>PHIM,SPHl,PHl,;PHI7b»0PHI. 
DIMENSION  P<5),W(5) 

cnMhON/C4/VFIC(12^5),SS('2UOO  )»D(j(28'i2).a(10,lU,  J),PHIW(50) 
C0MM0N/C5/M0BI5IU) ,  NHL  (50  ) ,  KC  ( 50  ) ,  KL  ( 28  ) ,  BSL  (  20  )  ,  DXE(  7  ) ,  7  P I 
*C0M.f0N/C6/fcvY,nX,  i)Y,£M,  Eft,  EKB*EKR,NP',  PP*NB,  NUUXi  KODE,  MODE* 
C0MM0N/C7/C7EH0»  PH  I »  PH  I  rE->  DPH  i ,  XL,  NS ,  N TM ,  K ,  J »  UR,  Q  l ,  OAB,  UAN 


,  SPH  1 
,U 


Sl»P7i830 
SUP7184U 
SUP7185U 
SUP71860 
SUP71B70 
SUP71380 
SUP71890 
SUP7190U 
SUP719i0 
SUP/1920 
SUP71930 
§IJP7d94  0 
SUP7195U 
SUP7196U 
SUP71970 
SUP7198U, 
SUP71990 
SIJP720  0U, 
SUP720H), 
SUP72020 
SUP/2030 
SUP72G4U' 
.SUP7205U 
5UP/2G67I 
SUP72Q70 
SUP7208U 
SUP72G90 
SUP7210U 
SUP7211 0 
SUP72120 
SUP?2i3U 
SUP/214U 
SUP7215U 
SUP7216U 
SUP7.2170 
SUP72180 
SUP/2190 
SUP  7  2.2  0  0 
SUP/ 221.6 
SUP7222U 
SUp7?230 
SUP72240 
SUP72250 
SUP/2260 
SUP72.270 
SUP72280 
SIJP72290 
SUP  / 230  iJ 
SUP7231U 
SUP  7232l! 
SUP72330 
SUP72340 
SUP7235 1) 
SUP7236U 
SUP/2370 
SUP72381I. 
SUP7505U 
SUP/240U 


NOW, NOT  SUP72410 
,  IFR, 1WLSUP72420 


COMHON/G8/RhO 

DATA  P/ 0.95308992, 0.76923465,0 <>5,0.23076535, 8. 046918  03/ 

1  ,  W/0.11846344, 0.23931434* 0.28444444, 8.23931434,  0-.  11846344/ 

PI  =  TPI/2.H 

I F  (  PKB.GT.  0  .  O')  GO  TO  10  ■ 

VPIC  =  (-1.0, 0.0) 

GO  TO  30  , 

10  VPIC  =  CZERfl 
DO  20  I  =  1,5 
ARG ■=  EKB*P ( I) /2s  U 
F  =  -(0.5*APG/E4)**2 
ZJ  =1.0 
FI  =  1.0 
AE  =  1.0 
DO  15  K  =  1,20 
AE  =  AE  *  F/F 1**2 
FI  =  FI  ♦  1.0 

IF(A8SCaE).I  E.1 .E-5J  GO  TO  20 
1-5  7 J  =  ZJ  *  kb 

2 II  Vfh«  =  VPIC  -  7  J*W  (  1  )  »CMPLX.(  COS  (  ARG  )  ,  -S I N  (ARG  )  ) 

3n  no  oo  NP  =  2.NBOX 
K I  i  KC(NP) 

KZ  =  kC(NP  +  1  )  -  1 
ISO;  4  0  K  =  K  l  ,KZ 
4 it  VPll'(K)  =  C/FRO 
NO  =  NP  *  L 
Jlil  i'  0  I  =1,-5 

X  =  F  LOA  f  ( N;: )  -  0.5  +  P(I  ) 

APO  =  EKB»X 

PHI  =  U(  1  )*i.‘MPl  X ( -CoS ( ARG ) , S I N ( ARG ) ) *2. 0/P  I 
CAM  H'SIf  S(  ARG/F'i,  4  ) 

K  =  KC(NP) 

DO  7  0  HP  =  3  .HU 

FOX  =  (FLOAKrtP)  -  0.5  )/X 

C  =  SQRfll.O  -  tOX**2) 

AF  =  2.U*ATAN(E'JX/(1.0  ♦  O) 

S  =  2.0*EOX»C 

C  =  2 . 0*C*C  -  1.0 

SO  =  0.0 

VIN  =  BSL*AI 

F  =  1.0 

FI  =  t.O 

DO  50  L  =  1,N 

VIM  =  BSL(L+L)«S/f I  -  VIN 

SN  =  2.0*S«C  -  SO 

SO  =  S 

S  =  SN 

T  =  -F 

50  FI  =  FI  +  1.0 
DPHI  =  PWl»VIN*F 
VPIC(K)  =  VPIC(K)  ♦  OPHI 
VP!C(KU  >  =  VP  I C ( K ♦ 1 )  -  DPHI 
IFCHP.E0.1)  VPIC(K)  =  VPIC(K)  ♦  DPHI 
70  K  =  K  +  1 

80  VPIC(K)  ='  VPIC(K)  +P  T»BSL  *PH  I  /  2 . 6 
RETURN 
FNP- 

T  FORTRAN  NL  ^TuU, DECK 

H-  IMCODF  I R '  F 

r.OSLS  BSl.S 


SUP.72430 
SUP72440 
SUP72450 
SUP7246U 
SUP72470 
S.UP72480 
SUP7249U 
SUP72500 
SUP72510 
SUP7252C 
SUP72530 
SUP72540 
SUP72550 
SUP72560 
SUP/2570 
SUP72580 
SUP72590 
SUP72600 
SUP-7263  0 
SUP72ft2U 
...  SO  P'7  2  630 
SUP7264.U 
S0P7265U 
SUP72660 
SUP7267U 
S0P/268li 
SUP72691 
SUP72701 
SUP/2711 
SUP/2721 
SUP72731 
SOP/2741 
SUP72751 
SUP7276I 
SUP7.277I 
SUP72783 
SUP  727  9.1 
SUP7280 1 
SUP728l‘l 
SUP7282I 
SUP-7283' 
SUP72B4 
SUP/285. 
SUP7286i 
SUP7287 1 
SUP7288I 
SUP7289I 
SUP7290I 
SUP7291 
SUP7292 
SUP7293 
SUP7294 
SUP7295 
SUP729& 
SUP7297 
SUP729B 
SUP7299 
SUP7300 
SUP/3U1 
SUP7302 
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C 

il  « 


SUBROUTINE  kSLSLAKG,  N) 

C0HhfiH/C5/Hi‘B  (  5  B )  #  NftL(5U  > ,  KC  (.5 0  ) ,  KL  (  28  ) » 3SL  (  2  0  ?  *  DXE  ( 7 ) , TP I , U 
COMMON/ C8/RO0  - 

DO  1  1=1/20 

1  8 SI  ( 1 )  =  0.0 
ASQ  =  ARG**2 
I FTASU.LT. 0 i 01 )  GO  TO  50 
N’  =  HINl(17.n,  (ARO  +  10.0  )  j  ' 

F  =  ?*H  +  4 

BS1  (H+2)  =  r4.;6*r*(Frlo0>)/ASQ-(F-1.0)/F)*0.3E-30 
PF  =  0.0 
J  =  0 

no  10  j  =  j,n 
H  =  N’  -  I  +  1 
F  -  2*H  +1 

BSL  <M).=  <4.0*(F-t.U  )/ASQ-1.0/F-1.0/(Fr2.0)  )*8SL  ( M+l  )-BSL  ( M+2  )/F 
10  Pf  =  PF  +  2*  0* ( F'-2 .  U)*BSL,(H*1 ) 

PF  =  PF  *  BSU(f) 

F  =  0.0 

IF(ABSt(PF).GT.1.0)  F  =  A8SCPF  )*1.E-10 
N  =  N  *  2 
no  30  I  =  1,N 

IF  (F.GE. ABS(BSI.f(  D  )  )  GO  TO  20 
BSL ( I  )  =  BSl  (  I  )/PF 
GO  10  30 
20  BSL C I >  =  0.6 
30  CONTINUE 

nO  40  I  =  It N 

IF(ABS(BSL(N) ).GT.1.0  E-7)  RETURN 
40  N  =  N  -  1 
RFTORN 

&0  B.Sj:(  2  )  =/'  0 . 125*  A  SO 

BSl (I)  =  1V«  -  <>.U*BSL(2> 

N  =  2 

RETURN 

END 

7;  FORTRAN  NLST OU, DECK 

Sr  INCODE  IBMF 

CARLE  ARl-E 

FUNCTION  ARLE(T'iB) 

COH.HON/Cl/K}iOX(20U  0  ) , xF ( 5  )  , YE ( 3  ) ,  AR ( 3 > , XI, X2, X3, X4 , Yl , Y2 , BET* , NBS 
COMHON/C6/X, Y,nX,nY,EH,.fcK,EKB,EKR,NP,KP,N8,NB0X,K0DE,MUDE,NBH,NbT 
C0MM0N/C8/RH0 

IF OX-0 .5*DX.G6.X1-X2 )  60  TO  10 
YT  =  ( Y-Yl)/T)Y  +  0.b-(X/nx-0.5)/TUB 
XR  =  Yl*TOB 

YH  =  AMAX1  (  0. 0,AH.IN3(  1.0,  YT-1. 0/TUB)) 

Y.f  =  AM  INI  (J  . ft,»  AiMAXl  (  6 . 0,  YT  ) ) 

XL  ~  AMAX1  ( (i  ,/H  >  AiKi  N1  (3.0,  XR-TOB.)  ) 

XR  =  AM  I  N1  ( 1>.,i) ,  AHAX1  (  0 . 0  ,  XR  )  ) 

ARI.E  =  AMAXl  (.  H  .  V*  (  Y(*  (  XR  +  XL  )  ♦  YB»  (  XR-XL  )  ),0.0) 

IF(MP.EU.I)  ARLE  =  PEARCE 
RETURN 

in  ARIF  =  AMIN3  (1.U,  AMAXl/dl.O,  (  Y-Y2  )/DY  +  0 .5  )  ) 

RETURN 

END 

<•:  FORTRAN  NL'-T  jU  •  DfcCK 

I'H  OOF  j  H '  F 
CPhLP  Pul>i 

SUBROUTINE  fH  IB 
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SUP7303U 

SUP/3n4a 

SUP73H5U 

SUP73060 

SUP73070 

SUP73083 

SUP730$|, 

SUP73100 

SUP73110 

SUP7312U 

SUP73130 

SUP73140 

SUP7315U 

SUP73160 

SUP7317U 

SUP733.80 

SUP733.9U 

SUP  7,3  2  0  0 

SUP7321 0 

SUP7322U 

SUP73230 

SUP/324U. 

SUP73250 

SUP73260 

SUP73270 

SUP7328U 

SUP73290 

SUP7330  0 

SUP73O10 

SUP73321I 

SUP/3330 

SUP73340 

SUP73356 

SUP73360 

SUP7337  6 

SUP7338U 

SUP733V0 

SUP73406 

SUP7341(I 

SUP73420 

SUP73430 

SUP7344J) 

SUP/345  0 

SUP73460 

SUP/347G 

SUP73480 

SUP/3490 

SUP/35  0  0 

SUP/ 3510 

SUP/352U 

SUP73530 

SUP/3540 

SUP73550 

SUP73560- 

SUP73570 

SUP73580 

SUP/3590 

SUP7360U 

SUP73610 

SUP 73620 


m 


r  » 


in 

2n 


one  1  j 


211 


3n> 


5n 


J 

$ 

COEGE 


COMPLEX  CZEKG»VPIC,SS/DQ.*Q>PHIW>SPHI,PHI» PHI YE* DPHI  SUP/3630 

C0HM0M/C4/VP I C(  1275) ,  SS(2000  ) ,  DCl<28,2 ) »  Q(li,  18'.  3)  iPHlHTSS  )»SPHI  SUP7505U 
CGMMON/C5/MOB(50),NBL(5G),KC(5Q>,KL(28>,BSL(28>*DXEt7>iTPJ,U  SUP73650 

COMM6H/C6/X»Y,DX>nY,EH,EK,EKBiEKR,NP,HP»HB*NBOX»kObE,MbDE>NBH,NBT  SUP7366U 
CQHH0N/C7/CZ£R0,PHl,PHlTE,DPHI>XL,NS,NTM,K,J,QR,QI,QAB,QAHiIFRiTWLSyp7367Q 
C0MH0K/C8/RH0  SUP/3680 

DO  20  1=2, NP  ,  '  SUP73690 

NU=NP-I+1  SUP73700 

JL=MAX0  Cl » HP- I *1 )  SUP73710 

JR=MINO(HOB(NU)»MP*I-1)  SUP7372G 

Nd=NBL<NU)+JL  -  SUP73730 

DO  20  J= JL, JR  SUP/3740 

K=KC<  I  )+IABS(MP-J>  SUP7375Q 

DPH1=VPIC(K)  SUP73760 

IF  ( J.GT. I-uP*l.OR. J.EQ.l)  GO  TO  10  SUP/3770 

K=KC( I )+MP+ J-2  SUP/3780 

DPHl=DPHl+VPinU)  SUP73790 

SPH 1 =SPH I +  DRH I  *SS(  HJ )  SUP73800 

N J=N J*1  SUP/3810 

RETURN  SUP73820 

END  *  SUP73830 

FORTRAN  NLSTltU,  DECK  SUP/3840 

I MCQDE  IBrF  SUP 73850 

nolJ  SUP/3860 

SUBROUTINE  ROIJ  SUP7387U 

COMPLEX  CZEKO.VPIC,SS,D0,O,PHIM,SPHI,PHI,PHITE,DPHI  SUP73880 

COMMON/Cl/Kf  0<(2(liiO  ),XE15)i  YE(3),AR(3),X1,X2,X3,X4,Y1,Y2,BETA,NBS  SUP7389U 
COMhON/C3/CO(10,28,2)*N1 ( 10 , 2 > , NF < 10 > 2 ) , NTMAX (2 )> PXY (28 ) , FN< 28 )  SUP/3900 

C0MM0N/C4/ VP  » C<  12/5  ),  SS  (  20  0  0),BCl  (28,2),  Q(1D,1U,  3  ),PHIH(50)>SPH  I  SUP75U5U 
COMMON/CS/MOB  <  50  ) »  NfiL  ( 5f« ) » KC  { 50  )  ,  KL  ( 28  3  , BSL( 20  ),DXE(7), TPl»U  SUP73920 

C0MN0N/C6/X, Y>  DX, l)Y, EM >  EK , EKB, EKR, NP , MP> N8, NBOX, KODE, MODE*  NBH, NBT  SUP73930 
rOMMON/C7/CZERn,PHl,PHlTE,DPHi,XL>NS,NTM,K,  J,aR,Ol,UAB,OAN,IFR,TRLSUP7394'0 
C0KM0N/C8/.RU0  SUP73950 

NTM  =  NT MAX ( NS )  SUPZ396,U 

DPUI=PHJ  SUP73970 

1  F  <  MP.EU.l  )*  nPH  I  =o  . 5*l)PH I  SUP73980 

CUM  =  M*UXE<  K*»n»: )  SUP73990 

110  2  0  K  =  1  ».N  I  M  SUP740  0  0 

K  K  =  p  L ( K )  SUP74010 

DUdwNS)  =  *  Q  ( K ,  N  !> )  +  nPHI*CON*CMPLX(PXY(KK)*F-N(K)»-FK*PXY(K) )  SUP-74020 

IF ( K OIJF •  Nfc •  3 )  O' i  10  5J  •  SUP/4030 

110  30  K  =  1  >  NTM*  SUP/4040 

DO(K,l)  =  <l)i  (  K ,  1 )  -U»0PHI«PXY(K)/DX  SUP74050 

ARC  =  F K *  (  X .5  -  XI)  SUP74060 

DPHI  =  -  DPH I  *CM?L  a(COS(ARG)#tSIN(ARG))'  SUP/4070 

XL  =  0.0  S  UP  7  4  0  8  O' 

CALL  UEGE  SUP7409Q 

I F  (>K ODE .  NE . 5  )  RETURN  SUP741Q0 

XL  =  X4  SUP7$110 

IF(MP.EU.1)PH,!TE  =  0 .  5*PH  IT  E  SUP/4120 

DPHI  =  PHITE  SUP74130 

CALL  QEGE  SUP/414U 

RETURN  SUP74150 

END  SIJP74160 

FORTRAN  NLSTulI,  DECK  SUP/4170 

INCODE  IBMF  SUP74180 

OFGE  SUP74190' 

SUBROUTINE  UEGE  SUP 7420  0, 

COMPLEX  CZERO, VP  TO,  SS,DO.<J»;PHTR»  SPH  I ,  PHI ,  PH  1 TE,  DPHI  SUP74210 

C0MM0N/G4/VP1C(  l?  /  5.) ,  SS  ( 20  0  0  ),  DO  ( 28 , 2>) ,  Q  ( 10 , 10 , 3 )  ,  PH  I  W  ( 50  )  /SPH  I  SUP  75 1)  50 
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l 

y-.  ^ 


10 


COHMON/C5/KOB(50),NBL(5d),KC(50),KL(28)>BSL(20)»DXE(7j#'tFl*U  SUP74230 
CQHH0N/C6/X, Y» DX,  I) Y,  EH /  EK ,  EKB/ EKR/ NP .» HP >  NB >  NBOX.  KODE ,  MODE*  NBW»  NBT  SUP7424t) 
C0HH0N/C7/C2ER0.*  PHI »PHITE/DPH1* XLV  NS»  NTH*  K»  J#  QR*  01 / GAB/ QAN>  iFR, TWLSUP74250 
C0HH0H/C8/RH0  SUP74260 
PX=1.0/-DX  '  SUP74270 
IF(XL.E£J.0.0)  GO  TO  30  SUP74280 
YX-Y/XL  ,  SLIP74290 


K=1  '  SUP7430U. 
DO  20  N=l,7  '  SUP74310 
P Y=PX  SUP74320 
DO  10  H=1,N  "  *  SUP74330 
DOCK, 2)  =  DG(k>2)  -U*DPH I *P Y  SUP74340 
IFCK.EQ^HTM)  RETURN  -  SUP7435U 
K=K+1  SUP 7 4 36 U 
P Y=PY*YX  SUP/437 0  c 


2n  PX=PX*XL 
RETURN 


30  K-0 


SUP7438U. 
SUP/4390 
SUP 744 00 


HO  4(1  n=1i>7 

: K  =  K  +  N 

IF(K.GT.NTN)  RETURN 
PU(K,2)'  =  Hi' (  K  /  2  )  -U*DPH I *PX 
40  PX=HX*Y 
RETURN 
END 


SUP/441  U. 
SUP7442U 
SUP7443U 
SUP7444U 
SUP74450 
SUP74460 
SUP/4470 


FORTRAN  MLS T slH^Di-CK 


SUP7  4.480 


1NCODE  IBUF  SUP74490, 

HASH  SUP7450  0 

SUBROUTINE  RASH  SUP7451U 

COMPLEX  C2ERO,  VP  I  G»  SS*  DO»  U*  PHIH*  SPH  i  *  PH  PH  I  TE»  DPH  I'  SUP74520 

COMhON/Cl/KHOXf  2000  }/  XE ( 5 ) /YE ( 3 ) / AR ( 3 ) /XI /X2/X3 / X4>  Yl* Y2* BETA/ NBS  SUP74531? 
COMMON/C3/Cn(  HI  *2h»  2  )/NT(10/2)/  NF(10>2  ) »  NTHAX  ( 2 )  /  PXY(2B  ) « FN(-28 )  SUP7454  0 

C0H.H0N/C4/VP  I U  ( 1 275  )  >  SS  ( 200  tl  )  »  DO  (28»  2 )  *  0  ( 10 » 10  »  3 ) » PH  I  H(  50  )  /  SPH.I  SUP75050 

C0MM0N/C5/M0R  ( 50  ) ,  NUL  ( 50  )  /  KC  ( 50  )7  KL ( 28  )  /  BSL  ( 20  )./’DXE  (  7  )  /  T” I  /  U  SUP74560 


CUMM0N/C6/X. Y»  RXi UY, EH»  EK/EKB  »EKR*  NP 
C0HH0N/C7/CZERn>  PH  I » PH  I TE, DPH I / XL  *  NS 
C0MM0N/C8/Rh0 
XP  =  X 

Iff  NS.FQ»2 )  XP  =  X  -  X3 
PXY  =1.0 
DO  10  N  =  2,2« 

10  PXY(N)  =0.0 
TM  =  NTHAX (NS) 

II)  =  SORT  (2. 0*TM )  -  0.5 
I F ( IDiEQ.O )  GO  TO  60 
K  =1 

IFCXP.EO.O.u )  GO  TO  40 
YX  =  Y/XP 
PX  =  XP 

DO  30  N  =  1/11/ 

K  =  K  +  1 
PXY(K)  -  PX, 

4)0  20  M  =  1 ,  N 
K  =  K  +  1. 

20  PXY(K)  =  PXY (K-t )*YX 
30  PX  =  PX«XP 
GO1  10  60 

4 ij  px  =  l.n 

no  5  0  N:  =  i,  in 
PXY(K)  =  PX 
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MP,NB/NBOX/KODE,HODE/NBW,NBT  SUP74570 
NTH/K/J/QR/QI  /  GAB /Q AN/  1 FR7  THLSUP7 458U 

SUP74590 
SUP74600 
SUP74610 
SUP74620 
SUP74630, 
SUP7464U 
SUP74650 
SUP7466U 
SUP74670< 
SUP74680 
SUP74690 
SUP7470U 
SUP74710 
SUP7472U 
SUP74730 
SUP74740 
SUP7475U 
SUP74760 
SOP74770 
SUP7478li 
SUP74790 
SOP/4HOO 
SUP  7  4  Bllj 
SUP7482  0 


•*  -  -S' 


f 

r 

r 

[ 


m 


l 

l 

[ 

[ 

L 

L 

L 

L 

[ 

t 

l 

t 


5n 

60 


7  0 
an 


pout 


K  =  K  ♦  N  ♦  1  S0P74831 

PX  =  PX*Y  SUP74843 

NTH  =  NKHO.HE.NS-)  T  SUP74R5C 

2  =  CO ( MODE , 1 >  NS )  SUP74860 

DZ  =  0.0  SGP74B70 

lF(NTH.Pp.l)  00  TO  80  SUP74880 

no  70  K  =  2, NTH  ,  SUP74890 

7  =  2*  PXY(  K,'*C0(  MODE,  K,  NS )  SUP7498.0 

KK  =  KL(K)  SUP74910 

02  =  DZ  +  PXY(KK)*Cu(HODE,K,NS)»FN(K)  SUP74920 

SS(NB)  =  U*i:HPI  X(i)Z,EK*Z>  -  SUP74930 

RETURN  SUP74941) 

END  SUP74950 

FOkTKAN  NLKTUll.DhCK  SUP74960 

POUT  SUP/4980 

SUBROUTINE  POllT  M  nD  )  SUP74990 

COMPLEX  CZEhO, VPIC,SS,nU,U,PHIW,SPHI,PHI,PHITE>DPHI  SUP75000 

11 1  MENS  I  ON  CARDS  (25,50)  «i*****« 

0  I  MENS  I  ON  Sw(5,A), SURF ( 2 , 3 ) , CUD ( 7 ) , C ( 50 )  SUP/5U10 

COMMON/Cl/KPOX(/iOOO  ) , XE ( 5 ) , YE ( 3 ) , AR ( 3 ) , XI, X2, X3, X4, Yl, Y2, BET A, NBS  SUP/5020 
COMMON/C2/AS ,  UMACH,  FMACH ( 6  ) ,  NFRL  Q,  FREO  ( 10  ) ,  NHODE,  NSURF ,  L-V.P I Ci  LSSVPSUP/5030 
COMHON/C3/CCnO,2«,?),NT  ( iO , 2  ) , NF ( 1 0 , 2  ) , NTMAX ( 21 , PX Y ( 28 ) , FN (28 )  SUP/50411 

COMMON/C4/VPI C ( 1-275 ) , SS( 2000 ) ,  DU(28,2 )  ,0(40 ,10 , 3 ),PHIW(50),SPHI  SUP75050 
C0MM0N/C5/H<'B  ( 50  ) »  N|»L  (50  ) ,  KC  ( 50  ) » KL ( 28  ) ,  BSL  (20 )  ,DXE(  7 ) ,  1  PI ,  U  SUP75060 

COMMON/CA/X,  Y,nx<i  HY,EH,FK,EK8,EKR,NP,MP,NB,NB0X,K0DE,M0DE^NBW,NBT  SOP/50 70 
COMMON/C7/C/ERO,P>il,PHlTE,DPHI,XL,NS,NtH,K, J, OR, 01 , OAB, UAN, I FR, THLSUP/5080 


COMMON/CB/Ri-n 

DATA  (SW(1, 1 ), 1=3 .6)/26HHAP  OF  MACH  BOX  OVERLAY  ON, 


SUP 75 090 
SUP/5100 


1 

26HWING, 

TAIL,  AND  DIAPHRAGM 

9 

SUP75110 

2 

26H 

(S')  "  WING 

SUP/5120 

3 

26H 

(S)  -  TAIL 

9 

SUP/5'130 

4 

26H 

(,)  -  WAKE 

9 

s6p7514  U 

5 

26H 

(i)  -  DIAPHRAGM 

/ 

SUP/5150 

8HTAIL 


iiHHiNG  4  tail  / 
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n A 1 A  (SlIRF'd  .  Hi  I  =1,3  )/8HWIKG 
data  COU/lHS,lHt,  LHSMHS,iH*,lH,  ,1H./ 

GO  10  UO  ,20  ,  (IP  ,40  ,50  ,60  ,70  >,JND 

WRITE (6,11 )f  M,AS,RH0,XE(1),XE(4),X1,X4,Y1, Y2, Y2iY2,X2,X4,AR(l 
1  AR(2),NBW,i'B  r,  iBS„NBS, 


SUP/5160 
SUP’/ 51 7.0 
SUP/518IJ 
SUP/51  90 
SUP 752 00 

11  FORMAK3H1/ / /  ///3'i X, 43H|jAC/NAA>  HISSILfc  SUPERSONIC  AIRLOADS  PR0GRAHSUP7521G 

1  /  /  /  3  7  X ,  3  0  H  K 1. 1 R  H  T.  CONDITIONS  AND  OE0ME7RY/1H0//15X,  13HMACH  NUHBERSUP75220 

2  =,F8.6,4X,  1  A.ISPE>:D  OF  SOUND  =  F10.3,4.H  L/T  ,  4X,  4HRH0= ,  614 . 8//1H0/  SUP/5230 

X54X,4HWlNG,iax,  SUP7524U 

3  4HlAlL///2?X,loHL.E.  S  I  ATI  ON  (  L  )  *  2F22  *3//22X»,16HR00  T  CHORD  (  L  )  ,  SUP/5250 

4  2F22.3//  22XMAHL.E.  SPAN  ( L  ) ,  2F22 . 3//22X ,  16HT,  b;  SPAN  ( L ) /SUP75260 

5  2F22.3//  22X, 3  AH  f I P  CHORD  ( L  )  , 2F22 . 3//22X, 1AHTOT AL  AREA  ( L»L ) , SUPV5270 

6  2F22.3//  22X,16HCH0Ri)HISE  bOXES  ,  119,  122/72 2X, 

/lAHSPANKTSE  BOXES  ,119,122) 

W  R I  T E  ( A ,  1 2  )  UB  0 X »  D.X ,  ll',Y 

1?  FORMA  f  MHO  / » 11  X, 23 HT OTAL  CHGRDWISE  BOXES  =,13,  5X,11HB0X  CHORD  =, 

1  1  P3  El'2.5, 2m  L,  5X,inHBOX  SPAN  =,  1P1E12.5, 2H  L/  ) 


13 


WRITE(6,?1) 

NB  =  1 

no  17  NP  =  ,1  ,  NB-iX 
MB  =  MOU(NP) 

IF ( MB . OTo 50 )  0)  10  80  0 

DO  13  HP  =  i,MB 
K  =  KBOX(NB) 

C(MP)  =  CUD(K)' 

NR  =  NB  +  i 

IF  (NP.  0^7.  A)  fit)  TO  15 
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SUP75280 
SUP75290 
SliP/530  0 
SUP75310 
SOP75320 
SUP75330 
SUP75340 
SUP75350 
SUP75360 
SUP 75370 
SUP7538U 
SUP75390 
SUP75400 
SUP/54130 
SUP/5420 


WR!TE(6,14XSHU,NP),I=l,5)i  <C(HP),HP=1,HB> 
FORMAT (iQX*5A6»50Al ) 

GO  10  17 

WRITE(6„16)  (C(MP),HP=1,H8) 

FORMAT! 40 X; 50 A1 ) 

CONTINUE 
GO  TO  1000 

NRIlE(6,8fll)  1 


SUP75430= 

SUP75449 

SUP75450* 

SUP75460’ 

SUP75470 

SUP7548U 

S0P75485 

SUP7-549U 


I  FORMAT (9X,52HHHEN  HOB  EXCEEDS  50  THE  HAP  PRINTING  IS  DISCONTINUED  SUP/5500. 
1/ 73 MO ,  48H  CALCULATIONS  PROCEED  IN  NORMAL  MANNER  )  SUP7551U 

GO  10  1000  .  SUP/5520 

1  RRItE (6,51 )F  Rt  Q( IFR) , NBOX, EKR  , EH  SUP75530 

WRITF(6,21)EM;EKB  SUP75540 

1  FORMATdH  ,28X,  48MPLANAR  VELOCITY  POTENTIAL  INFLUENCE  COEFFICIENTSSUP75550 
3  /1H0,30X,ll!HriACH  NO,  =,F8.5,10X,  6HKBAR  =, P9 . 5/1H0 » 13X, 1HI , 5X» 1H JSUP/5560 

2  10X,  5HNTJBAP  *  5X,  51!HUBaR»11X.,  1-4HKEAL  VP  I C ( I »  J ) , 8X, 14H T  MAG  VP  1 C  (  I:,  J  JSUP7557U 


3-  /in  ) 

K  =  0 

DO  22  I  =  l.NttOX 
DO  22  J  =  1,1 
fliRUU  «  I  -  3  . 

RARMU  =  J  -  1 
K  =  K'*1 

?  WRITF(6,23)  J , J, 8 ARNU, BARHU,VP 1 C (K ) 

5  FORMAT  (9X,  2  J  6 , 4'X,  2Fl  0 .3 ,2X,  1P2E22 .5  ) 

go  in  moo 

1  WR.I  1 E  ( 6»5l )  FRF'JMFR )', NBOX, EKR  , EH 
HR1TF(6,31>H0UF 


SUP75580 

SUP75590 

SUP/56IJU 

SUP/561U 

SuR;7562li 

SUP7b63.U; 

SUP/5640 

SUP/5650- 

SUP  75660 

SUP  75  67 1)5 

SUP756BU 

SUP75690 


I  FORMAT, (1H  ,23  X ,  69HIJPPER  VELOCITY  POTENTIALS  AND  SOURCE  SJRENGYHSSUP75700 

1  FOR  MODE  NO. 13/1  MO , 9X, 3 HN, 6X, 1 HH, 5X, 2HNB, 7 X, 1 OHR  PH K N, M 1, 7X, 10H I SUP7571U 

2  PH l  ON,  M ) iO 9HN  SS (  N, H ) , 8X, 9H I!  SS(N,H  )/1H  )  SUP75720 

GO  TO  1000  SUP/5730 

WRI TE ( 6,41 )  MP ,  Ml>,  '«B,  PH  I  ,SS( MB  )  SUP75740 

I  F  OR  HA  T  ( 4Xg,3 1 7  >  1  P2F;1 7.6 » 3X,  i;P2E17  •  51  SUP/5750 

GO  10  ,1  00  0  SUP‘7-5760' 

1  IF0MS.F0.3)  GO  fO  100  SUP/-5  770 

»  WRIT  F(6, 51  )F  REQGIFP )  i  HBOX,  EKR  ,,EM  SUP7578li' 

!  FORMAT  (1H1,o3X,44HM1SS!LE  SUPER’SONIC  AIRLOADS  PROGRAM  (CONT-D)//  SUP/5790 
X1H0  SUP7580  0 

1  8X,  27H0SC  VL  L  AT*1R  Y  FREQUENCY  (CPS  ) ,  F12 . 5 , 14  X,  1 2, 25H  BOXES  IN  CH0RDSUP7581U 

2  D  I  RFC  T I  ON  /HiO  ,  8 X,  30HREDUC ED  FREQUENCY  (SEMI  CHORD ) , F9.5, 14X, ,  SUP75820 

3  23HFRFE  STREAH  MACh  NUMBER , F9. 3 , /1H  .)  *  SUP75830 

HRITE(6,101)  SURF ( 1 »  NS )  SUP/5840 

I  FORMAT  <  2aX,45H  I  wPiJT  MODE  SHAPE  POLYNOMIAL  COEFFICIENTS  FOR  ,A6/  SUP/5850 

3  2?X,62HREn-RENCEn  10  THE  SURFACE  LEADING  EDGE-CENTERLINE  I  Nc?£RSECSUP75860 
2TION  7/2 X',,4 HMOflc » 2 0 X »  7HC0EPFS.  )  SUP/5870 

DO  102  I=l,MMODfc  S  U  P  7  5  8,8  0 

NTM=NT(I,NS)  SUP7589U 

?  WRITE (6, 10,3)  l ,  ( Cu(  I ,  o,  NS  ) ,  J=1 ,  NTH )  SUP/590U 

1  FORMAT(lHOi  I  4 , 4  X,  1P7E13 .4/ (;9X,  1  P7E13 . 4  )  )  SUP75910 

1  WRITE ( 6,'5l  )F  REO d  FR  ) ,  NBOX, EKR  ,EM  SUP/5920 

WRI1E(6,61)(S0RF( l, NS), 1=1,2)  SUP/593II 

I  FORMAT ( 1H  ,35X,23HGENERALIZED  FORCES  FOR  , 2 A6/1H0 , 6X , 4HDEF L, 3X ,  SUP7594U 

1  4  H  l  0  A  D.  >  1 0  X  >  9  H  R  E  A  L  PART , 10 X , 9H I  MAG  P  ART ,  1 0  X  ,  9HABS  VALUE,10X,  SUP/5950 

2  33HPHASE  AA’GLE  //)  SIJP/596U 

GO  10  3000  SUP/5970 

V  WR I  I E ( 6 , 71 ) 0 , K , 9R , Q I , QAB , QAN  SUP/5980 

IF  ( NSIJRF  .fcO.  11  Gli  TO  632 

IF  (NS.NE.3)  GO  TO  3000  #¥#*«*** 

1  KKK=2#NM0UE 


, 2 A6/TH0  » 6X , 4HDEF L,3X, 
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r 

[ 

NNN=2*K 

NNNN=NNN-1 

CARDS(J,NNNN)=QP 

******** 

***•***«, 

: 

i 

i 

r 

r  6969 

CARDS(.J*NNN)=UI 

IF  (,J.NE«NHO0E)  GO  TO  1000 

IF  (K.NE.NHODE)  GO  TO  1000 

PUNCH  6969,  {(C«RuS(II,JJ),JJ=l,KKK),  11=1, NMODE) 

FORHAT (6E12.5 )  5 

******** 

«•••«*** 

******** 

i 

71 

FORMAT (1H0, 1 9, 1 7 , 2* , lP3El9i 5, OFF  16 .3, 4H  DEG ) 

SUP7599U 

* 

i 

I 

jv 

1  91 

i  0  U  <1 

L 

FORMA  r(lHl,30X>44HHlSSILF.  SUPERSONIC  AIRLOADS  PROGRAM  (CQMT-D)//) 
RETURN 

SUP76000 

SUP76010 
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